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ABSTRACT 


Ship  Synthesis  Model  for  Coast  Guard  Cutters 

by:   Michael  Jay  Goodwin 

Submitted  to  the  Department  of  Ocean  Engineering  on 

May  9»  1975*   in  partial  fulfillment  of  the  requirements 

for  the  degree  of  Ocean  Engineer  and  the  degree  of 

Master  of  Science  in  N-/al  Architecture  and  Marine  Engineering 

In  order  to  determine  the  optimum  design  of  a  Coast 
Guard  cutter,  it  is  necessary  to  perform  a  large  number 
of  feasibility  studies.   This  is  an  impossible  task  if 
done  using  hand  calculations.   A  mechanized  version  of 
these  hand  calculations  is  developed  which  will  permit 
the  calculations  to  bo  performed  by  a  digital  computer. 
The  final  program  listing  is  given  together  with  the 
development  of  the  empirical  relations  used. 

The  program  is  applicable  to  search  and  rescue  and 
patrol  type  cutters  with  lengths  between  150  and  400  feet. 
Twin  screw  propulsion  is  assumed. 

A  check  on  the  accuracy  of  the  program  is  made  using 
three  U.S.  Coast  Guard  cutters  now  in  service  and  the 
U.S.  Navy's  Patrol  Frigate  design.   The  program  is  found  to 
be  much  more  accurate  than  existing  programs  when  applied 
to  Coast  Guard  vessels.   The  program  predicts  the  full  load 
displacement  within  four  percent  for  the  Coast  Guard  designs 
tried. 
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CHAPTER  I 
INTRODUCTION 

Although  digital  conputers  have  been  used  by  engineers 
for  more  than  twenty  years,  their  potential  as  design 
aids  is  only  beginning  to  be  realized.   Naval  architects 
made  use  of  computers  soon  after  they  became  available 
for  solving  specific  algorithms  such  as  displacement  sheet 
calculations  and  speed-power  estimates.   Such  algorithms 
ar-3  beneficial  in  the  preliminary  and  later  design  stages. 
However,  with  the  growth  of  systems  engineering  in  ship 
design,  m.ore  emphasis  is  being  placed  on  the  conceptual  or 
earliest  phase  of  a  design.   At  this  stage  of  the  design, 
th^  computer  can  be  used  more  effectively  than  at  any  other 
stage . 

Because  of  the  complexity  of  the  naval  architectural 
calculations  involved  in  later  design  stages,  the  principle 
characteristics  of  a  ship  must  be  chosen  early  in  the  design. 
The  output  of  the  conceptual  design  phase  is  this  set 
of  characteristics.   Without  the  use  of  a  computer,  the 
best  that  can  be  done  is  to  develop  a  feasible  set  of 
characteristics  when  a  optimum  set  is  desired.   Systems 
engineering  is  placing  more  emphasis  on  obtaining  this 
optimum. 

To  choose  an  optimum  or  near  optimum  set  of  character- 
istics, a  great  many  feasibility  studies  are  required.   It 
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8 
is  impossible  to  do  more  than  a  few  feasibility  studies 

without  the  aid  of  computerized  models  because  of  the 

time  required.   At  this  early  stage,  the  ship  must  be 

designed  as  a  whole.   This  precludes  dividing  the  work 

into  a  number  of  parallel  paths  with  several  people 

working  simultaneously  on  a  single  feasibility  study. 

jc)ach  study  is  essentially  a  one  man  job. 

V/ith  hand  calculations  a  good  engineer  can  perform 
about  one  feasibility  study  per  day.   As  a  result,  great 
reliance  must  be  placed  on  the  designers  intuition.   In 
the  past  this  has  led  to  series  of  designs  which  may  or 
may  not  have  been  optimum.  £ach  new  design  was  a  slight 
perturbation  of  a  previous  successful  design.   While  this 
procedure  produces  satisfactory  designs,  bad  features 
tend  to  be  duplicated  and  few  new  features  are  added. 
The  engineer  was  not  really  at  fault  because  no  tool 
was  available  to  perform  all  the  calculations  necessary 
to  obtain  an  optimum  design. 

Computers  permit  the  mechanization  of  these  hand 
calculations.   The  software  required  will  be  referred  to 
here  as  a  ship  synthesis  model.   The  model  aids  the  designer 
by  freeing  him  from  the  routine  calculations  which  take 
up  most  of  his  time  during  the  early  phases  of  design. 
The  algorithms  referred  to  previously  also  free  the 
designer  from  tedious  calculations. 

"hat  makes  a  ship  synthesis  model  different  is  that 
it  goes  beyond  merely  reducing  the  time  and  effort  required 


to  do  calculations  previously  done  by  hand.   Using  the 
model,  hundreds  or  even  thousands  of  feasibility  studies 
can  be  made.   vVhile  this  still  does  not  guarantee  that 
the  optimum  solution  is  obtained,  the  engineer  will  have 
much  more  confidence  in  his  design  than  he  would  if  only 
a  few  feasibility  studies  had  been  made.   The  designer 
is  free  to  concentrate  on   the  more  important  aspects  of 
the  early  stages  of  ship  design  such  as  systems  level 
tradeoff  decisions.   These  include  hull  form,  speed, 
endurance  and  payload  options  among  others.   These  aspects 
were  seldom  considered  in  detail  before  the  advent  of 
ship  synthesis  models. 

Equally  important  is  the  fact  that  computers  permit 
a  consistency  among  feasibility  studies.   With  several 
people  doing  independent  feasibility  studies  or  even 
with  a  single  person  doing  them  all,  it  is  difficult  to 
insure  that  the  same  assumptions  will  be  made  for  each 
study  and  that  no  bias  is  applied.  The  model  guarantees 
consistency.   And  consistency  is  important  because  without 
it  two  feasibility  studies  can  not  be  compared  and  an 
optimum  solution  is  difficult  to  select. 

Several  models  have  been  available  for  a  number  of 
years  and  have  proven  their  value.   The  U.S.  Navy's 
destroyer  model,  DD0?v  '  ,  on  which  the  program  developed 
in  this  thesis  is  patterned,  was  started  in  the  1960's 


*  Numbers  in  parenthesis  refer  to  References  at  the 
end  of  this  thesis 
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and  updating  work  has  been  continuing  since  then.   This 

war  the  first  true  Hsys terns"  model.   Models  also 
have  been  developed  which  attempt  to  model  all  surface 
ship  types  with  a  single  program.   For  accurate  results, 
each  different  ship  type  really  requires  its  own  model. 

In  addition  to  the  destroyer  model,  the  U.S.  Navy 
has  developed  individual  synthesis  models  for  aircraft 
carriers,  underway  replenishment  ships,  submarines, 
hydrofoils,  and  several  other  ship  types.   The  U.S.  Coast 
Gu-ird  presently  has  no  model  which  gives  good  results  for 
cutters.   The  destroyer  model,  DD07»  does  not  apply  to 
ships  as  short  as  150  feet  in  length  which  is  one  Coast 
Guard  requirement.   Some  all  inclusive  models  have  been 
tried  but  with  poor  success,  with  predictions  of  displace- 
ment as  much  as  30  percent  off  when  data  on  existing 
designs  is  inputed. 

The  program,  which  will  be  referred  to  as  the  cutter 
model,  developed  in  this  thesis  is  an  attempt  to  satisfy 
the  Coast  Guard's  need.   The  cutter  model  applies  to 
search  and  rescue  and  patrol  type  cutters  with  lengths 
between  150  and  ^00  feet.   It  is  not  applicable  to  buoy 
tenders,  tugs  or  ice  breakers.   An  effort  has  been  started 
by  the  Coast  Guard  to  develop  a  model  for  buoy  tenders 
and  tugs  but  this  work  is  not  yet  complete. 

Zven   in  the  early  stages  of  design  where  the  synthesis 
models  are  of  use,  a  great  many  calculations  are  required 
to  prove  the  feasibility  of  a  ship  design.   This  has  led 
designers  to  employ  a  type  of  regression  analysis  where 
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data  from  existing  ships  is  scaled  up  or  down  in  order 
to  estimate  values  for  a  new  design.   This  estimation 
procedure  is  employed  in  ship  synthesis  models  also. 
More  detailed  calculations  could  be  made  but  at  the  cost 
of  greatly  increased  computer  time  since  a  large  number 
of  feasibility  studies  will  be  run. 

The  shortcut  estimating  techniques  employed  place 
more  emphasis  on  the  relative  accuracy  among  studies 
rather  than  on  the  absolute  accuracy  of  the  result. 
Because  of  this,  all  estimating  relationships  in  the 
model  mu3t  show  the  dependence  of  the  result  on  all 
significant  input  variables.   At  present,  this  is  an  im- 
possible task  because  not  enough  is  known  about  this 
dependence.   As  models  become  more  widely  used  much 
development  work  v/ill  have  to  be  done  in  this  area.   An 
estimate  of  the  variation  with  the  most  significant 
input  variable  is  the  best  that  can  be  done  at  present. 

Because  the  estimating  relationships  are  rather 
crude,  estimates  are  only  made  at  the  ship  system  level. 
Ship  system  level  information  consists  of  such  elements 
as  total  outfit  weight,  total  volume  required  for  machinery 
and  required  electrical  load  for  ventilation.   Although 
total  arrangements  area  is  computed,  no  attempt  is  made 
to  perform  an  arrangement  of  spaces  or  to  compute  trim. 
The  engineer  usually  has  sufficient  flexibility  in  these 
areas  to  arrive  at  a  feasible  design  if  sufficient  space 
is  available. 


12 

One  problem  in  developing  a  model  for  Coast  Guard 

cutters  is  that  there  are  few  designs  available  from 
which  to  obtain  data.   This  causes  difficulty  in  deter- 
mining rational  scaling  laws.   In  many  cases  the  use  of 
naval  architectural  principles   indicates  the  correct 
scaling  procedure.   In  other  cases  all  that  can  be  done 
is  to  search  for  the  most  logical  correlation. 

The  design  produced  by  a  synthesis  model  should 
be  tailored  to  meet  the  conditions  required  of  a  feasible 
solution.  First  there  must  be  a  balance  between  weight 
and  displacement.   Second,  internal  space  available  must 
equal  internal  space  required.   This  requirement  is  becoming 
very  important  as  most  Coast  Guard  vessels  are  volume 
rather  than  weight  limited.   Third,  the  energy  available 
must  equal  the  energy  required.   This  includes  both 
propulsion  and  electrical  power  requirements.  Finally, 
the  distribution  of  weight  and  volume  must  be  such  as  to 
satisfy  design  criteria  for  transverse  stability,  girder 
strength  and  sea  keeping.   There  are  also  other  require- 
ments implicit  in  the  solution  such  as  hull  type,  material 
of  construction  and  space  requirements. 

In  jeneral,  the  inputs  required  for  a  ship  synthesis 
model  are  a  description  of  the  "payload"  to  be  carried, 
a  statement  of  vehicle  mobility  performance  and  a  selection 
of  options  to  be  used.   The  payload  description  must  be 
more  detailed  than  a  simple  weight  breakdown.   The  relative 
vertical  centers  of  gravity,  area  requirements,  electrical 
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pov/er  requirements  and  cost  must  also  be  given.   Payload 

for  the  cutter  model  also  includes  people,  and  thus  the 
size  of  the  crew  is  an  input.   The  description  of  vehicle 
mobility  performance  includes  both  maximum  calm  water 
speed  and  endurance  speed  and  range  requirements.   Options 
include  length,  Cp,  Cx,  as  well  as  machinery  type  and 
margins.   There  is  also  an  option  in  the  cutter  model 
to  input  the  maximum  sustained  horsepower,  in  which  case 
the  program  will  calculate  the  maximum  sustained  speed. 
A  full  description  of  the  program  inputs  for  the  cutter 
model  is  contained  in  Chapter  III  (Main  Program  Description) 
and  in  Appendix  A. 

The  output  generated  by  the  model  must  be  sufficient 
to  demonstrate  the  feasibility  of  the  design.   V/herever 
possible,  the  output  should  permit  the  designer  to  verify 
the  feasibility  v/ithout  doing  additional  calculations. 
On  a  few  of  the  studies  run  and  in  particular  on  the  design 
selected,  the  design  team  will  have  to  lay  out  a  rough 
arrangements  drawing  to  verify  that  sufficient  space  is 
-available.   Sufficient  information  should  be  given  in  the 
output  to  permit  this  to  be  done  easily. 

The  output  generated  should  also  allow  the  designer 
to  rapidly  determine  the  design  that  best  satisfies  his 
optimization  criterion.   This  criterion  usually  takes  the 
form  of  a  total  acquisition  cost  or  life  cycle  cost. 


*!2he   definition  of  coefficients  and  dimensions  used 
in  this  thesis  may  be  found  in  Reference  (3). 
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Predictions  of  acquisition  cost  are  often  included  in  a 

ship  synthesis  model. 

finally,  the  output  should  provide  a  basis  for  beginning 
the  preliminary  design  phase.   Generally  the  preliminary 
phase  is  handled  by  a  number  of  design  teams  working  in 
parallel.   The  output  produced  must  provide  information 
to  initialize  each  of  These  parallel  developments. 

Sample  output  listings  for  the  cutter  model  are 
given  in  Appendix  3. 

It  is  important  to  emphasize  that  the  cutter  model 
produces  a  feasible  solution  only  and  not  an  optimum 
solution.   An  optimum  solution  implies  that  the  best  of 
a  very  large  number  of  feasible  ships  has  been  chosen. 
It  is  felt  that  the  design  process  works  more  efficiently 
if  the  designer  acts  as  a  guide  for  the  computer  rather 
than  having  the  computer  do  the  entire  optimization  on  a 
single  run.   In  this  way,  many  unnecessary  designs  can 
avoided  in  areas  where  designs  are  less  attractive  than 
others  already  tried.   Such  a  process  takes  more  real 
tire  but  less  computer  time.   It  also  gives  the  designer 
a  much  better  feel  for  the  sensitivity  of  the  optimum 
to  changes  in  characteristics. 

The  remainder  of  this  thesis  is  arranged  as  follows: 

Chapter  II  gives  an  overview  of  the  cutter  model 
and  explains  the  major  options. 

Chapter  III  goes  into  detail  on  each  of  the  sub- 
routines in  the  model  and  explains  the  derivation  of  all 
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of  tho  empirical  formulas  used.  This  chapter  may  be 

skipped  by  anyone  not  interested  in  the  details  of  the 
program. 

Chapter  IV  gives  an  evaluation  of  the  model  including 
test  runs  made  using  existing  ships. 

Chapter  V  contains  recommendations  and  conclusions. 

Appendices  A  and  3  contain  input  and  output  formats 
respectively. 

Appendix  C  is  a  listing  of  the  program. 

The  remaining  appendices  give  background  information 
used  in  this  work. 

It  is  recommended  that  at  least  Chapters  I,  II,  IV,  V 
and  Appendix  A  be  read  before  attempting  to  use  the  program. 
Chapter  III  is  aimed  primarily  at  the  person  who  wishes 
to  modify  the  program.   However,  the  important  aspects 
of  each  subroutine  are  covered  in  the  introduction  to 
each  of  the  subsections  of  Chapter  III. 


CHAPTER  II 
GENERAL  PROGRAM  DESCRIPTION 

This  chapter  gives  an  overall  description  of  the 
solution  method  used  by  the  cutter  model  but  does  not 
discuss  the  detailed  calculations.   The  detailed  calculations 
are  explained  in  Chapter  III.   This  chapter  should  also 
serve  as  a  guide  to  understanding  how  each  subroutine 
discussed  in  Chapter  III  fits  into  the  total  program. 

The  organization  of  the  computer  program  can  be 
viewed  in  two  different  ways.   The  first  is  the  program 
control  organization  which  is  shown  in  Figure  1,   Sub- 
routine XclCUTS  controls  the  entire  solution  of  each 
feasibility  study  case.   Control  is  returned  to  the  main 
program  only  when  data  for  a  new  case  is  required. 
Subroutine  XiCUTS  calls  other  subroutines  as  needed  and 
gradually  fills  all  the  data  storage  locations  needed 
for  output.   This  subroutine  also  controls  the  iteration 
loops  for  weight  and  vertical  center  of  gravity. 

The  logical  organization  of  the  program  is  shown 
schematically  in  Figure  2.   This  is  an  easier  organization 
scheme  to  follow  and  shows  more  clearly  how  the  program 
solves  each  input  case.   It  will  be  noted  that  subroutine 
XECUTE  has  not  been  shown  in  Figure  2.   As  explained 
above,  X^CUTS  merely  guides  the  computer  around  the 
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logical  organization  chart.  Error  messages  are  also 
printed  by  subroutine  XECUT2  if  any  constraint  of  the 
model  is  violated. 

The  program  begins  with  data  input  in  the  MAIN 
program.   The  primary  purpose  of  the  MAIN  program  is  to 
read  data  cards  and  all  input  is  done  in  this  routine, 
borne  bookkeeping  is  also  done  in  MAIN  such  as  summing 
payload  weights  and  costs.   The  final  purpose  of  this 
routine  is  to  estimate  initial  values  of  vertical  center 
of  gravity,  ?CGj  displacement,  and  cubic  number. 

Once  all  data  for  one  feasibility  study  has  been 
read  by  MAIN  and  the  initial  values  of  KG,  displacement  and 
cubic  number  have  been  determined,  control  is  transferred 
to  subroutine  X2CUTE.   This  subroutine  directs  the  execution 
of  the  program  until  a  solution  is  obtained  or  until  • 
a  constraint  is  violated,   "tfhen  either  of  these  occur 
control  reverts  to  the  MAIN  program  where  data  for  the 
next  feasibility  study  is  read.   If  there  is  no  more 
data  the  program  stops. 

The  first  subroutine  called  by  XECUTiS  is  an  under- 
water dimensions  subroutine,  UWDIM.   Subroutine  UWDIM 
uses  the  input  values  of  length,  free  surface  correction, 
CD,  Cx  and  Gi'./3eam  as  well  as  the  estimated  KG  and  dis- 
placement to  calculate  the  required  beam,  draft  and 
waterplane  coefficient,  CWD.   To  the  level  of  accuracy 
of  this  program  these  completely  determine  the  under- 
water shaoe  of  the  hull. 


20 

Two  options  are  available  with  regard  to  speed  and 

power  calculations.   Under  the  first  option  the  maximum 
sustained  speed  of  the  vessel  is  input.   If  this  option 
is  specified,  XiCU'IE  calls  subroutine  HPCA1C  which, 
together  with  subroutine  CRVAL,  estimates  the  horsepower 
required  to  maintain  the  maximum  sustained  speed. 

Option  two  allows  the  user  to  specify  the  maximum 
sustained  horsepower.   In  this  case  the  program  will 
compute  the  maximum  sustained  speed.   Subroutine  HPCALC 
has  been  written  to  take  speed  as  an  input  and  return 
a  value  of  horsepower.   Therefore,  an  iterative  procedure 
is  required  if  this  second  option  is  used.   Control  of 
this  iteration  takes  place  in  subroutine  XSCUTE. 

Regardless  of  which  of  the  above  options  is  used 
the  endurance  speed  is  input.   Subroutine  HPCALC  is 
called  to  determine  the  power  required  at  endurance  speed. 

A  Taylor  Standard  Series  horsepower  estimate  is 
used  in  subroutine  HPCALC.   Stored  values  of  Cr  are 
used.   The  purpose  of  subroutine  CRVAL  is  to  choose  the 
correct  Cr  value  from  a  stored  array.   Reference  (7) 
explains  this  method  for  estimating  horsepower. 

The  electrical  power  routine,  SPLANT,  is  the  next 
one  called.   This  routine  determines  the  maximum  electrical 
load  at  sea,  the  average  electrical  load  and  sizes  the 
generators.  The  average  load  is  needed  for  estimating 
generator  fuel  consumption.   The  input  electrical  loads 
along  with  crew  size  and  cubic  number  are  used  in  this 
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routine. 

The  weight  of  fuel,  lubricating  oil  and  potable 
water  is  next  calculated  in  subroutine  MACHLQ.   The 
endurance  range  given  as  input  together  with  crew  size 
and  the  endurance  power  calculated  earlier  are  used  by 
this  routine.   Lub  oil  requirements  are  either  input, 
in  the  case  where  horsepower  is  specified,  or  are  estimated 
based  on  specified  machinery  type.   Potable  water  require- 
ments depend  entirely  on  the  crew  size.   If  make-up  feed 
water  is  required  it  should  be  specified  as  a  cargo  load 
item.   The  specific  fuel  consumption  is  either  input,  in 
the  case  where  horsepower  is  specified,  or  is  estimated 
based  on  machinery  type  and  horsepower. 

Subroutine  MB3IZ2  is  next  called  if  the  sustained 
speed  is  specified.   This  routine  estimates  the  machinery 
box  length  and  minimum  depth.   If  the  horsepower  input 
option  is   specified  these  dimensions  are  input. 

The  next  subroutine  is  one  of  the  most  important 
in  the  program.   Subroutine  VOLUME  determines  both  the 
sheer  line  and  the  deckhouse  volume.   These  complete 
the  definition  of  internal  volume  and  arrangements  area 
since  the  underwater  volume  has  been  determined  in  sub- 
routine UV/DIw.  Because  most  ships  for  which  this  program 
applies  are  volume  rather  than  weight  limited,  it  is 
important  to  insure  that  there  is  sufficient  space  for 
all  ship  functions  and  equipment  which  must  be  installed. 

The  volume  is  divided  into  three  categories.   First, 
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there  is  the  engine  room.   The  dimensions  for  the  engine 

room  have  been  estimated  earlier.   Second  is  tankage 
volume  and  finally  arrangements  volume.   Arrangements 
volume  is  more  useful  if  converted  to  an  arrangements 
area. 

As  a  first  step  in  determining  the  space  available 
in  each  of  these  categories,  the  sheer  line  must  be  deter- 
mined.  Subroutine  VOLUi-IZ  calls  subroutine  SHEER  to 
determine  an  acceptable  sheer  line.   The  sheer  line  must 
meet  several  criteria  including  adequate  freeboard  forward 
and  aft,  adequate  depth- amidships  for  the  engine  room 
and  for  adequate  strength,  and  an  acceptable  sheer  line 
curvature . 

V/ith  the  sheer  line  fixed,  the  total  volume  of  the 
hull  is  estimated.   The  volume  of  the  engine  room  is 
also  estimated  and  then  subtracted  from  the  total  volume. 
Next,  an  estimate  is  made  of  the  fraction  of  the  remaining 
volume  vvh^ch  can  be  used  only  for  tankage  because  it  is 
unsuitable  for  use  as  arrangements  volume. 

A  required  tankage  volume  is  then  estimated  based  on 
the  weights  of  fuel,  lub  oil,  potable  water  and  aircraft 
fuel.   This  required  tankage  volume  is  compared  to  the 
volume  which  must  be  used  for  tankage.   The  larger  of  the 
two  is  subtracted  from  the  hull  volume  that  remains.   The 
remaining  hull  area  is  available  for  use  as  arrangements 
volume.   This  volume  is  converted  into  an  arrangements  area. 

A  required  arrangements  area  is  computed  based  on 
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empirical  data  from  past  designs  and  on  input  area  values. 

By   adjusting  the  size  of  the  superstructure  or  by  adding 
a  raised  deck,  the  available  arrangements  area  can  usually 
be  made  to  equal  the  required  arrangements  area.   If 
too  large  a  deckhouse  is  required  or  if  a  raised  deck 
longer  than  the  length  of  the  ship  is  needed,  an  error 
message  is  generated  and  a  new  input  case  is  tried. 

Once  the  volume  of  the  hull  has  been  determined, 
the  weight  of  all  the  light  ship  weight  groups  and  of  the 
load  items  can  be  calculated.   This  is  done  in  subroutine 
WEIGHT.   Although  many  of  the  items  calculated  by  this 
routine  are  later  printed  as  output,  the  variable  of 
immediate  importance  is  the  full  load  weight,  WFULLD. 
This  value  is  compared  with  the  initial  guess  of  displace- 
ment.  If  there  is  a  difference  of  greater  than  2  tons 
between  the  two,  a  new  estimate  of  displacement  is  made 
and  the  program  returns  to  subroutine  UWDIi-i  for  a  new 
estimate  of  dimensions.   Once  a  balance  between  WFULLD 
arid  displacement  is  achieved  the  program  proceeds  to  the 
-  estimation  of  vertical  center  of  gravity. 

Subroutine  VCG  makes  this  estimate.   Most  centers 
are  taken  as  a  fraction  of  the  average  depth  of  the  hull 
or  of  the  amidships  depth.   The  estimated  vcg  is  compared 
to  the  initial  guess  and  a  new  KG  estimate  is  made  if 
there  is  no  agreement  within  a  tenth  of  a  foot.   Once 
a  balance  is  obtained,  the  program  proceeds  to  a  cost 
estimate . 


The  lead  ship  cost  is  estimated  using  a  Coast  Guard 
cost  estimating  procedure  known  as  "Flanagan's  r.ethod">°' 
This  procedure  has  been  programmed  as  subroutine  COST, 

An  estimate  of  the  speed  the  vessel  will  be  able 
to  sustain  in  the  North  Atlantic  Ocean  is  made  in  sub- 
routine J£A3?D.   This  routine  is  based  on  work  done  at 
the  i.aval  ^hip  Engineering  Center  as  reported  in  References 
{k)   and  (5). 

Finally i  subroutine  OUTPUT  is  called.   The  purpose 
of  this  routine  is  to  print  the  values  determined  by  the 
other  routines.  The  input  data  is  also  printed  so  that 
a  compact  record  of  the  feasibility  study  is  produced. 
After  the  output  is  printed,  the  program  returns  to  the 
nAIN  program  to  read  more  data  and  start  a  new  case. 

This  concludes  the  general  description  of  the  program. 
A  much  more  detailed  description  is  given  in  Chapter  III. 
The  reader  v/ho  is  not  interested  in  the  details  of  the 
program  can  skip  Chapter  III  with  no  loss  of  continuity. 


CHAPTER  III 
D2TAIL2D  DESCRIPTION  OF  THE  PROGRAM 

3.1  Introduction 

This  chapter  is  composed  of  sixteen  sections  including 
this  introduction.   The  other  fifteen  sections  each 
explain  in  detail  one  subroutine  of  the  cutter  model. 
For  the  location  of  these  subroutines  in  the  overall 
program  refer  to  Figure  2,  Logical  Program  Organization. 

An  attempt  has  been  made  to  keep  each  section  in  the 
same  format.   A  brief  description  of  the  routine  will 
be  given  first  followed  by  a  list  of  required  inputs. 
Statements  will  then  be  explained  one  by  one.   A  list 
of  the  outputs  produced  by  the  subroutine  will  be  given 
next.  Finally,  a  nomenclature  list  will  be  given.   The 
nomenclature  list  applies  only  to  the  subroutine  being 
considered.   In  some  cases  the  same  variable  will  have 
different  names  in  different  subroutines. 

In  both  the  input  and  output  listings  a  code  is 
employed  to  indicate  whether  the  variable  is  transferred 
by  a  labeled  common  statement,  LC(AA),  or  as  a  subroutine 
call  dummy  argument,  3CDA.   The  letters  in  parenthesis 
after  the  labeled  common  code  indicate  which  labeled 
common  was  used  to  transfer  the  argument. 
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3.2     MAIN   Program 

3.2.1  introduction 

The  MAIN  subroutine  performs  several  functions. 
First  and  most  important  is  the  reading  of  all  input 
data  used  by  the  model.  Each  feasibility  study  case 
starts  with  the  reading  of  one  or  more  data  cards  by 
the  MAIN  program.   The  actual  execution  of  each  case 
is  controlled  by  subroutine  XECUTS.   This  routine  is 
discussed  in  Section  3»3»   However,  once  execution  on  a 
case  is  complete,  control  is  returned  to  the  MAIN  sub- 
routine.  A  check  is  then  made  to  see  if  there  is  more 
input  data.   If  not  the  MAIN  subroutine  directs  the 
computer  to  stop. 

Several  bookkeeping  tasks  are  performed  in  this 
subroutine  also.   These  include  such  jobs  as  adding 
together  the  deck  area  and  electrical  load  requirements 
of  the  individual  payload  items.   The  final  function 
performed  is  the  setting  of  initial  values  for  displace- 
ment, KG  and  cubic  number.   These  values  are  used  for  the 
-  first  iteration  and  are  then  replaced  by  better  estimates. 

3.2.2  Inputs 

Block  Data  subroutines  have  not  been  used  in  this 
model.   As  a  result  there  are  no  variables  with  assigned 
values  at  the  beginning  of  the  MAIN  program.   All  of  the 
variables  required  are  read  from  data  cards  or  assigned 
values  in  the  routine. 

A  flow  chart  for  this  routine  is  given  in  Figures 
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3a  and  3b.   A  listing  of  the  program  is  given  in  Appendix 

C.   The  nomenclature  list  at  the  end  of  this  section 
gives  a  description  of  each  variable  name  used, 
3.2.3  Statement  Descriptions 

A  line  by  line  discussion  of  the  program  statements 
follows.   The  REAL  and  COMMON  statements  have  not  been 
listed  but  may  be  founc*  in  Appendix  C. 

REAL  Statements  -  An  explicit  type  definition  has 
been  employed  so  that  the  variable  names  will  be  more 
descriptive  of  their  actual  meanings. 

COMMON  Statements  -  Because  a  great  many  variables 
are  required  to  define  each  case,  COMMON  statements  have 
been  used  to  transfer  values  between  subroutines.   This 
also  saves  computer  storage  space.   All  COMMON  statements 
used  are  labeled  common. 

Cr  array  input: 

READ(tf,120)  CR1 

120  FORMAT (16F5.0; 
KEAD(8,121)  CR2 

121  F0RMAT(9F5.0) 
READ(tf,122)  CR3 

122  FORMAT (6F^.O) 

The  three  arrays  which  contain  the  Taylor's  Standard 
Series  Cr  values  are  read  only  once  at  the  beginning  of 
the  program. 

INDEX  input: 

2  READ( 8,124)  INDEX 
124  FORMAT(Il) 

The  variable  named  INDEX  has  two  functions.   First, 
as  described  below,  it  tells  the  program  when  to  stop. 


MAIN      SUBROUTINE 
FLOW     CHART 
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It  is  also  used  to  discern  between  full  hull  ships  and 

fine  hull  ships  for  the  initial  guess  of  displacement. 

Stop  execution: 

IF(INDEX.EQ.O)  STOP 

This  is  the  only  place  where  the  program  is  directed 
to  stop  execution.   However,  if  a  data  card  is  left  out, 
the  program  will  stop  when  there  is  no  more  data  cards 
to  read. 

Heading  input: 

RBAD(8,123)(HEAD(I),I=li20) 

123  F0RMAT(20A<O 

Variable  HEAD  contains  any  alphanumeric  information. 
It  will  be  printed  without  change  at  the  beginning  of  the 
output  for  all  test  cases  using  the  data  which  follows. 

NOARM  input: 

READ(S,1U0;  NOARM 

100  F0RMAT(I2) 

NOARM  is  the  total  number  of  armament,  aircraft 
and  cargo  payload  data  cards  which  follow.   A  maximum 
value  of  20  is  allowed. 

Armament,  aircraft  and  cargo  input: 

READ(«,101)((TITLE(J,I),J=1,5),CARGOW(I),CARGOC(I), 

GRyWT(I),GRyCG(l),AKOWT(I),AMOCG(Ij,ACWT(I), 
ACCG(I),ELLD(I),DHAR(I),HLAR(I),GR?CST(I),I= 
1, NOARM) 

101  F0RMAT(5A4,12F5.0) 

Up  to  20  armament,  aircraft  and  cargo  payload  descrip- 
tion cards  are  allowed.   The  data  may  be  lumped  together 
and  appear  on  only  one  card  or  it  may  be  broken  up  in  any 
way  desired.   However,  if  a  weight  is  given,  a  vertical 
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center  of  gravity  must  also  be  given  on  the  same  card 

immediately  following  the  weight.   The  first  20  spaces 

on  each  card  are  reserved  for  an  alphanumeric  description 

of  the  payload  items  on  the  card.   This  is  used  as  a  line 

title  in  the  output  listing.   Of  the  four  weight  categories 

on  the  card  only  the  group  7  weights  are  included  in  the 

light  ship  weight.   From  a  calculation  standpoint  it  is 

unimportant  whether  ammunition  or  aircraft  are  listed 

separately  or  considered  cargo.   However,  if  they  are 

listed  as  cargo  they  will  be  classified  incorrectly  in  the 

output  listing. 

All  vertical  center  of  gravity  values  should  be 
inputed  as  a  multiplier  of  the  hull  depth  at  amidships. 
This  number  is  easier  to  estimate  than  the  actual  center 
in  feet.   As  an  example,  if  a  payload  item  has  its  center 
10  feet  above  the  Daseline  and  the  depth  at  amidships 
is  20  feet  the  value  entered  is  0.5* 

The  deckhouse  area  should  include  only  the  area 
that  is  required  to  be  in  the  deckhouse.   If  the  area 
could  be  in  either  the  deckhouse  or  the  hull  it  should 
be  included  in  the  hull  area.   The  group  7  cost  should 
include  only  the  cost  of  materials  and  should  be  estimated 
for  the  time  when  the  ship  will  be  built.   Because  of 
format  limitations  the  value  input  should  be  the  true 
cost  divided  by  1000.   All  weights  are  in  tons,  all  areas 
in  square  feet,  all  electrical  loads  in  KW  and  all  costs 
in  dollars.   This  applies  to  all  the  input  items. 
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Summing  the  payload  data: 

WCARG0=0. 

WTG7=0. 

HfTAMMOsO, 

WTAC=0. 

EL0AD=0. 

BLOAD=0. 

AREADH=0. 

AREAHL=0. 

CSTG/=0. 

DO  10   I=1,N0ARM 

WCARG0=WC ARG0+CARG0W( I ) 

WTG7=WTG7+GR7WT( I ) 

WTAMMO=WTAMMO+AMOWT(IJ 

WTAC=WTAC+ACWT(I) 
BL0AD=BL0AI>4-ELLDCI  ; 
GR7CST(I ;=GR7CST(I )*10O0. 
CSTG7=CSTG7+GR7CST(I ) 
AREADH=ARi£ADH+DHAR  ( I ) 
10  AREAHL=AREAHL+HLAR(I) 

Only  the    totals   or   the   payload   items   entered  are 
needed  Tor  calculations.      These   totals  are   computed  in 
the   above   statements. 

NOELT  input* 

READ(b\102;   NOJiLT 
102  FORMAT(IZ) 

This  variable   serves   the   same   purpose   as  NOARM.    A 

maximum  value   of  20  is   allowed. 

Electronics   input « 

RBAD(»,103)((TITtE(J,I+20)tJ=l,5).GR4WT(I),GR4CG(I), 

ELD(I),DHA(I),HLA(I).GR4CST(I)fI=l,N0ELT; 
10}  FORMAT(3A4,6F5.0) 

The  comments  made  about  armament  input  also  apply 

here.   Group  4  weights  are  included  in  the  light  ship 

weight.   All  items  normally  in  group  4  should  be  included 

since  no  allowance  is  made  in  the  program  for  navigation 

equipment. 
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Summing  electronics   input  data: 

WTG4=0. 
CSTG4=0. 
DO  11    I=1,N0ELT 
WTG4=WTG4+GR4WT( I ) 
ELOAD=ELOALM-ELD(  I ) 
Gtt*CST(I)=GR4CST(I)*1000. 
CSTG4=CSTG4+GR4CST( I ) 
AREADH=AREADH+DHA( I ) 
11   AR£AHL=AREAKL+HLAU; 

The   comments   about  summing  armament  data  also  apply- 
here.      Note    that  the   area  in   the   hull   and  deckhouse    is 
a  combined   total   or  all   the   payload  inputs,    both  electronic 
and  armament. 

Storing  initial  costs: 

G4CST=CSTG4 
G7CST=CSTG7 

Because    the  value    or  variables  CSTG4  and  CSTG7  are 
modified  during  the   execution   of   the   program,    a  record 
is  required   of   their   initial  values   so   that  they  will 
be   available    if  another  test  case   using  the   same   payload 
is  run. 

Manning  input: 

READ(8,10^;   N00FF,N0CP0,N0CREW 
104  PORMATC3I5J 

N0FF=FL0AT(N00FF) 
NCP0=FL0AT(N0CP0j 
NENL=FLOAT(NOCREW) 

The  number  of   officers,    chief  petty  officers   and 
crew  members  are   input  and  then  converted  to  real  values 
so   that  they  can  be   used  in  non-integer  equations. 

Cost   indices   input: 

READ(tf,110;D0LHR,GlIND,G2IND,G3IND,G5IND,G6IND 
110  F0RMAT(6F5.0) 
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The  labor  rate  for  the  year  the  ship  will  be  built 

must  be  estimated  as  must  the  inflation  corrections  for 
material  costs.   These  corrections  are  for  the  period 
i'rom  19^9  to  the  year  the  ship  is  to  be  built.   Section 
3.14  gives  some  guidance  as  to  the  values  to  use.   The 
indices  can  also  include  special  material  costs  such  as 
making  the  hull  out  of  !IY-80  steel  instead  of  mild  steel. 
Military  mission  consumables  input  1 

READC « , 1 05 jENDDAY , WACFUL , CACFUL 

105  F0RMAT(3F5.0) 

The  number  of  endurance  days  for  dry  provisions  should 
be  given.   All  other  endurance  days  are  taken  as  a  suitable 
fraction  of  this  number.   The  weight  and  center  of  gravity 
of  aircraft  fuel  must  also  be  given. 

HPMS=0. 

This  is  done  so  that  there  is  no  difficulty  in 
calling  subroutine  XECUTE. 

JOPT  input j 

READC a, 106)  JOPT 

106  F0RI.1AT(I2) 
IF(J0PT.EQ.2)G0T0  1 

This  variable  denotes  the  option  to  be  used  in  entering 
vehicle  performance  data.   If  JOPT  equals  1,    the  maximum 
sustained  speed  is  input.   If  JOPT  equals  2,  the  maximum 
sustained  horsepower  and  machinery  characteristics  are  input. 

Vehicle  performance  input  for  JOPT  equal  lx 

READ(«,107JVSUS, VEND, RGEND.PCE, PCM 
10?  F0RMAT(2F5.0,F10.0,2F5.0) 
GOTO  9 
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The  maximum  sustained  speed  is  input  along  with  the 

endurance  speed  and  endurance  range.   The  propulsive 

coefficients  at  the  endurance  speed  and  at  maximum  speed 

must  also  be  given. 

Vehicle  performance  input  for  JOPT  equal  Zt 

1  KEAD(a,10b)HP3S5StDHMtXLM,WTG^5FCMHP.SFCHHPtVENDtRGEND, 
PAL0AD,PCE,PCM,CSTG2,WTL0 

100  F0KMAT(7F5.0,F10.0,5F5.0) 
HPMS=HPAK>*1000, 

CSTG2=CSTG2*1000. 
GR2CST=CSTG2 

This  option  is  used  primarily  when  a  specific  machinery- 
plant  is  desired.   Besides  the  maximum  sustained  horsepower, 
the  length  and  minimum  depth  of  the  engine  room  must  be 
given.   The  total  group  2  weight  and  specific  fuel  con- 
sumptions at  maximum  and  half  power  are  required.   Tne 
endurance  speed  and  range  must  be  given  as  before.   Pro- 
pulsion auxiliaries  electrical  load  and  propulsive  coef- 
ficients at  full  and  endurance  speed  are  required  as  are 
the  group  2  material  cost  and  lub  oil  weight. 

To  fit  into  the  format  used,  the  horsepower  and 
cost  values  must  be  divided  by  100U.  The  cost  of  group 
2  must  be  stored  for  later  use  as  were  the  costs  of  group 
M-  and  group  7. 

Options  input: 

9  RE AD ( « , 1 09  JMTYPE , FR3C , LEN , C P , C X , GMMUL , DBMAR , EELCF 
109  F0RMAT(l5,6F5.0,FlO.Oj 

Machinery  type  must  always  be  specified  even  if  the 

option  of  inputmg  machinery  characteristics  was  used. 

If  machinery  information  was  given,  any  number  except 
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zero  can  be  used  for  KTYPE.   A  value  or  zero  will  cause 

the  program  to  assume  the  next  input  value  is  INDEX. 
This  can  be  used  to  return  and  stop  the  program  or  to 
input  a  new  set  of  data  for  another  ship.   If  the  machinery 
particulars  have  not  been  specified,  MTYPE  equals  zero 
has  the  same  effect.   However,  MTYPE  equals  1  will  cause 
the  program  to  assume  z.   diesel  plant  is  to  be  installed. 
A  value  of  MTYPE  equal  to  2  will  result  in  calculations 
lor  a  CODOG  plant.   These  are  the  only  two  machinery 
types  for  which  estimating  relationships  have  been  included 
in  the  program,   if  MTYPE  is  any  other  number,  the  results 
are  unpredictable. 

A  free  surface  correction  can  be  specified.   This 
value  is  interpreted  as  a  rise  in  the  full  load  KG  of  the 
ship  in  an  amount  equal  to  FRSC  in  feet.   A  negative  value 
can  also  be  used.   In  this  case  FRSC  can  be  interpreted 
as  a  lowering  of  the  average  KG  value  built  into  the 
program. 

The  length  between  perpendiculars  is  given  next.   The 
values  for  Cp  and  Cx  have  the  standard  naval  architectural 
definitions.     GMMUL  if  the  fraction  GM/Beam,  A  value 
of  0.1  is  commonly  used.   The  next  variable,  DBMAK,  is 
a  design  and  builders  margin.   This  is  given  as  a  fraction, 
i.e.,  a  10  percent  margin  is  input  as  0.1.  DELCF  is  the 
ACr  value  used  in  the  calculation  of  horsepower.   A  value 
of  0.0004  is  common. 

IF(MTYPE.SQ.0J  GOTO  2 
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One  or  more  option  input  cards  may  follow  the  first. 

The  program  will  assume  that  each  one  is  a  new  input 
case  with  all  other  data  the  same.   To  input  an  entirely- 
new  set  of  data,  an  option  input  card  with  a  value  of 
MTYPE  equal  to  zero  should  be  included.   No  other  values 
need  be  given  on  the  card.   This  will  cause  the  next 
input  card  tc  be  interpreted  as  a  new  value  for  INDEX. 

DtfKAR=l.+DBMAR 

CSTG2=GR2CST 

CSTG4=G4CST 

CSTG7=G7CST 

DBKAR  is  changed  to  a  multiplier  for  weight  and  the 

input  costs  are  reset  to  their  initial  values. 

Initial  KG  estimate: 

KGTRY=LEN/( . 04b>*LEN+4. U25) 
KGTRY=KGTRY-2. 

The  initial  KG  value  is  estimated  based  on  the  values 
for  current  Coast  Guard  vessels.   These  are  shown  in 
Figure  4.   Reducing  the  initial  estimate  by  2  feet  was 
found  to  be  advantageous.   It  helps  to  insure  that  the 
first  iteration  will  not  exceed  the  KG  constraints  of 
the  program. 

Initial  displacement  estimate: 

DPTRY=(.006562-.U01ia«*LEN/100. )*LEN**3/35. 
IF(INDEX.EQ.2)DPTRY=.0023*LEN**3/35. 

Two  initial  values  of  displacement  are  available. 

For  an  INDEX  value  other  than  2  the  displacement  is  taken 

as  that  of  conventional  Coast  Guard  vessels.   This 

displacement  curve  is  shown  on  Figure  4.   However,  if  the 
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vessel  has  a  high  speed-length  ratio,  greater  than  1.3, 

it  will  not  be  possible  to  find  the  horsepower  in  some 
cases  if  the  conventional  ship  curve  is  used.   This  is 
due  to  the  limitations  on  the  Cr  arrays.   By  specifying 
an  INDEX  equal  to  2  the  program  assumes  a  fine  form  at 
least  for  the  first  iteration. 

Cubic  number  estimate: 

CN=DPTRY*lbO 

An  initial  cubic  number  is  assumed  for  use  in  estmating 
the  electrical  load  on  the  first  iteration.   Values  for 
current  Coast  Guard  vessels  range  from  140  to  180  times 
the  displacement. 

CALL  XECUTE(HPMS) 

Control  is  transferred  to  subroutine  XECUTE  to 

calculate  the  vessels  characteristics. 

GOTO  y 
END 

The  GOTO  statement  causes  the  program  to  loop  back 

and  try  to  read  another  options  input  data  card.   The 

END  statement  stops  compilation  of  the  routine. 


3. 2. ^     List   of   Outputs    of  MAIN 

ACCG(20) LC(EE) 

ACWT(20) LC(EE) 

AMCCG(20) LC(EE) 

AM0WT(20) LC(EE) 

AREADK LC(BB) 

AREAHL LC(3B) 

BLOAD LC(DD) 

CACFUL LC(HH) 

CARG0C(20) ...LC(EE) 

CARG0W(20) LC(EE) 

CN LC(BB) 

CP LC(GG) 

CR1(1152) LC(TT) 

CR2(972) LC(TT) 

CR3(2l6) LC(TT) 

CSTG2?* LC(tfM) 

CSTG4 . ..LC(MM) 

CSTG7 LC(MM) 

CX LC(GG) 

DBMAR LC  (WW) 

DELCF LC(CC) 

DHA(20) LC(KK) 

DHAR(20) LC(KK) 

DHM?T LC(JJ) 

DOLHR LC(NN) 

DPTRY LC(GG) 


ifO 
Subroutine 

ELD(20) LC(KK 

ELLD(20) LC(KK 

ELOAD LC(DD 

ENDDAY LC(AA 

FRSC LC(GG 

G1IND LC(NN 

G2IND LC(NN 

G3IND LC(NN 

G5IND LC(NN 

G6IND ..LC(NN 

GMMUL LC(GG 

GR4CG(20) LC(EE 

GR^CST(20) LC(KK 

GR^WT(20) LC(EE 

GR?CG(20) LC(EE 

GR7CST(  20 ) LC  (KK 

GR7WT(20) LC(£E 

HEAD(20) LC(FF 

HLA(20) LC(KK 

HLAR(20) LC(KK 

HPMS SCDA 

JOPT LC(CC 

KGTRY LC(GG 

LEN LC(BB 

MTYPE LC  (CC 

NCPO LC(AA 


List  of   Outputs   of  MAIN  Subroutine    (cont) 


kl 


NENL • LC(AA) 

NOARM LC(FF) 

NQELT LC(FF) 

NOFF LC  ( AA) 

PALOADt? .LC(DD) 

PCE LC(CC) 

PCM LC(CC) 

RGEND LC(CC) 

SFCHHPTT LC(CC  ) 

SFCMHPT? LC(CC) 

TITLE(5.^0) LC(FF) 


VEND LC(CC) 

VSUST LC(CC) 

WACFUL LC(LL) 

WCARGO LC(LL) 

WTAC LC(LL) 

WTAMMO LC(LL) 

WTG2?T LC(WW) 

WTG4 LC(WW) 

WTG7 LC(WW) 

WTLO LC(LL) 

XLM LC(JJ) 


*  Output  only  if  JOPT  «  1 
**  Output  only  if  JOPT  «  2 
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3.2.5  Nomenclature  List 

ACCG(20)       aircraft  vcg/DIO,   D10  »  depth  amidships 

aircraft  weight,  tons 

ammunition  vcg/DIO 

ammunition  weight,  tons 

total  input  deckhouse  arrangements  area,  sq  ft 

total  input  hull  arrangements  area,  sq  ft 

average  electrical  load,  KW 

maximum  beam,  feet 

total  armament  electrical  load,  KW 

aircraft  fuel  vcg/DIO 

cargo  vcg/DIO 

cargo  weight,  tons 

total  vcg  of  cargo  inputs,  feet 

construction  services  cost,  dollars 

total  vcg  of  aircraft,  feet 

total  vcg  of  ammunition,  feet 

total  vcg  of  crew,  feet 

total  vcg  of  ships  fuel,  feet 

total  vcg  of  lub  oil,  feet 

total  vcg  of  personal  effects,  feet 

total  vcg  of  personnel  stores,  feet 

cubic  number,  cubic  ft 

prismatic  coefficient 

Cr  arrays  for  Taylor's  Standard  Series 

resistance  estimate 


ACWT(20) 

AMOCG(20) 

AM0WT(20) 

AREADH 

AREAHL 

AVGKW 

B 

BLOAD 

CACFUL 

CARG0C(20) 

CARG0W(20) 

CCARGO 

COST 

CGAC 

CGAMMO 

CGCREW 

CGFUEL 

CGLO 

CGFE 

CGPS 

CN 

CP 

CR1(1152) 

CR2(972) 

CR3(2l6) 


Nomenclature 

CSTG1 

CSTG2 

CSTG3 
CSTG4 

CSTG5 
CS  TG6 
CSTG? 
CWP 

cx 

DAVG 

DBMAR 

DCST 

DELCF 

DHA(20) 

DHAR(20) 

DHM 

DHV 

DOLHR 

DPTRY 

ELD(20) 

ELKW 

ELLD(20) 

ELOAD 

ENCVOL 

ENDDAY 

FRSC 
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List  (cont) 

total  cost  of  weight  group  1,  dollars 

total  cost  of  weight  group  2,  dollars 

total  cost  of  weight  group  3»  dollars 

total  cost  of  weight  group  k,   dollars 

total  co3t  of  weight  group  5»  dollars 

total  ccjt  of  weight  group  6,  dollars 

total  cost  of  weight  group  7,  dollars 

waterplane  coefficient 

midship  section  coefficient 

average  hull  depth,  feet 

design  and  builders  margin 

design  cost,  dollars 

ACf   used  in  resistance  calculation 

electronics  deckhouse  area,  sq  ft 

armament  deckhouse  area,  sq  ft 

minimum  machinery  box  depth,  feet 

deckhouse  volume ,  cu  ft 

labor  rate ,  dollars/hour 

displacement  estimate,  tons 

electronics  electrical  load,  KW 

generator  rated  load,  KW 

armament  electrical  load,  KW 

total  electronics  electrical  load,  KW 

total  volume  of  hull  and  deckhouse,  cu  ft 

number  of  endurance  days  for  dry  provisions 

free  surface  correction,  feet 


Nomenclature 

G1IND 

G2IND 

G3IND 

G4CST 

G5IND 

G6IND 

G7CST 

GMMUL 

GR2CST 

GR4CG(20) 

GR4CST(20) 

GRWT(20) 

GR7CG(20) 

GR7CST(20) 

GR7WT(20) 

H 

HEAD(20) 

HLA(20) 

HLAR(20) 

HPMS 

INDEX 

JOPT 

KGTRY 

LEN 

LRD 

MCST 


List  (cont) 

cost  index  for  weight  group  1 

cost  index  for  weight  group  2 

cost  index  for  weight  group  3 

total  material  cost  weight  group  4,  dollars 

cost  index  for  weight  group  5 

cost  index  for  weight  group  6 

total  material  cost  weight  group  7»  dollars 

GM/Beam  ratio 

total  material  cost  weight  group  2,  dollars 

group  k   vcg/DIO 

material  cost  weight  group  4,  dollars 

group  4  weight,  tons 

group  7  vcg/DIO 

material  cost  weight  group  7#  dollars 

group  7  weight,  tons 

draft,  feet 

descriptive  heading  (alphanumeric) 

electronics  hull  area,  sq  ft 

armament  hull  area,  sq  ft 

maximum  sustained  shaft  horsepower 

control  variable 

input  option  variable  for  machinery 

initial  KG  estimate,  feet 

length  between  perpendiculars,  feet 

raised  deck  length,  feet 

miscellaneous  costs,  dollars 


Nomenclature 

MTYPE 

NCPO 

NENL 

NOARM 

NQELT 

NOFF 

NOCPO 

NOCREW 

NOOFF 

PALOAD 

PCE 

PCM 

RGEND 

SFCHHP 

SFCMHP 

SHFE 

SHPM 

TITLE(5.^0) 

TOTCST 

VEND 

VSUS 

WACFUL 

WCARGO 

WFULLD 


List  (cont) 

machinery  type  option 

number  of  chief  petty  officers  (real) 

number  of  enlisted  men  (real) 

number  of  armament  input  cards 

number  of  electronics  input  cards 

number  o£   officers  (real) 

number  of  chief  petty  officers  (integer) 

number  of  enlisted  men  (integer) 

number  of  officers  (integer) 

propulsion  auxiliary  electrical  load,  KW 

propulsive  coefficient  at  endurance  speed 

propulsive  coefficient  at  maximum  speed 

endurance  range ,  nautical  miles 

specific  fuel  consumption  at  half  horsepower, 

lbs/SHP-hr 

specific  fuel  consumption  at  maximum 

horsepower,   lbs/SHP-hr 

shaft  horsepower  at  endurance  speed 

shaft  horsepower  at  maximum  speed 

line   titles . (alphanumeric) 

total  lead  ship  cost,   dollars 

endurance  speed,  knots 

maximum  sustained  speed,  knots 

weight  of  aircraft  fuel,  tons 

total  cargo  weight,  tons 

full  load  weight,  tons 


Nomenclature 

WLSHIP 

WTAC 

WTAMMO 

WTCRSW 

WTFUEL 

WTG1 

WTG2 

WTG3 

WTG4 

WTG5 

WTG6 

VTG7 
WTLO 
WTPE 
WTPS 
TOM 


List  (cont) 

light  ship  weight,  tons 
total  aircraft  weight,  tons 
total  ammunition  weight,  tons 
total  crew  weight,  tons 
total  ships  fuel  weight,  tons 
total  group  1  weight,  tons 
total  group  2  weight,  tons 
total  group  3  weight,  tons 
total  group  4  weight,  tons 
total  group  5  weight,  tons 
total  group  6  weight,  tons 
total  group  7  weight,  tons 
weight  of  lub  oil,  tons 
weight  of  personal  effects,  tons 
weight  of  personnel  stores,  tons 
machinery  box  length,  feet 
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3.3  Subroutine  XECUTE 

3. 3.1  Introduction 

Subroutine  XECUTE  controls  the  execution  of  each 
trial  case.  All  of  the  remaining  subroutines  except 
two  are  called  by  this  routine.  Error  messages  are  also 
printed  by  this  routine  whenever  a  constraint  of  the 
program  is  violated. 

In  three  cases,  iteration  control  statements  have 
been  included  in  XECUTE.  The  first  set  of  these  statements 
controls  the  convergence  of  the  speed-power  estimate 

when  maximum  sustained  horsepower  is  specified.  A  modified 

(2) 

Newton-Raphson  technique   '  is  used  with  the  slope  ap- 
proximated by  a  secant  to  the  speed-power  curve. 

The  second  set  of  iteration  control  statements 
controls  the  convergence  between  estimated  displacement 
and  full  load  weight.  A  similar  modified  Newton-Raphson 
technique  is  used. 

The  last  set  deals  with  the  convergence  of  the 
vertical  center  of  gravity  estimate.  For  this  set  there 
is  a  tendency  for  the  values  to  oscillate  and  a  Newton- 
Raphson  technique  could  not  be  used.  The  average  of  the 
estimated  and  calculated  values  of  KG  is  used  instead. 
While  this  eliminates  the  oscillations,  convergence  may 
be  slow  in  some  cases. 

A  flow  chart  for  this  subroutine  is  shown  in  Figures 
5a,  5b,  5c,  and  5<i« 


XECUTE   SUBROUTINE 
FLOW   CHART 
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o  p  i  s  dptr y 

DP7T?y=  \A/FULLJ> 


Figure      5c 


51 


c 


f-(     CALL     CoSUueN) 


CALL  S£ASPS>(AVS,F} 


(CALL  OOTPOT 


c 


&E.TUZ.KJ 


3 


F  igu  re     5d 


3,3*2  Input  List  for  Subroutine  XECUTE 

DPTRY LC(GG) 

HPMS SCDA      Only  if  JOPT  *  2 

JOPT LC(CC) 

KGTRY LC(GG) 

LEN LC(BB) 

MTYFE LC(CC) 

PCE LC(CC) 

PCM LC(CC) 

YE LC(CC) 

VSUS LC(CC)      Only  if  JOPT  =  1 

WFULLD LC(WW) 

3*3*3  Statement  Descriptions 

REAL  and  COMMON  statements  will  not  be  listed  with 

individual  subroutines  but  may  be  found  in  Appendix  C. 

EXCKGaO. 

JNDEX=»1 

LNDSX-1 

These  statements  initialize  values  for  the  excess 

KG  variable,  the  VCG  loop  counter  and  the  weight  loop 

counter  respectively, 

7  CALL  UWDIM(R) 

This  subroutine  is  called  to  calculate  the  beam 

and  draft  of  the  ship.   The  variable  R  transfers  the 

value  of  any  excess  KG, 

IF(LEN.LT.O.)  GOTO  999 

If  a  beam  to  draft  ratio  of  greater  than  4,0  is 
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required  to  satisfy  the  stability  requirements,  LEN  is 

set  equal  to  -1000.  in  UWDIM.   This  signals  XECUTE  to 
print  an  error  message  and  return  control  to  MAIN, 

IF(J0PT.EQ.2)  GOTO  1 

The  next  group  of  statements  control  the  speed- 
power  estimate.      Two  procedures  are  used  depending  on 
which  performance   input  option  is  specified. 

CALL  HPCALC(VSUS,SHP,PCM) 
IF(LEN.LT.O.)   GOTO  998 
GOTO  2 

These  statements  are  used  if  the  sustained  speed 

is  specified.     Subroutine   HPCALC  is   designed  to  return 

a  value   of  shaft  horsepower,  SHP,   if  speed  and  propulsive 

coefficient  are   specified.      In  this  case  VSUS  is  the 

speed  and  PCM  the   propulsive  coefficient  used. 

LEN  is  used  here  also  as  a  signal  that  an  error 

message   is  required.     A  negative  length,  LEN,   will  occur 

if  the  limits   of  the  Cr  array  are  exceeded.      The  ranges 

of  the  Cr  array  are i 

.48  <  Cp  <   .70 

2.   <  B/T  <  4. 

for  .001<  Cv  <.003 

•  5  <  v/7l  <  2. 
for  .003<  Cy  <.006 

•  5  <  v/Vl  <1.3 

The  GOTO  2  statement  causes  the  routine  to  skip 
the  speed  estimate  used  when  JOPT  «  2. 
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1  VMAX  «  (HPMS*200./LEN)**. 33333 

If  the  maximum  sustained  horsepower,  HPMS,  was 

specified,  the  following  statements  are  used  for  estimating 

the  maximum  sustained  speed.   The  statement  given  above 

provides  a  rough  first  guess  at  the  speed, 

k   CALL  HPCALC( VMAX, SHP, PCM) 
IPCLEN.LT.O.)  GOTO  998 
SHP1=SHP 

VMAX1».95*VMAX 

CALL  HPCALC(VMAX1,SHP,PCM) 

IF(LEN.LT.O)   GOTO   998 

SHP2=SHP 

Two  points  on  the   speed-power  curve  are  calculated 
in  the  above  statements,   one  point  at  VMAX  and  one  point 
at  ,95  x  VMAX.      The  constant,    ,95#  was  chosen  arbitrarily 
to  give  a  secant  to  the   curve  which  has  a  slope  close   to 
the  actual  slope  at  VMAX.      The  same  limits  of  the  Cr  array 
apply  here. 

DELHP-  ( SHP1 -SHP2 )/( VMAX-VMAX1 ) 

This  is   the  approximate  value   of 

*lgHP>        at  VMAX 

dv 

VMAX1=VMAX- ( SHP1  - .  8*HPMS  )/DELHP 

By  definition  the  maximum  sustained  speed  is   the 
speed  that  can  be  made  when  the  shaft  horsepower  is   80 
percent  of  the  maximum  sustained  horsepower.     This  assumes 
the  ship  has  a  clean  bottom  and  is  operating  in  calm 
water.     Therefore,   the  function  to  be  solved  isi 

f(V)   -  SHP(V)   -    .8  x  HPMS  «   0 

The  Newton-Raphson  method  is  an  iterative  procedure 


vk 
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in  which » 

f(Vk) 
vk+l  "  Vk  "  f(Vk) 

where  f(Vk)  -  V^ 

In  the  present  case : 
Vk  *  YMAX 

vk+l  s  VMAX1 

f(Yk)  »  SHP1  -  .8  x  HPMS 
f*(Vk)  =  DSLHP 

DELHP  is  only  an  approximation  to  f*(Vk)i  however, 
it  produces  very  good  convergence • 

IF(ABS(VMAX1-VMAX).LE.VMAX1».005)  GOTO  3 

The  iteration  stops  when  the  difference  between 

any  two  successive  speed  estimates  is  less  than  half  a 

percent  of  the  best  estimate.   The  above  iteration  technique 

is  very  fast  and  convergence  to  these  limits  usually 

occurs  in  3  to  4  iterations. 

YMAX  »  VMAX1 
GOTO  k 

These  statements  are  used  if  the  two  speed  estimates 

are  not  within  the  above  limit.   The  new  velocity  estimate 

is  used  and  another  iteration  is  run. 

3  VSUS  -  VMAX1 

SHPM  -  HPMS 
GOTO  5 

The  above  pair  of  values  are  saved  for  output  and  for 

use  in  other  subroutines. 

2  SHPM  =  SHP/.8 

In  the  case  where  VSUS  is  specified,  the  value  SHP 
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is  equal  to  SHPM  tines  .8  as  before. 

IF(MTYFE. EQ. LAND. SHPM. GT. 14000.)  GOTO  997 
If  a  diesel  plant  is  specified,  horsepowers  of 
greater  than  14000  are  outside  the  range  of  validity 
of  the  model.  An  error  message  is  printed  if  this  occurs. 

5  CALL  HPCALC(VE,SHP,PCE) 
IF(LEN.LT.O.)   GOTO  998 
SHFE-SHP 

The  endurance  horsepower  is  calculated  without  regard 
to  fouling  or  sea  conditions.   This  statement  is  executed 
for  both  performance  options. 

CALL  EPLANT 

This  subroutine  is  called  to  calculate  the  average 
electrical  load  at  sea  and  to  size  the  ships  generators* 

CALL  MACHLQ 

This  subroutine  is  called  to  calculate  the  weight 

of  fuel,  lub  oil  and  potable  water  required. 

IF(J0PT.EQ.2)  GOTO  6 
CALL  MBSIZE(B) 

Subroutine  MBSIZE  calculates  the  length  and  minimum 

depth  of  the  machinery  box.   If  JOPT  »  2  these  values 

are  given  as  input. 

6  CALL  VOLUME 
IF(LEN.LT.O)  GOTO  996 

Subroutine  VOLUME  sets  the  sheer  line  and  balances 
the  required  area  and  available  area.  If  insufficient 
area  is  available,  the  value  of  LBN  will  be  negative. 

CALL  WEIGHT 

All  the  weights  not  previously  estimated  or"  inputed 
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are  calculated  in  this  routine.   The  primary  value  returned 

is  WFULLD,  the  full  load  weight. 

IF(ABS(WFULLD-DPTRY).LT.2.)  GOTO  8 

The  calculated  full  load  weight  is  compared  to  the 

displacement  estimate  made  earlier.     If  the   two  do  not 

agree  within  2  tons,   a  new  estimate   of  displacement  is 

made  and  the  program  re  jurns  to  call  UWDIM.     If  there   is 

agreement,   the   program  proceeds  to  calculate  vertical 

center  of  gravity. 

IF(LND2X.GE.2)    GOTO  9 
WT1« WFULLD 
DPlcDPTRY 
DPTRY-WPULLD 
GOTO  10 

Two  values   of  weight  and  two  values  of  displacement 

are  needed  for  the  Newton-Raphson  procedure  explained 

below  to  work.      On  the  first  iteration,   the  above  simplified 

procedure  is  used. 

9  FACTOR" (DPTRY-WFULLD-DP1+WT1 )/(DPTRY-DPl) 
DPlrsDPTRY 
WT1=WFULLD 
DPTRY*DPXR  Y-  ( DPTRY-VfFULLD  )/FAC  TOR 

In  this  case i 

f(DPTRY)  »  DPTRY  -  WFULLD(DPTRY)  «  0 
Again  a  secant  approximation  is  made  to  f* (DPTRY).   This 
is  the  variable  FACTOR. 

DPTRYk+1  -  DPTRYk  -  fflggfo 

10  IF(LNDEX.GT.20)  GOTO  995 

In  case  the  iteration  does  not  converge,  this  statement 
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provides  a  safety  valve   to  stop  the  loop.     An  error 

message   is  printed. 

LNDEX=LNDEX+1 
GOTO  7 

A  new  iteration  is  run  beginning  with  a  call  to  UWDIM. 

8  CONTINUE 
LNDEX-1 
CALL  VCG(CFULLD) 

Once   the  full  load  weight  is  sufficiently  close  to 
displacement,   the  program  proceeds  to  calculate  KG.     The 
weight  counter  is  reset.     The  value  CFULLD  is  the  calculated 
KG.     This  value  is  compared  with  the  previous  estimate, 
KGTRY. 

IF(ABS(CFULLD-KGTRY).LT..l)   GOTO  11 

This  iteration  stops  when  the  two  KG  values  are  less 
than  0.1  feet  apart. 

IF(JNDEX.GT.2l)  GOTO  99^ 

Again  a  safety  valve   is  provided  in  case   the  iteration 
does  not  converge. 

IF (CFULLD.LT. KGTRY)   GOTO  12 
13  KGTRY=(CFULLD+KGTRY)/2. 

JNDEX-JND2X+1 

GOTO   ? 
12  IF(R.LT..001)   GOTO  13 

EXCKGaKGTRY-CFULLD 

In  some  cases   the  minimum  beam  to  draft  ratio  may 
provide  more  stability  than  is  required  to  satisfy  the 
GM  criterion.     The   ship  is  still  feasible  if  this  occurs, 
but  may  be   over  designed.     When  this  occurs  R  will  have 
a  value  greater  than  0.     In  fact,   it  will  be  equal  to  the 
excess  KG  when  KGTRY  is  used  in  UWDIM.      The  input  value 
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of  KGTRY  is  modified  in  UWDIM  by  adding  R.   Therefore, 

the  value  of  KGTRY  at  this  point  in  the  program  is  the 
maximum  value  of  vertical  center  of  gravity  for  which 
a  minimum  beam  to  draft  ratio  will  still  be  valid.   If 
CFULLD  is  less  than  this  value  there  is  no  point  in  doing 
another  iteration.   Instead,  the  program  continues  and 
eventually  prints  the  excess  KG  available  as  an  output. 
In  all  other  cases  a  new  estimate  of  KGTRY  must  be 
made.  As  mentioned  earlier,  this  estimate  is  taken  as 
the  average  of  the  initial  estimate  and  final  calculated 
value.   This  is  done  to  prevent  the  two  values  from 
oscillating.  As  an  example,  in  the  case  of  the  378* 
WHEC,  a  KGTRY  of  1^  feet  returned  a  CFULLD  of  16  feet. 
A  KGTRY  of  16  feet  returned  a  CFULLD  of  14  feet.   This 
oscillation  may  not  occur  in  othor  cases,  but  the  averaging 
estimate  used  will  converge  in  every  case,  although  pos- 
sibly at  a  slow  rate. 

11  CONTINUE 

CALL  COST(LEN) 

Subroutine  COST  calculates  any  remaining  cost  items 
not  given  as  input. 

CALL  SEASPD(AVSP) 

This  subroutine  estimates  an  average  speed,  AVSP, 
in  the  North  Atlantic  Ocean  which  includes  the  effects 
of  sea  state. 

CALL  0UTPUT(AVSP,EXCKG,2) 
Subroutine  OUTPUT  contains  the  write  statements  for  printing 
all  the  output  and  input  data.   The  value  2  causes  OUTPUT  to 
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print  both  the  input  and  output  data* 

RETURN 

The  present  input  case  has  been  completed  and  control 
is  returned  to  the  MAIN  program  for  more  data. 

The  error  message  statements  are  all  similar  in  format 
and  generally  self  explanatory.   Only  one  typical  one 
will  be  discussed. 
99^  CALL  0UTPU?(0.,0.,1) 

Subroutine  OUTPUT  is  called  with  a  value  of  1.  This 

directs  OUTPUT  to  print  only  the  input  data. 

WRITE(5,106) 
106  FORMAT (///'  NO  BALANCE  BETWEEN  ASSUMED  AND  CALC  KG 

WITHIN  0.1  FEET  COULD  BE  MADE  IN  20  ITERATIONS*///) 
GOTO  1000 
1000  WRITE(5,102) 
102  F0RMAT(' PROGRAM  PROCEEDING  TO  NEXT  INPUT  CASE*///) 
RETURN 
END 

The  program  prints  the  appropriate  error  message 

and  returns  for  more  data.   This  concludes  subroutine 

XECUTE. 

3.3.4  Output  List  for  Subroutine  XECUTE 

DPTRY LC(GG) 

EXCKG SCDA 

KGTRY .LC(GG) 

SHPE LC(CC) 

SKPM LC(CC) 
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3.3*5  Nomenclature  List 

All  variables  have  the  same  definition  as  given  in 

the  MAIN  program  nomenclature  list  except  for  the  following* 


AVSP 

CFULLD 

DELHP 

DPI 

EXCKG 

FACTOR 

JNDSX 

LNDEX 

R 

RE 

SFCH 

SFCM 

SHP 

SHP1 

SHP2 

VE 

VMAX 

YMAX1 

WT1 


average  speed  in  N.A.  Ocean,  knots 

vertical  center  of  gravity,  KG,  feet 

derivative  of  horsepower  function 

used  to  store  DPTRY  value,  tons 

excess  KG,  feet 

derivative  of  weight  function 

vcg  counter 

weight  counter 

excess  KG,  feet 

same  as  RGEND  in  MAIN 

same  as  SFCHHP  in  MAIN 

same  as  SFCMHP  in  MAIN 

dummy  used  for  shaft  horsepower 

dummy  used  for  shaft  horsepower 

dummy  used  for  shaft  horsepower 

same  as  VEND  in  MAIN 

dummy  used  for  speed,  knots 

dummy  used  for  speed,  knots 

used  to  store  WFULLD,  tons 
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3.*f  Subroutine  UWDIM 

3,4. 1  Introduction 

Three  of  the  important  variables  determining  the 
underwater  hull  dimensions  are  given  as  input.  These  are 
length,  Cp  and  Cx.  Two  others,  the  beam  and  draft,  are 
calculated  in  UWDIM.   The  waterplane  coefficient,  G^, 
is  also  calculated,  but  only  as  a  linear  function  of 

cp. 

Displacement  is  an  input  to  this  subroutine  and 
therefore  the  value  of  beam  times  draft  is  fixed  since t 

B  x  T  »  Displacement  x  35«/(Lbp  x  cp  x  cx)    Ref.(3) 

Reference  (3)  may  be  referred  to  for  a  detailed 
definition  of  the  naval  architectural  variables  used 
in  this  routine. 

The  beam  to  draft  ratio,  B/T,  is  determined  by  the 
GM/B  criterion.   The  GM  available  is  given  by  the  formula i 

GM  a  KB  +  BH  -  KG  -  Free  surface  correction   Ref.(3) 
Values  for  KG,  free  surface  correction  and  GM/B  are 

inputs  to  this  routine.   A  value  for  KB  is  estimated 

(3) 

using  Morrish's  approximate  formula* 

KB  =  T  -  1/3  x   (T/2  +  V/A) 
where      T  is   the   draft  in  feet 

V  is   the  volume   of  displacement  ■ 

B  x  T  x  Lpp  x  C     x  Cx   ,   cu  ft 
A  is   the  waterplane   area  a  B  x  LDp  x  C^,   sq  ft 
Substituting  for  V  and  A  gives i 

KB  »  T  x   (5/6  -    (Cp  x  Cx)/(3  x  C^)) 
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A  value  for  BM  is  estimated  using* 

BM  =  LfrP  x  B^  x  Ca 

where  Ca  is  a  transverse  moment  of  inertia  coefficient 
which  is  assumed  to  be  a  linear  function  of  Cwp» 

A  number  of  coefficients  are  used  in  this  routine 
which  are  defined  by  the  above  equations.   These  are 
listed  below, 

CI  =  DPTRY*35./(LEN*CP*CX)   «  B  x  T 

C2  =  .833  -  CP*CX/(3*CWP)   =  KB/T 

C3  «  LEN*CALPH/(DPTRY*35.)   =  Bty/B3 

Ck   ■  KG  +  FRSC 

C5  =  GM/B  =  GMMUL 

Figure  7  gives  a  typical  plot  of  the  equation* 

KB  +  BM  -  GM  =*  C2  x  T  4-  C3  xCl3  -  CI  x  C5 

^  T 

The  curve  shown  was  calculated  using  data  for  a 
3?8f  WHEC. 

The  value  of  KG  is  also  plotted  in  Figure  7.   In 
the  discussion  below,  KG  refers  to  the  value  of  KG  +  FRSC. 
A  variable  R  is  now  introduced  with  a  value  of i 

R  ■  KB  +  BM  -  GM  -  KG^y 
When  R  ■  0,  the  stability  criterion  that  GM  »  C5  x  B  is 
just  satisfied.   This  is  the  desired  solution.   This 
point  is  shown  as  point  3  in  Figure  7. 

However,  the  solution  is  only  valid  provided  point  3 
lies  within  satisfactory  limits  of  beam  to  draft  ratio. 
These  limits  are  specified  by  the  speed-power  estimation 
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kv 
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to  be  i 

2.  <  B/T  <  ^. 
or  since  B  «  Cl/T 

Vc!7S  <  T  <  V^l72 

It  can  be  seen  from  Figure  7  that  if  the  value  of 
KG  is  lower  than  the  value  of  KB  +  BM  -  GM  at  point  1 
there  would  be  no  draft  for  which  R  ■  0,   The  minimum 
GM  that  could  be  provided  would  still  be  greater  than 
C5  x  B. 

The  solution  method  used  in  UWDIM  proceeds  as  follows i 

(1)  The  value  of  B/T  at  point  1  is  calculated. 
Then  a  check  is  made  to  see  if  point  1  lies  between  the 
acceptable  limits  of  B/T, 

(2)  If  point  1  lies  within  the  limits  of  B/T,  the 
value  of  R  is  calculated.   If  R  is  positive,  no  solution 
exists  to  Re  0  but  the  ship  may  still  be  feasible.  In 
this  case  the  minimum  GM  is  greater  than  C5  x  B 
this  occurs  the  KG  estimate  is  raised  until  R  ■  0  at 
point  1  and  the  program  returns  to  subroutine  XECUTE. 

If  it  is  found  later  that  the  actual  KG  of  the  ship  is  ^.t 
lower  than  the  new  KG  estimate,  the  difference  is  printed 
as  an  excess  KG, 

If  the  value  of  R  at  point  lis  negative ,  the  program 
executes  the  procedure  described  in  step  4, 

(3)  If  point  1  does  not  lie  within  acceptable 
limits  of  B/T,  the  value  of  R  at  point  2  is  calculated. 
If  R  is  positive  at  this  point  the  procedure  is  the  same 
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as  that  outlined  in  step  2.   If  R  is  negative,  execution 

proceeds  to  step  4. 

(*}•)  If  the  value  of  KG  is  greater  than  the  value 
of  KB  +  BM  -  GM  at  point  1  or  point  2,  the  value  of  R 
at  point  k   is  calculated*   If  the  value  of  R  is  negative 
at  this  point,  no  draft  can  be  specified  within  the 
limits  of  B/T  that  will  satisfy  the  stability  criterion 
and  the  program  will  have  to  go  on  to  another  input  case. 

If  R  is  positive  at  point  4,  the  solution  must  lie 
between  point  4  and  point  1  or  point  2,  whichever  is 
lower.   It  is  assumed  that  if  point  1  is  within  the 
limits  on  B/T  that  the  solution  does  not  lie  between 
point  1  and  point  2. 

(5)  To  solve  for  the  draft  at  which  R  ■  0,  a  synthetic 
division  technique  is  used  together  with  the  Newton- 

Raphson  iteration  technique.   This  is  sometimes  known  as 

(2) 

the  Birge-Vieta  method.  The  equation  for  R  must  be 

modified  slightly  to  use   this   technique. 

f (T)   »  T%/C2  =  T^   -  C4T3/C2  -CjC^T2/^  +  c3cl"Vc2  *  ° 
In  the  Birge-Yieta  method  three  columns  are  used  as 
shown  below i 

111 

ai      h      ci 

a2       J>2        c2 
ll      \l      CJ 

The  values  in  the  first  column  are   the   coefficients 
in  the   equation  for  f(T),  namelyi 

al  ■  -C^/^,   a2  «  -0^5/02 .  a^  «  0,  a^  *  O^C^/0Z 
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The  values  in  columns  two  and  three  are  calculated 

using  the  equations i 

b^  a  a^  +  z  x  b,  .. 

Ci  ■  b^  +  z  x  ci-l 

where  z  is  the  current  estimate  of  the  draft $   T. 

The  new  estimate  of  draft  is  them 

znew  ■  z  '   Vc3 
New  values  are  computed  for  columns  two  and  three  using 

znew.   The  procedure  starts  with  the  draft  at  point  b 

and  continues  until  the  difference  between  zn0w  and  z 

is  sufficiently  small.  This  iteration  procedure  is  very 

fast  and  usually  only  three  or  four  iterations  are  required, 

A  flow  chart  for  this  subroutine  is  shown  in  Figures 

6a  and  6b. 

3.^.2  Input  List  for  Subroutine  UWDIM 

Ck LC(GG) 

C5 LC(GG) 

CP LC(GG) 

CX LG(GG) 

DISPL LC(GG) 

FRSC LC(GG) 

LEN LC(BB) 
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3. 4.3  Statement  Descriptions 

DPTRY  «  DISH/I .  014 

DISPL  has  the  value  assigned  to  DPTRY  in  other 
subroutines.   It  is  assumed  that  the  total  displacement 
is  1.4  percent  greater  than  the  molded  displacement. 
DPTRY  has  the  value  of  molded  displacement  in  this  sub- 
routine • 

CWP=.8065*CP+.2645 

This  relationship  is  plotted  in  Figure  8. 

CALPH».12*CW?-.04l4 

This  relationship  is  plotted  as  Figure  8-8(b)  in 
Reference  (1).  This  relationship  applies  to  destroyer 
ship  types  but  was  not  checked  against  Coast  Guard  designs. 

C1«DPTRY*35./(LE-N*CP*CX) 
C2=.833-CP*CX/(3.*CWP) 

C3»LEN*CA1PH/(DPTRY*35. ) 

C4«C4+FRSC 

The  derivation  of  these  coefficients  is  described 
in  the  introduction  to  this  section. 

BMIN=SQRT( (C5+SQRT(C5**2+12.*C1*C2*C3) )/(6.*C3) ) 

HMIN=C1/BMIN 

BTR=BMIN/HMIN 

IF(BTR.LT.2..0R.BTR.GT.4.)  GOTO  ! 

The  ratio  of  B/T  at  point  1  in  figure  7  is  calculated 
using  these  statements.   The  formula  for  BMIN  is  derived 
by  setting  dR/dE  ■  0  and  solving  for  B. 
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R=C2*HMIN+C3*BMIN**3-C4-C5*BMIN 

IF(R.LE.O)  GOTO  2 

C^«*C^-FRSC+R 

B1»BMIN 

H-HMIN 

RETURN 

These  statements  compute  the  value  of  R  if  point 

1  lies  within  the  acceptable  B/T  limits.  If  R  is  positive, 

the  program  returns  the  value  of  beam  and  draft  at  point 

1  and  increases  the  estimate  of  KG  by  R. 

1  HMAX»SQRT(Cl/2.005) 
BMIN=C1/HMAX 

R=C2*HMAX+C3*EMIN**3-C^-C5*BMIN 
IF(R.LE.O)    GOTO  2 
C4-C4-FRSC+R 

Bl-BMIN 

KaHMAX 

RETURN 

These  values  are  calculated  for  point  2  in  Figure  7. 

Actually  they  are   slightly  off  of  point  2  to  insure  an 

acceptable  value   in  the  speed-power  routine. 

2  HMIN=SQRT(Cl/3.995) 
BMAX-C1/HMIN 

R*C2*HMIN+C3*BMX**3-CJ+-C5*BKAX 
IF(R.LT.O.)    GOTO  6 

These  values  are  calculated  for  point  4  in  Figure   7« 

If  R  is  negative  at  this  point,  an  error  message  must 

be  printed. 

A(l)=l. 

A(2)=-CVC2 

A(3)=-C1*C5/C2 

A(4)=0. 

A(5)»C3/C2*C1**3 

B(l)-1. 

C(l)-1. 

Z-HMIN 

3  DO  k  I  =  2,5 
B(I)«A(I)+Z*B(I-1) 

*  C(I)=B(l)+Z*C(I-l) 
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The  columns  of  the  Birge-Vieta  method  are  calculated 

in  these  equations, 

ZNEW=Z-B(5)/C(^) 

IF(ABS(ZNEW-Z).LE..01*ZNEW)   GOTO  5 
Z«ZN2W 
GOTO  3 

A  revised  estimate  of  draft  is  calculated  and  the 

iteration  is  repeated  unless  the  difference  between  the 

old  and  new  values  is  less  than  1  percent  of  the  new 

estimate. 

5  H-ZNEW 
BlsCl/H 
C^Clf-FRSC 
RbO. 
RETURN 

The  results  are  returned  for  point  3  in  Figure  7. 

6  LEN  »  -1000. 
RETURN 

END 

These  statements  signal  subroutine  XECUTE  to  print 
an  error  message  and  proceed  with  a  new  data  case. 

3.^.4  Output  List  for  Subroutine  UWDIM 

Bl LC(GG) 

CWP LC(GG) 

H LC(GG) 

3.^.5  Nomenclature  List 

All  variables  have  the  same  definition  as  given 
in  the  MAIN  program  nomenclature  except  for  the  following i 
A(5)  column  1  of  the  Birge-Vieta  method 

B(5)  column  2  of  the  Birge-Vieta  method 


Nomenclature  List  (cont) 


7* 


Bl 
BMAX 

BMIN 
BTR 

C(5) 

CI 

C2 

C3 
Ck 

CALPH 

DISPL 

KMAX 

HMIN 

I 

R 

Z 

ZNEW 


Beam,  same  as  B  in  MAIN 

dummy  value  of  beam 

dummy  value  of  beam 

beam  to  draft  ratio 

column  3  in  of  the  Birge-Vieta  method 

coefficient,  see  text 

coefficient,  see  text 

coefficient,  see  text 

KG  or  KG  +  FRSC,  feet 

GMMUL  in  MAIN 

coefficient,  see  text 

DPTRY  in  MAIN 

dummy  value  of  draft 

dummy  value   of  draft 

index  variable 

KB  +  BM  -  GM  T  KG  -  FRSC 

dummy  value  of  draft 

dummy  value  of  draft 
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3.5  Subroutine  HPCALC 

3.5,1  Introduction 

Subroutine  HPCALC  performs  a  Taylor's  Standard 
Series  power  estimation. w '   The  procedure  used  is  very 
similar  to  the  hand  calculating  procedure.   However,  two 
problems  are  of  interest  with  regard  to  doing  the  calculations 
on  a  computer.  First,  the  Cr  curves  plotted  in  Reference 
(7)  must  be  converted  to  array  form  and  stored  in  the 
program.  Subroutine  CRVAL,  described  in  Section  3.6, 
chooses  the  correct  value  of  Cr  from  this  stored  array. 
Second,  the  surface  area  of  the  hull  must  be  estimated. 

One  approximate  formula  for  wetted  surface  is$ 

S  =  1.7*L«T  +  V/T 

where  L  is  the  ship  length 
T  is  the  draft 

7  is  the  volume  of  displacement 
This  form  was  used  as  a  starting  point  and  an  empirical 
relationship  was  derived  to  bring  the  wetted  surface 
more  closely  in  line  with  that  specified  on  the  curves 
given  in  Reference  (7).   To  do  this  a  number  of  points 
were  chosen  on  the  curves  given.   These  points  were 
compared  to  the  above  formula  and  some  correcting  factors 
for  B/T,  Cv  and  Cp  were  added  to  reduce  the  difference. 
The  formula  given  below  gives  values  within  5  percent 
of  the  values  calculated  from  the  curves  of  Reference  (7) 
over  the  entire  range. 

S  -  (1.7'L«T+V/T)(.0053(B/T)2-.02(B/T)+3#Cv+.08'CD+.926) 


76 

This  formula  is  probably  not  the  best  method  for 

estimating  the  surface  area.   One  method  would  be  to 
convert  the  curves  to  arrays  and  interpolate  to  find  the 
correct  area.  Such  a  method  requires  much  more  time  and 
compute  storage  to  produce  a  result  only  slightly  better 
than  this  formula.  Most  of  the  error  between  a  power 
calculated  by  this  routine  and  by  hand  is  probably  due 
to  error  in  the  surface  area  prediction* 

A  flow  chart  for  this  routine  is  given  in  Figure  9» 

3*5.2  Input  List  for  HPCALC  Subroutine 

B LC(GG) 

CP LC(GG) 

DELCF LC(CC) 

DPTRY LC(GG) 

H LC(GG) 

LEN LC(BB) 

PC SCDA 

VK SCDA 

3»5»3  Statement  Descriptions 

DISPL=DPTRY/l.0l4 

DISPL  is  the  molded  displacement.   It  is  assumed 
that  the  total  displacement  is  1.4  percent  greater  than 
the  molded  displacement. 

VOL«DISPL*35. 

BDR=B/H 

RN-VK*LEN*1 32050. 

CF».075/(AL0G10(RN)-2. )**2+DELCF 

SLRAT=VK/(SQRT(L2N)) 

CV-Y0I/(LEN)**3 


HPCALC   SUBROUTINE 
FLOW    CHART 
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These  formulas  calculate,  respectively,  volume  of 

displacement,  beam  to  draft  ratio,  Reynold's  number, 

Cf ,  speed  length  ratio,  and  Cv.  Cf  is  the  coefficient 

of  frictional  resistance  calculated  based  on  the  ITTC  line, 

IF(BDR.LT.2..0R.BDR.GT,4.)  GOTO  1000 
IF(CP.LT..48.QR.CP.GT.,7)  GOTO  1000 
IF(SLRAT.LT..5.0R.SLRAT.GT.2.)  GOTO  1000 
IF(CV.LT..001.0R.CV.GT..006)  GOTO  1000 
IF(CV,GT..003.AND.SLRAT.GT.1.3)  GOTO  1000 

These  statements  check  to  make  sure  none  of  the 
constraints  of  the  Cr  array  are  violated.   If  any  is 
violated,  the  program  will  write  an  error  message  and 
go  to  the  next  input  case, 

CALL  CRVAL(BDR,CV,CP,SLRAT,CR) 

This  subroutine  is  called  to  determine  the  value 

CR,  the  residuary  resistance  coefficient.   The  routine 

is  described  in  section  3«6. 

CT»CF+CR*.001 
Sa(1.7»I£N*H+V0L?H)*(.0053*BDR»BDR-.02*BDR+3.*CV+ 

,08*CP+.926) 
A«. 00438*1, 9905*S 
EHPBH«A*VK**3*CT 

The  total  resistance  coefficient,  wetted  surface 
and  finally  the  bare  hull  effective  horsepower  are  cal- 
culated in  these  statements, 

EHPAPP= ( 1 .467- . 2*SLRAT )*EHPBH 

The  effective  horsepower  with  appendages  has  been 

based  on  data  for  the  378*  WHEC,   This  data  is  shown 

in  Figure  10  together  with  the  above  approximation, 

EHP-.9*EHPAPP 
IF(SLRAT.LT.1.4,AND.SLRAT.GT..9)EHP«(1.74-.6*SLRAT)» 

EHPAPP 
IF(SLRAT.LE..9)  EHP»1.2*EHPAPP 
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dHJ 

dH3 
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These  equations  approximate  the  EHP/EHPipayior  curve 

for  the  378*  WHEC  as  shown  in  Figure  11. 

SHP-EHP/PC 
RETURN 

Finally,  the  shaft  horsepower  is  calculated  using 

the  input  propulsive  coefficient. 

1000  LEN»  -1000. 
RETURN 
END 

If  a  constraint  of  the  Cr  array  is  violated  the 

value  of  LEN  is  made  negative  to  signal  subroutine  XECUTE 

to  write  an  error  Eessage. 

3.5^  Output  for  Subroutine  HPCALC 

The  only  output  is  SHFS,  which  is  returned  as  a 
call  statement  dummy  argument. 

3.5»5  Nomenclature  List 

All  variables  have  the  same  definition  as  given  in 

the  MAIN  program  nomenclature  except  for  the  followingi 

A  dummy  coefficient 

BDR  beam  to  draft  ratio 

CF  frictional  resistance  coefficient 

CR  residuary  resistance  coefficient 

CT  total  resistance  coefficient 

CY  Taylor's  volumetric  coefficient 

DISPL  molded  displacement 

EHP  effeotive  horsepower 

EHPAPP  effective  horsepower  with  appendages 

EHPBH  effeotive  horsepower  bare  hull 
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-r--r 

.  i.  .  i. ...    ._ 
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jo|Abx 


dH3 


dH3 


Nomenclature  List  (cont) 

PC  propulsive  coefficient 

R£  same  as  RGEND  in  MAIN 

RN  Reynold's  number 

S  wetted  surface,  sq  ft 

SFCH  same  as  SFCHHP  in  MAIN 

SFCM  '  same  as  ^FCMHP  in  15AIN 

SHPS  shaft  horsepower 

SLRAT  speed  length  ratio 

VE  same  as  VEND  in  MAIN 

VK  speed  in  knots 

VOL  volume  of  displacement 
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3.6  Subroutine  CRVAL 
3.6,1  Introduction 

Although  this  routine  is  physically  one  of  the 
longest  in  the  cutter  model,  the  job  it  performs  is  one 
of  the  simplest.  Its  only  function  is  to  choose  from 
stored  data  arrays  the  values  of  the  residuary  resistance 
coefficient,  Cr»  which  straddle  the  input  point  and 
then  to  interpolate  to  find  the  value  of  Cr  at  the  input 
point.  The  input  point  is  specified  by  four  variables. 
The  value  of  each  of  these  four-beam  to  draft  ratio, 
Cp,  Cv,  and  speed  length  ratio-is  specif iod  as  input 
to  the  routine. 

The  constants  stored  in  the  three  arrays  which 
contain  the  Cr  data  were  determined  using  the  curves 
given  in  Reference  (7).  Although  the  three  arrays  - 
CR1 ,  CR2,  and  CR3  -  are  each  one  dimensional,  the  data 
they  contain  is  four  dimensional,  varying  with  each  of 
the  four  variables  which  specify  the  input  point.  Each 
of  the  arrays  contains  values  for  beam  to  draft  ratios 
of  2.25,  3.0,  and  3«75«  For  each  beam  to  draft  ratio 
Cp  was  varied  from  0.^8  to  0.70  in  steps  of  0.02.  For 
each  of  these  values  of  B/T  and  Cp,  all  the  possible 
combinations  of  Cv  and  speed  length  ratio  given  below 
were  stored. 

In  array  CR1  is  stored  data  for  Cv  values  of  0.001 
and  0.002  and  speed  length  ratios  of  0.5  to  2.0  in  steps 
of  0.1.   Array  CR2  corresponds  to  Cv  values  of  0.003, 
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0.004,  and  0.005  plus  speed  length  ratios  from  0.5  to 

1.3  in  steps   of  0.1.     Data  for  a  Cv  value   of  0.006  and 
speed  length  ratios   of  0.5  to  1.0  in  steps   of  0.1   is 
stored  in  CR3# 

In  each  of  the   three  arrays,   the  first  element 
corresponds  to  the  lowest  value   of  each  of  the  four  input 
variables.     Speed  lengcn  ratio  is  increased  most  rapidly 
followed  by  Cp,  C     and  finally  beam  to  draft  ratio. 

The  limits  of  the  Cr  arrays  given  above  are  extended 
somewhat  by  extrapolation.     Beam  to  draft  ratios  between 
2.0  and  4.0  are  accepted.      The  limiting  value   on  Cv  has 
been  extended  by  0.001  for  speed  length  ratios  between 
1.0  and  2.0.      The  maximum  ranges  that  are  accepted  are 
given  below. 

2.0  <  B/T  <  4.0 

-1*8  <  cp  <  *70 
.001  <  Cv  <  .006  for  0.5  <  V/Vl  <  1.3 
.001  <  Cv  <  .003  for  0.5  <  V/Vl  <  2.0 
Sixteen  values  of  Cr  are  required  to  straddle  the 
input  values  of  B/T,  Cp,  Cv,  and  V/^L.   After  these      s 

sixteen  numbers  have  been  chosen,  the  value  of  C„  at  the 

r 

input  point  is  found  by  interpolation.  A  total  of  fifteen 
interpolations  are  required,  eight  to  select  the  correct 
speed  length  ratio,  four  to  select  the  correct  C  ,  two 
to  select  the  correct  beam  to  draft  ratio,  and  a  final 
one  to  select  the  correct  Cv. 

Even  with  the  extensions  listed  earlier,  the  constraints 
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on  the  combined  Cr  array  cause  the  major  limitations 

of  the  cutter  model.   The  Cr  array  data  is  listed  in 
Appendix  D.  A  flow  chart  for  this  routine  is  shown  in 
Figures  12a  and  12b, 

3.6.2  Input  List  for  Subroutine  CRVAL 
BTR. SCDA 

CP SCDA 

CRK1152) LC(TT) 

CR2(972) LC(TT) 

CR3(2l6) LC(TT) 

CV SCDA 

VL SCDA 

3.6.3  Statement  Descriptions 

Il(J,K,L.M)a(K-l)*384+(J-l)*192+(L-l)*l6+M 

I2(J,K,L.M)«(K-1)*32MJ-3)*108+(L-1)*94-M 
I3(J.K,L,M)«(K-1)*72+(L-1)«6+M 

These  are  statement  functions  which  convert  the 

four  dimensions  -  J,  K,   L,  and  M  -  into  a  single  dimension 

for  use  with  arrays  CR1,  CR2,   and  CR3  respectively. 

IF(BTR.GE.3.)  Kl=2 
IF(BTR.LT,3.)  Kl«l 
K2.K1+1 

An  index  is  chosen  in  these  statements , for  the 

beam  to  draft  ratio  below  the  input  point  and  above  the 

input  point. 

CP1-.68 
3  IF(CP-CP1+.0001)  1,2,2 

1  CFUCP1-.02 
GOTO  3 

2  L1«IFIX((CP1-,46)*50.1) 
L2-L1+1 
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In  these  statements,  an  array  index  is  chosen  for  the 

C  below  and  above  the  input  point.   The  extra  0.0001 

in  statement  3  and  the  50.1  instead  of  50  in  statement 

2  are  used  to  account  for  small  errors  in  the  computer 

calculations. 

CVla.005 
I*   IF(CV-CV1+. 00001)  5.6,6 

5  cvi*cvi-.ooi 

GOTO  Jf 
13  Jl-Jl-1 

CV1-CV1-.001 
GOTO  Ik 

6  J1»IFIX(CV1»1001.) 
IF(Jl.EQ.5.AND.VL.GT.l.)  GOTO  13 

IF(J1.EQ.2.AND.VL.GT.1.3)  GOTO  13 
Ik   J2=J1+1 

Array  indices  for  Cv  are  chosen  below  and  above 

the  input  point.   The  statements  connected  with  statement 

13  extend  the  limits  of  the  array  by  extrapolation. 

VL1-1.9 

7  IF(VL-VL1+.0001)  8,9,9 

8  VL1-VL1-.1 
GOTO  7 

9  M1=IFIX((VL1-.4)*10.1) 
M2=M1+1 

Again  array  indices  are  chosen,  this  time  for 

speed  length  ratio  values. 

J-Jl 

INDEX-1 

The  value  of  Cv  is  held  fixed  each  time  through 

the  following  statements.   When  C  values  at  each  of  the 

two  Cv  values  have  been  calculated,  an  interpolation 

is  done  to  find  the  Cr  value  at  the  correct  Cv  value. 

11  IF(J.G2.3)  GOTO  10 
I«I1(J,K1,L1,M1) 
CRA=CR1(I) 

. 
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I=I1(J,K2,L2,M2) 
CRH*CR1(I) 
GOTO  50 

Variables  CRA  to  CRH  are  used  for  the  eight  Cr  values 

that  remain  when  Cv  is  held  constant.  In  this  set  of 

statements  array  CR1  is  used* 

10  IF(J.EQ.6)   GOTO  12 
I-I2(J,K1,L1,M1) 

CRA«CR2(I) 


I»I2(J,K2,L2,M2) 
CRK«CR2(I) 
GOTO  50 

This  group  of  statements  is  identical  to  the  proceeding 

group  except  that  array  CR2  is  used. 

12  I»I3(J.K1,L1,M1) 

CRA«CR3(D 

# 

I»I3(j.K2,L2,M2) 
CRH»CR3(I) 

Again  these  statements  are  identical  with  array 

CR3  being  used, 

50  V=.(VL-VL1)*10. 
C»(CP-CP1)*50. 
BTR1«FL0AT((K1-1))». ?5+2,25 
Ba(BTR-BTRl)/. 75 
C1«(CRB-CRA)*V+CRA 
C2»(CRD-CRC)*V+CRC 
C3«(CRF-CRE)»V+CRS 
C4=(CRH-CRG)»V+CRG 
C5*(C2-C1)*C+C1 
C6c:(C^-C3)*C+C3 

C7-(c6-c5)*bk;5 
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These  are  the   interpolation  statements  required 

to  first  calculate   the  value   of  Cr  at  the  input  speed 

length  ratio,   then  at  the  input  C_  and  finally  at  the 

input  beam  to  draft  ratio. 

IF(INDEX.GE.2)    GOTO  51 

A=C7 
J=J2 
INDEX-2 

GOTO  11 

When  IND2X  ■  1   the  lower  value   of  Cv  is  used  in  the 

interpolation  above.      The  next  step  is  to  do  the  interpolation 

at  the  upper  value   of  Cy.     The  value   of  Cr  at  the  lower 

Cv  is  saved  as  variable  A. 

51  D=C? 

CRa(D-A)*(CV-CVl)«1000.+A 

RETURN 
END 

One  final  interpolation  is  performed  to  calculate 

the  desired  value  of  CR  and  the  routine  returns  control 

to  subroutine  HPCALC. 

3.6.^  Output  for  Subroutine  CRVAL 

The  only  output  is  the  calculated  value  of  CR  which 
is  returned  as  a  subroutine  call  dummy  argument. 

3.6.5  Nomenclature  List 

All  variables  have  the  same  definition  as  given  in 
the  MAIN  program  nomenclature  except  for  the  followingi 
A  Cr  value  at  lower  Cv  value 

B  fraction  of  BTR  interval 

BTR  input  value  of  beam  to  draft  ratio 

BTR1  lower  value  of  beam  to  draft  ratio 
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C  fraction  of  C  interval 

Cl  Cr  at  correct  V/^L,  lower  BTR,  lower  C 

C2  Cr  at  correct  V/JL,   lower  BTR,  higher  Cp 

C3  Cr  at  correct  V/Vl,  higher  BTR,  lower  C_ 

Ck  Cr  at  correct  V/JZ,   higher  BTR,  higher  C 

C5  Cr  at  correct  V/^L  and  C  ,  lower  BTR 

C6  Cr  at  correct  V/JZ  and  C  ,  higher  BTR 

C7  Cr  at  correct  V/^U,  Cp  and  BTR 

CP1  lower  C  value 

P 

CR  Cr  value  at  input  point 

CRA  Cr  at  lower  V/^L,  Cp  and  BTR 

CRB  Cr  at  higher  V/JLt   lower  Cp  and  BTR 

CRC  C  at  lower  V/^L  and  BTR,  higher  Cp 

CRD  Cr  at  higher  V/JZ   and  C  ,  lower  BTR 

CRE  C  at  lower  V/^L  and  C  ,  higher  BTR 

CRF  Cr  at  higher  V/JZ   and  BTR,  lower  C 

CRG  Cr  at  lower  V/JZ,   higher  BTR  and  C 

CRH  Cr  at  higher  V//L,  BTR,  and  Cp 

CV  input  value  of  Cy 

CY1  lower  value  of  Cv 

D  Cr  value  at  upper  Cy  value 

I  array  index 

INDEX  Indicates  whether  Cv  is  at  lower  or  upper  value 

J  used  for  either  Jl  or  J2 

Jl  index  for  lower  value  of  Cv 

J2  index  for  higher  value  of  Cv 

Kl  index  for  lower  value  of  BTR 


K2  index  for  higher  value  of  BTR 

LI  index  for  lower  value  of  C 

L2  index  for  higher  value  of  C 

Ml  index  for  lower  value  of  V/Vl 

M2  index  for  higher  value  of  V/JZ 

V  fraction  of  V/^L  interval 

YL  input  speed  length  ratio 

YL1  lower  value  of  speed  length  ratio 
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3.7  Subroutine  EPLANT 

3.7.1  Introduction 

This  routine  performs  two  functions.   The  first  is 
to  estimate  an  average  plant  electrical  load  for  use 
in  calculating  fuel  requirements.   The  second  is  to 
size  the  generators.  Generator  size  is  used  primarily 
for  estimating  weights,  but  it  is  also  a  useful  output 
which  provides  information  needed  for  beginning  the 
preliminary  design  phase. 

Three  ships  were  used  in  developing  electrical 
load  relations.   These  were  the  378'  WHEC,  the  210*  WMEC, 
and  the  150'  WPB.   The  electrical  loads  for  these  vessels 
were  categorized  as  shown  in  Table  1.  This  table  also 
gives  the  loads  in  KW  that  were  used. 

The  first  two  categories,  electronics  load  and  arm- 
ament load,  are  inputs.  Propulsion  auxiliary  load  includes 
only  equipment  needed  for  the  direct  operation  of  the 
main  propulsion  machinery.  For  example,  lub  oil  service 
pumps  are  included  in  this  category  but  lub  oil  transfer 
pumps  are  included  under  other  auxiliaries.  Air  condition- 
ing and  ventilation  also  includes  electric  heaters. 
Hotel  is  primarily  galley  and  lighting. 

Estimating  relationships  were  developed  for  each 
of  the  last  four  categories  in  Table  1.   If  the  propulsion 
auxiliary  load  is  not  input,  a  value  of  1  KW  is  assumed 
for  diesel  plants  and  71  KW  for  C0D0G  plants.  Air 
conditioning  and  ventilation  is  calculated  as  a  function 
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Table 

1   1 

ELECTRICAL  LOAD  REQUIREMENTS 

Electrical  Load 

150*   WPB 

210*   WMEC 

378'    WHBC 

Electronics 

4.8 

21.99 

117.84 

Armament 

42.7 

0.0 

46.64 

Propulsion  Aux. 

1.02 

.71 

71.03 

Air  Cond.   &  Yent. 

21.1 

^9.13 

308.34 

Hotel 

20.9 

35.4 

170.51 

Steering  gear 
Aux.   Machinery 
Shops 
Other  Aux. (total) 

.5^             ^8              37.38 

18.26            40.08              70.72 

1.0                2.2                  4.98 

19.8               46.76           113.08 

Total  110.32     153.99     827.44 

All  loads  are  in  KW 


of  cubic  number  times  the  number  of  accommodations.  Hotel 
load  is  a  function  of  the  number  of  accommodations  only. 
Other  auxiliaries  load  is  assumed  to  be  a  function  of 
cubic  number  only.   These  relationships  are  plotted  in 
Figures  14,  15  and  16. 

The  total  electrical  load  required  is  equal  to  the 
sum  of  the  six  categories  above  with  a  two  percent  growth 
margin  on  all  equipment  plus  a  100  percent  margin  on 
electronics.  An  additional  25  percent  is  included  above 
the  growth  margin  to  allow  for  deterioration  and  starting 
loads. 
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It  is  assumed  that  three  generators  will  be  installed* 

Each  generator  is  sized  so  as  to  be  able  to  carry  the  full 
electrical  load*   This  represents  the  current  Coast  Guard 
practice  for  new  designs.  Less  capacity  has  been  installed 
on  past  designs. 

The  smallest  generator  that  will  be  installed  is 
assumed  to  be  100  KW.   Other  possible  plant  sizes  lie 
between  100  KW  and  250  KW  in  steps  of  50  KW,  between 
250  KW  and  1000  KW  in  steps  of  250  KW,  and  above  1000  KW 
in  steps  of  500  KY/.   The  smallest  generator  that  is  able 
to  supply  the  total  load  with  the  margins  stated  above  is 
chosen. 

A  flow  chart  for  this  routine  is  shown  in  Figure  1?« 

3.7.2  Input  List  for  Subroutine  EPLANT 

BL0AD LC(DD) 

CN LC(BB) 

EL0AD LC  (DD) 

J0PT LC(CC) 

MTYPE LC(CC) 

NCP0 LC(AA) 

NENL LC(AA) 

N0FF LC(AA) 

PAL0AD LC(DD) 
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3.7«3  Statement  Descriptions 

IF(J0PT.EQ.2)  GOTO  1 

The  estimate  of  PALOAD  is  skipped  if  a  value  is 

given  as  input. 

PAL0AD=1. 

IF(MTYPE.EQ.2)  PAL0AD«71. 

Propulsion  auxiliaries  electrical  load  is  calculated 

for  diesel  and  CODOG  plants  respectively, 

1  OAL0AD=.0002695*CN 
NAC=NOFF+NCPO+NENL 
HL0AD=.917*NAC 

AYLOAD=15.  +  .376*NAC*''CN/100000. 

KW»EL0AIH-BL0AlM-PAL0AD*0AL0AD+HL0AD+AVL0AD 

AVGKW=KW-BLOAD 

KW1=1.25*(1.02*KW+EL0AD) 

The  individual  load  items  are  calculated  and  summed 

to  give  an  estimate  of  the  total  load.   The  average  load 

is  taken  as  the  total  load  minus  the  armament  load* 

Margins  are  then  added  to  the  total  load  for  sizing  the 

generators. 

ELKW=100. 

2  IP(ELKW.GE.KWl)  GOTO  3 
IF(ELKW.LT.2^0.)  GOT  k 
IF(ELKW.LT.900.)  GOTO  5 
ELKW-ELKW+500. 

GOTO  2 

4  ELKWsELKW+50. 
GOTO  2 

5  ELKWaELKW+250. 
GOTO  2 

3  RETURN 
END 

These  statements  choose  the  generator  rating,  ELKW, 

which  is  the  next  size  larger  than  the  total  load  with 

margins • 
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3.7.^  Output  List  for  Subroutine  EPLANT 

AVGKW LC  (DD) 

ELKW LC  (DD) 

3.7.5  Nomenclature  List 

All  variables  have  the  same  definition  as  given 
in  the  MAIN  program  nomenclature  except  for  the  following i 


AVLOAD 

HLOAD 

KW 

KW1 

NAC 

OALOAD 

RE 

SFCH 

SFCM 

VE 


air  conditioning  and  ventilation  load,  KW 

hotel  load,  KW 

total  load  without  margins,  KW 

total  load  with  margins,  KW 

number  of  accommodations 

other  auxiliaries  load,  KW 

same  as  RGEND  in  MAIN 

same  as  SPCHHP  in  MAIN 

same  as  SFCMHP  in  MAIN 

same  as  VEND  in  MAIN 
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3.8  Subroutine  MACHLQ 

3,8,1  Introduction 

At  this  point  in  the  program,  the  shaft  horsepower 
and  size  of  the  electrical  plant  have  been  determined  and 
the  weight  of  fuel,  lub  oil  and  potable  water  can  be 
calculated.  Subroutine  MACHLQ  performs  these  calculations,' 
In  conjunction  with  the  fuel  weight  calculations,  the 
specific  fuel  consumption,  SFC,  at  endurance  power  is 
estimated. 

For  diesel  plants  it  is  assumed  that  the  SFC  at 
rated  power  is  0.^2  lbs/shp-hr  and  at  one  half  rated 
power  is  0.46  lbs/shp-hr.  A  linear  relation  is  used 
for  other  power  ratios.   The  SFC  used  for  C0D0G  plants 
is  0.42  if  the  endurance  power  is  less  than  7000  SHP. 
This  assumes  that  the  diesels  will  be  sized  to  be  at 
maximum  power  at  endurance  speed.  If  the  endurance 
power  requirements  are  greater  than  7000  SHP,  the  gas 
turbines  must  be  used  at  endurance  speed.   Three  maximum 
power  levels  for  the  gas  turbines  are  used  -  25000  HP, 
35000  HP  and  50000  HP.   At  each  of  these  horsepower 
levels  the  SFC  is  taken  to  be  0.50  lbs/shp-hr.   At  half 
power,  the  SFC  is  .595*  A  linear  relation  is  used  for 
other  power  levels. 

If  a  specific  machinery  plant  is  input,  the  specific 
fuel  consumptions  at  half  and  full  power  are  given.  A 
linear  relation  is  used  for  other  power  levels. 
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The  specific  fuel  consumption  of  the  ship  service 

generators  is  assumed  to  be  about  0.48  lbs/shp-hr  or 

0.65  lbs/KW-hr.   The  generator  fuel  consumption  is  computed 

using  the  average  electrical  load. 

Fuel  consumption  for  hotel  services  is  taken  to  be 
0.32  lbs/man-hr.   This  assumes  some  waste  heat  recovery. 

The  lub  oil  requirements  for  diesels  is  assumed  to 
be  5  tons  plus  1  ton  per  thousand  horsepower.  For  C0D0G 
plants  this  same  relationship  is  used  for  the  diesel  part 
of  the  plant  (endurance  power  is  used  in  computing  the 
requirements).   However,  additional  lub  oil  for  the  gas 
turbines  is  included.  For  the  gas  turbines,  the  oil 
requirements  are  3.5  tons  plus  1  ton  per  ten  thousand 
horsepower. 

Potable  water  requirements  are  taken  to  be  50  gal/man 
or  0.186  tons/man.  A  flow  chart  for  this  routine  is 
shown  in  Figure  18. 

3.8.2  Input  List  for  Subroutine  HACHLQ 

AVGKW LC  ( DD )     RE LC  ( CC  ) 

JOPT LC  (CC  )     SFCH LC  (CC  ) 

MTYFE LC  (CC  )  SFCM LC(CC  ) 

NCPO LC(AA)  SHEE LC(CC) 

NENL LC(AA)  SHFK. . . . LC(CC) 

NOFF LC(AA)  VE LC(CC) 
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3,8,3     Statement  Descriptions 

NAC«N0FF+NCPO+NENL 

The  total  number  of  accommodations  is  calculated* 

WTPS=.l86*NAC 

Additional  items  will  be  added  to  the  weight  of 
personnel  stores  later  but  for  now  it  includes  only  the 
weight  of  potable  water, 

IF(J0PT.EQ«2)  GOTO  1 

This  skips  the  lub  oil  and  standard  SFC  estimates 
since  these  values  are  inputs, 

WTLO=.001*SHR5+5. 

IF(MTYPE.EQ,2)   WTL0».0001*SHPZS+.001*SHF£+8.5 

The  lub  oil  weight  for  diesel  and  CODOG  plants 

respectively  are  calculated  here, 

IF(MTYFE.EQ.2)   GOTO  2 

SFC= .  5-.  08*SHE3/SHPM 

SHPH-SHPM/2. 

IF(SHK4.GT.7000..^D.SHPE.LT.SHPH)SFC=.5-.l6«SHFE/SHFM 

GOTO  3 

These  statements  calculate  the  SFC  for  diesel  plants. 

In  case  the  shaft  horsepower  is  greater  than  7000  it  is 

assumed  that  four  diesels  will  be  installed.   If  the 

endurance  power  is  less  than  half  of  the  maximum,  only 

two  diesels  need  be  in  operation  at  endurance  speed, 

2  IF(SHFE.LT,?000,)  GOTO  k 

SFC=(1.38-.38*SHPE/50000.)».5 

IF(SHPM.LE. 25000,  )SFO(1.38-.33*SHFE/25000.  )«,5 

IF(SHPM. GT. 25000., AND.SHFM.LE. 35000,)  SFC-U.38- 

,38«SHPE/35000.)*,5 
GOTO  3 
If  SFC=,42 
GOTO  3 

For  CODOG  plants  the  SFC  is  calculated  using  the 
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statements  above, 

1  SFC=SFCM+(1.-SHPE/SHI^)*2.*(SFCH-SFCM) 

If  SFCH  and  SFCM  are  input,  the  above  linear  relation 

is  used  for  calculating  the  SFC  at  endurance  speed. 

3  FRMAIN=SFC*SH?E 
PRH0TL«.32*NAC 

FRGEN=«.  65  *AVGKW 
FRcFRMAI N+FRKOTL+FRGEN 

All  of  these  fuel  rates  are  in  lbs/hr. 

WTFUEL=FR*RE/(1788.*VE) 

RETURN 

END 

The  factor  1788.   is  derived  froxai 

22frOT»(l   -  FtpJL 
rm  rdeg  rsc 


where  1 


F^  ■  .05  is  a  tailpipe  allowance 

Fm  «  1.03  is  an  extra  margin 

F^e_*  1.05  is  a  degradation  allowance 

F8C  *  1.10  is  an  allowance  for  sea  conditions 

and  fouling 
The  value  2240  converts  pounds  to  tons. 

3.8.4  Output  List  for  Subroutine  MACHLQ 

WTFUEL LC  (LL) 

WTLO LC(LL) 

WTPS LC(LL) 

3«8.5  Nomenclature  List 

All  variables  have  the  same  definition  as  given  in 
the  MAIN  program  nomenclature  except  for  the  followingi 


Nomenclature  List  (cont) 
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FR 

FRGEN 

FRHOTL 

FRMAIN 

NAC 

RE 

SFC 

SFCH 

SFCM 

SHPH 

VE 


total  fuel  rate  in  lbs/hr 

generator  fuel  rate  in  lbs/hr 

hotel  fuel  rate  in  lbs/hr 

main  engine  fuel  rate  in  lbs/hr 

total  number  of  accommodations 

same  as  RGfEND  in  MAIN 

specific  fuel  consumption  at  endurance  power 

same  as  SFCHHP  in  MAIN 

same  as  SFCMHP  in  MAIN 

one  half  maximum  shaft  horsepower 

same  as  VEND  in  MAIN 
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3.9  Subroutine  HBSIZE 

3,9,1  Introduction 

Very  little  data  is  available  concerning  the  minimum 
allowable  depth  and  length  of  the  machinery  compartment. 
As  a  result,  only  very  simple  calculations  are  made  in 
this  subroutine  to  estimate  these  two  dimensions. 

The  product  of  length  of  machinery  box  times  the 
maximum  beam  of  the  ship  was  calculated  for  two  diesel 
designs,  the  210*  WMEC  and  the  150*  WPB,  and  for  one 
C0D0G  plant,  the  378*  WHEC,   A  value  of  about  1050  sq. 
ft.  was  obtained  for  both  diesel  designs  and  a  value  of 
2300  sq,  ft,  was  obtained  using  the  C0D0G  plant.   These 
values  are  divided  by  the  beam  of  the  ship  to  obtain 
the  required  engine  room  length  in  the  program, 

A  check  is  also  made  to  insure  that  a  minimum  length 
is  not  exceeded.   This  minimum  length  was  taken  to  be 
28  feet  or  56  feet  for  diesel  plants  below  or  above 
7000  BHP  respectively  and  45  feet  for  a  C0D0G  plant, 

A  rough  deck  height  approximation  was  derived  using 
data  from  the  same  three  vessels  mentioned  above.  The 
formula  used  gives  a  deck  height  of  about  17  feet  for 
the  150*  WPB  and  about  19  feet  for  the  210'  WMEC.   The 
deck  height  for  all  plants  with  horsepowers  greater  than 
7000  is  taken  as  22  feet, 

A  flow  chart  for  this  subroutine  is  shown  in  Figure 
19. 


MBSIZE      SUBROUTINE 
FLOW         CHART 
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3.9«2  Input  List  for  HBSIZE  Subroutine 

5 ..SCDA 

MTYPE LC(CC) 

SHPM LC(CC) 

3.9.3  Statement  Descriptions 

IF(MTYPE.EQ.2)  GOTO  1 

XLM»1050./B 

IP(SHFil.GT. 7000,. AND.XLM.LT. 56.)  XLM«56. 

IF(SHFM.LE«7000..AND.XLM.LT.28.)  ZLM-28. 

GOTO  2 

These  statements  set  the  engine  room  length  for 

diesel  plants.  (MTYPE  «  1  for  diesels) 

1  XLM=»2300./B 
IF(XLM.LT.^5«)  XLM=^5» 

These  statements  set  the  engine  room  length  for 
CODOG  plants. 

2  DHM=11.+.00157*SHPM 
IF(DKM.GT.22.)  DKM=22. 
RETURN 

END 

The  minimum  deck  height^  which  is  the  same  for  both 
plant  types,  is  set  by  these  statements. 

3.9.^  Output  List  for  Subroutine  MBSIZE 

All  variables  have  the  same  definition  as  given  in 
the  MAIN  program  nomenclature  except  for  the  following* 
RE  same  as  RGEND  in  MAIN 

SFCH  same  as  SFCHHP  in  MAIN 

SFCM  same  as  SFCMHP  in  MAIN 

VE  same  as  VEND  in  MAIN 
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3.10  Subroutine  VOLUME 
3,10.1  Introduction 

Many  of  the  ships  designed  today  are  volume  rather 
than  weight  limited.   Because  of  this  it  is  important 
that  a  ship  synthesis  model  determine  a  balance  between 
required  and  available  volume  as  well  as  between  weight 
and  displacement.  Subroutine  VOLUME  performs  this  function 
in  the  cutter  model. 

The  dimensions  of  the  hull  below  the  load  waterline 
are  determined  in  subroutine  UWDIM.   Only  the  sheer  line, 
deckhouse  volume  and  the  length  of  raised  deck,  if  any, 
remain  to  be  determined  to  define  the  volume  available 
in  the  ship.   The  enclosed  volume  of  the  hull  and  deck- 
house will  be  completely  determined  by  these  variables 
since  the  waterplane  coefficient  i6  calculated  in  sub- 
routine UWDIM  and  since  the  general  shape  of  the  hull 
is  implied  by  the  estimating  relationships  used  in  the 
model. 

Deckhouse  volume  is  constrained  primarily  by  stability 
considerations  and  by  external  deck  area  requirements. 
If  the  deckhouse  becomes  too  high,  it  will  not  be  possible 
to  make  the  vessel  stable.  Also  a  certain  amount  of 
free  deck  area  is  required  in  nearly  all  designs,  restricting 
the  horizontal  spreading  of  the  deckhouse.  At  the  other 
extreme,  the  ship  may  have  excessive  stability  if  too 
small  a  deckhouse  is  installed. 

These  deckhouse  volume  restrictions  have  been  included 


112 

in  the  cutter  model  by  restricting  deckhouse  size  to  be 

within  the  limits  of  past  practice.  Both  an  upper  and 
a  lower  limit  are  placed  on  the  size  of  the  deckhouse 
as  shown  in  Figure  20.  The  deckhouse  volume  of  Coast 
Guard  vessels  has  been  generally  greater  than  that  of 

Navy  designs.  Upper  and  lower  limits  of  0,00150  x  Ir  and 

3 
0,001  x  L  ,  respectively,  are  used  in  the  Navy's  destroyer 

model.  This  difference  is  probably  due  to  the  fact  that 

the  Navy  requires  more  external  deck  area  for  armament* 

The  deckhouse  volumes  shown  in  Figure  20  do  not 
include  uptakes.  It  is  assumed  that  sufficient  additional 
area  is  available  for  uptakes  since  uptake  area  was 
excluded  from  the  data  points  also.  However,  it  is 
important  to  check  a  design  be  sketching  a  layout  for  the 
external  decks  and  the  superstructure  before  acceptance. 

The  volume  in  the  hull  is  also  eatimated  in  VOLUME. 
The  underwater  volume  is  already  available.  The  first 
step  in  determining  the  above  water  hull  volume  is  the 
development  of  an  acceptable  sheer  line.  This  step  is 
performed  by  subroutine  SHEER  and  is  discussed  in  con- 
junction with  that  routine. 

Once  the  sheer  line  and  with  it  the  average  freeboard 
have  been  estimated,  the  above  water  volume  is  found 
by  multiplying  the  area  of  the  waterplane  by  the  average 
freeboard  and  by  a  factor  to  account  for  flare.  The 

total  volume  in  the  hull  is  the  sum  of  the  above  and  below 

water  volumes. 
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This  total  hull  volume  includes  the  volume  to  the 
main  deck  only.   The  hull  volume  cannot  be  reduced  below 
this  size  but  a  raised  deck  can  be  added  to  increase  the 
volume  of  the  hull.  For  vessels  of  the  size  of  Coast 
Guard  cutters  it  is  impractical  to  add  more  than  one 
raised  deck.  Also  the  deckhouse  size  can  be  varied  some- 
what to  change  the  available  volume. 

Required  hull  volume  is  broken  down  into  three  cat- 
egories.  These  are  the  machinery  box,  tankage  volume 
and  arrangements  volume.   The  object  of  subroutine  VOLUME 
is  to  alter  the  available  volume  so  that  there  is  sufficient 
space  for  each  of  the  three  categories  of  required  volume. 
The  model  assumes  that  all  tankage  volume  and  the  machinery 
box  volume  are  in  the  hull.  The  deckhouse  contains  only 
arrangements  volume. 

The  machinery  box  size  is  the  first  category  estimated* 
This  volume  is  found  by  multiplying  the  midships  area 
by  the  length  of  the  machinery  box  and  by  a  machinery 
box  prismatic  coefficient.  Machinery  box  volume  is 
subtracted  from  the  available  hull  volume  to  the  main 
deck  leaving  a  volume  which  can  be  used  for  tankage  and 
arrangements. 

Some  of  the  remaining  volume  cannot  be  used  for 
arrangements  because  of  hull  shape.   The  next  step  is  to 
estimate  the  fraction  of  the  remaining  volume  that  can 
be  used  only  for  tankage*   This  is  the  minimum  tankage 
volume  for  the  ship.  If  less  tankage  is  required  for 
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the  liquids  that  must  be  carried,  the  excess  volume  must 

be  used  for  voids  or  cofferdams  since  it  cannot  be  con- 
verted into  machinery  box  space  or  arrangements  volume* 
If  more  tankage  is  required  than  this  minimum,  some  of 
the  arrangements  volume  must  be  used  for  tankage. 

The  weights  of  all  the  liquids  which  must  be  carried 
are  calculated  in  subroutine  MACHLQ.   In  VOLUME,  these 
weights  are  converted  to  volumes  and  added  together. 
A  twenty  percent  allowance  is  added  to  this  figure  to 
account  for  required  cofferdams  and  voids.   The  volume 
required  for  liquids  is  then  compared  to  the  minimum 
tankage  volume  to  determine  the  amount  of  arrangements 
volume,  if  any,  that  is  required  for  tankage. 

All  of  the  volume  remaining  in  the  hull  after  tankage 
volume  and  machinery  box  volume  have  been  subtracted , can  be 
used  for  arrangements.  The  entire  superstructure  volume 
is  also  available  for  arrangements.   Arrangements  space 
requirements  are  more  commonly  given  as  areas  rather  than 
as  volumes,  so  the  next  step  is  to  convert  the  available 
arrangements  volume  to  an  available  arrangements  area. 
For  the  deckhouse  this  is  a  simple  task  since  the  deck- 
house is  nearly  wall  sided.   The  volume  is  divided  by 
the  deckheight  of  8.5  feet  to  find  the  area.  Por  the 
hull  no  such  simple  procedure  is  possible  because  of 
hull  curvature. 

The  Navy  performed  extensive  work  on  this  problem 
while  developing  their  destroyer  model,  DD07.   This  work 
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was  made  use  of  in  the  cutter  model  not  only  for  the 

arrangements  area  formulation  but  also  for  the  flare 
factor,  machinery  box  prismatic  coefficient  and  tankage 
volume  fraction  estimates  discussed  earlier.   The  relation- 
ships used  in  the  cutter  model  were  taken  directly  from 
the  original  notes'  ^  of  Mr,  Jim  Mills  who  was  largely 
responsible  for  the  development  work  done  by  the  Navy, 

The  major  deviation  between  the  cutter  model  and  the 
destroyer  model  lies  in  the  engine  room  location.   The 
destroyer  model  assumes  that  the  engine  room  is  centered 
at  amidships  but  Coast  Guard  practice  usually  places  the 
engine  room  aft  of  amidships.   No  correction  is  made  for 
this  difference  and  the  output  produced  for  test  runs 
with  existing  ships  indicates  that  there  is  little  effect. 

By  using  the  Navy's  formula,  the  arrangements  area 
in  the  hull  can  be  calculated.   This  area  plus  the  area 
of  the  deckhouse  must  be  compared  to  the  area  required 
for  all  the  ship's  functions  plus  the  area  specified  as 
input.   The  area  required  for  ship's  functions  such  as 
stores  area,  living  area  and  auxiliary  machinery  area, 
are  estimated  in  the  cutter  model  based  on  past  practice. 

Data  was  obtained  by  measuring  the  areas  of  five 
past  designs  and  assigning  the  areas  to  various  categories. 
This  data  is  summarized  in  Table  2,   Linear  relationships 
were  developed  from  this  data  as  shown  in  Figures  21 
through  25.   Berthing  and  sanitary  areas  are  based  on 
minimum  Coast  Guard  standards.   The  area  required  for 
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CATEGORY 


Office  Spaces 
iv.es sing  Facilities 
Crew  Special 
Officer  Staterooms 
Officer  Sanitary 
CPO  Staterooms 
CPO  Sanitary 
Crew  Berthing 
Crew  Sanitary 
CO  Stateroom  &  Cabin 
CO  Pantry 
Commissary  stores 
Other  Stores 
Workshops 
Passages 
Repair  Lockers 
Steering  Gear  &  Windlass 
Chain  Locker 
A/C  &  Fan  Spaces 
I.C.  &.   Gyro  Room 
Aux.  Machinery  Spaces 
Uptakes  (Hull) 
Pilothouse,  Chartroom  & 

C.I.C. 

Sub  total 


378 

327 

255 

210 

150 

WHEC 

WHEC 

WHEC 

WMEC 

WPB 

<*93 

465 

347 

130 

65 

2713 

2348 

1379 

1452 

339 

725 

609 

717 

324 

32 

1157 

1338 

966 

56? 

187 

243 

133 

120 

193 

30 

1040 

364 

336 

283 

110 

250 

100 

- 

75 

27 

3641 

2655 

2568 

1658 

375 

716 

444 

492 

384 

90 

697 

511 

273 

406 

- 

58 

128 

22 

- 

— 

600 

724 

531 

480 

35 

3015 

235^ 

2397 

1294 

435 

1161 

818 

752 

544 

- 

2567 

2588 

2138 

739 

229 

148 

18 

97 

110 

- 

736 

548 

645 

270 

294 

47 

30 

36 

42 

16 

369 

324 

39 

165 

- 

312 

143 

232 

180 

- 

i860 

1874 

414 

820 

543 

462 

1216 

64 

- 

1125 

856 

521 

422 

300 

24135 

20588 

15086 

10538 

3107 

EXTRAS 


Crew  Recreation  Rooms 
Cargo  &  Flume  Tanks 
Electronics  Sc   Stores 
Armament  &   Stores 
3alloon  Shelters 
Mission  Offices  &   Labs 
Mission  Workshops 
Bow  Thruster 

Total 


1445 

- 

682 

- 

71 

1219 

249 

- 

- 

241 

1618 

775 

1206 

319 

157 

844 

822 

730 

280 

- 

512 

156 

120 

- 

- 

746 

336 

482 

- 

- 

- 

- 

70 

105 

- 

344 

~ 

~ 

™ 

■• 

30863 

22926 

I8376 

11242 

3576 

Table 
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repair  lockers  is  taken  as  a  constant  100  square  feet. 

The  area  for  chain  lockers  as  40  square  feet  and  the  area 
for  uptakes  in  the  hull  as  either  6$   or  460  square  feet 
depending  on  whether  diesels  or  gas  turbines,  respectively, 
are  installed.  A  constant  value  of  nine  percent  of  the 
total  area  is  assumed  for  passages.   Passages  ranged 
from  6.84  to  13»17  percent  of  the  total  area  on  the 
five  designs  used. 

All  of  the  required  areas  are  summed  to  determine 
the  total  required  area.   This  area  is  then  compared  to 
the  available  arrangements  area  in  the  hull  and  super- 
structure. 

Comparisons  are  made  in  the  following  order.  Pirst, 
the  arrangements  area  required  to  be  located  in  the 
deckhouse  is  compared  to  the  area  of  the  largest  allowable 
deckhouse.  If  the  required  area  exceeds  that  available, 
the  design  is  infeasible  and  an  error  message  must  be 
printed.  Next,  the  total  required  arrangements  axea 
is  compared  to  the  total  available  area.   If  the  required 
area  is  the  smaller,  the  deckhouse  size  is  decreased. 
The  deckhouse  size  used  is  either  the  size  which  balances 
the  area  requirements  or  the  smallest  allowable  deckhouse 
size  or  the  area  required  to  be  in  the  deckhouse,  whichever 
is  greatest.  Should  the  total  required  area  exceed 
that  available  a  raised  deck  is  added. 

The  length  of  raised  deck  is  determined  by  an  iterative 
procedure.   On  each  iteration  the  area  required  to  be 
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included  in  the  raised  deck  is  computed  and  the  curve 

shown  in  Figure  26  is  used  to  convert  the  area  to  a 
raised  deck  length.  This  curve  is  based  on  the  main 
deck  of  a  378'  WHEC.   The  iteration  continues  until  two 
successive  iterations  fall  within  5  percent  of  the  length 
of  the  ship  or  until  twenty  iterations  is  reached.  In 
the  later  case  the  design  is  ruled  infeasible.  A  check 
is  also  made  to  insure  that  the  raised  deck  length  does- 
not  exceed  the  length  of  the  ship. 

If  the  raised  deck  length  lies  between  0«*  and  6.6 
times  the  length  of  the  ship,  the  raised  deck  is  extended 
to  0.6  times  the  length  of  the  ship. 

A  flow  chart  for  this  routine  is  shown  in  Figures 
2?a  and  27b. 

3.10.2  Input  List  for  Subroutine  VOLUME 

AREADH. LC(BB)  LEN LC(BB) 

AREAHL LC(BB)  LMB LC(JJ) 

B LC(GG)  MTYEE LC(CC) 

CP LC  ( GG )  NCPO LC  (  AA) 

CWP LC  (GG)  NENL LC(AA) 

CX LC(GG)  NOFF LC(AA) 

D LC(AA)  WACFUL LC(LL) 

DMB LC(JJ)  WTFDEL LC(LL) 

H LC(GG)  WTLO LC(LL) 

JOPT LC (CC  )  WTPS LC(LL) 
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3.10.3     Statement  Descriptions 

KNDEX-0 
NAON0FF+NCP0+NENL 

DSCKHT=0. 

LRD=*0. 

LRD1-0. 

:     KNDEX  is  an  iteration  counter  used  for  the  raised 
aec*  length  iteration.  The  variable  DECKHT  is  explained 
ln  the  SHEER  subroutine.  The  two  values  of  raised  dec* 

length  start  at  zero. 

13  CALL  SHBER(DHB.FO.FlO.F2O.DO.DlO.02O.aH.H.WJKHf) 

The  sheer  line   of  the  vessel     is  defined  by  the 
freeboard  at  stations  0,   10.   and  20.     These  values  are 
returned  as  F0.  F!0.  and  F20  respectively.     The  associated 
hull  depths  are  also  returned. 

AREA-LEN/6 .  *  (F0+F20+i+*Fl  0 ) 

favg-areaAsn 

DAVG*FAVG*H 
!  The  average   depth  is  calculated  using  the  projected 

lateral  area.     A  better  estimation  could  be  made  by  using 
the   local  beam  to  weight  the  freeboard  values  but  this 
was  not  done. because   the  value   of  F7  which  corrects  for 

flare   is  based  on  the   above  average  depth. 

F7-.018828MDAVG/H)**2+.18098*DAVG/H+.71^599 

This  is  the  flare  correction  factor.   The  formula 

r^m  the  Navv's  work  in  developing  their 
used  was  taken  from  the  Navy  s 

destroyer  model. 

HVAW«LEN*B»CW?*FAVG*F7 
HVBW«LEN*B*H*CP*CX 

TKV-HVAV/+HVBW 

1  The  total  hull  volume  to  the  »ain  deck  is  calculated 
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by  these  statements. 

AK«B*H*CX+B*F10 

X-LMB/LEN 

CPM— .  525*X*X+ .  055*X+ .  99675 

CPM=CPM+3.2*X*X*(CP-.6) 

VMB»LMB*AM*CPM 

AHVOM-THV-VMB 

The  volume   of  the  machinery  box  and  the  remaining 

available  hull  volume  are  calculated  in  these  statements. 

The  formula  for  CPM  ha*  been  taken  from  the  Navy's  work 

and  applies  to  machinery  boxes  centered  about  amidships* 

Q«(-.Q037115+.2139?27*D10-1.38263E-2*D10**2+ 
4.008058B-4*D10**3-5.^89^81E-6*D10***H- 
2.892153E-8*D10**5)*CP/.58«CV.8i>« 

(l.ll*LMB/LEN+.667) 

ATVOM»Q*AHVOM 
VFUEL»WTPUEL*^3. 

VLO-WTLO*39« 

VFW».186*NAC*36. 

YACFUL»WACPUL*^4 . 1 

VREQ«(VFUEL+YL0+VFW+YACPUL)*1.2 

RTAA»0. 

IF(VREQ.GT.ATVOM)   RTAA«(VREQ-ATV0M)/8.5 

AAVOM-AHVOM-ATVOM 

TANKVL-ATVOMf RTAA*8. 5 

The  available  volume  which  can  be  used  only  for 

tankage  and  the  required  tankage  volume  are  calculated 

and  compared  in  these   statements.     The  formula  for  Q 

was  again  taken  from  the  destroyer  model  work.     The 

required  tankage  arrangements  area  is  calculated  if  the 

required  tankage  volume   is  in  excess  of  that  available. 

P«(6.51 04lE-7*Dl  0**3-7 .  851484B-5*D1  0**2+3.  438375S-3* 
D10+. 0530399 )*(.191*CP+. 8393 )*(«191*CX+. 8445)* 

(-.105*U1B/LEK+1.0315)-  .  .  ':■  * 

AVAA-AAV 0M*P-RTAA 

The  available  arrangements  volume  is  converted  to 

an  available  arrangements  area  again  using  a  Navy  developed 

formula. 
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DHV«.00223*LEN**3 

DHAA-DHV/8.5 

CN= ( DAVG+LRD/LEN*8 . 5 ) *LEN«B 

RDHA»(AREADH+2^.  +  .00225«Cr*+2O7*LEN)*1.09 

IP(RDHA.GT.DHAA)  GOTO  1000 

TAAA=DHAA+AVAA 

The  area  in  the  largest  acceptable  deckhouse  is 
compared  to  the  area  required  to  be  in  the  deckhouse.' 
If  the  required  area  is  smaller,  the  program  continues. 
If  not,  an  error  messa^j  is  printed.   The  areas  required 
to  be  in  the  deckhouse  include  the  input  areas  plus  the 
commanding  officers  cabin  and  the  pilothouse.   An  allowance 
is  also  made  for  passageways. 

0DA«2.12*LEN-232. 
. 

PHDA«100.+.00225*CN 

The  statements  above  and  those  included  between 
them  are  used  to  calculate   the  required  deck  areas. 

TRAA=1 .  09*(  AREADH+AREAHI^0DA+XM?DA+CSDA+0SRDA+0HDA+ 
C  PODA+C  POHDA+CBDA+C  HDA+C  0DA+C3  TDA+  CSDA+WDA+ 
RLDA+CLDA+UTDA+SGDA+ACDA+XICDA+AHDA+PHDA) 

IP(TRAA.GT.TAAA)    GOTO  10 

The  total  required  arrangements  area  is  computed 

and  compared  with  the  area  available  in  the  hull  to  the 

main  deck  plus  the  largest  acceptable  deckhouse.  A 

decision  is  then  made  to  either  make  the  deckhouse  smaller 

or  to  add  a  raised  deck. 

DHSV-. 00151 *LEN»*3 
DHSAA=DKSV/8.5 
DHAMIN-AMAX1(DKSAA,RDHA) 
TAAA-DHAMIH+AVAA 
IF(TRAA.GT.TAAA)  GOTO  11 
12  DHV«DKSV 

DHAA«DHY/8.5 
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ENCVOL=DHV+THV 

RETURN 
11   DHAR=TRAA-AVAA 
DHSY=DHAR*8.5 
GOTO  12 

These  statements  compute   the  deckhouse  size  required 
if  the  deckhouse   is  to  be  made  smaller.     First,   the 
larger  value  of  the  minimum  deckhouse  size   or  the  required 
deckhouse  area  is  chosen  and  added  to  the  hull  voulume. 
This  new  available  area  is  then  compared  to  the  total 
required  area.      If  the  available  area  is  still  the  greater t 
the  difference   is  excess  area  that  cannot  be  reduced 
since  the  minimum  size  hull  and  deckhouse  are  already- 
being  used.     If  more  than  the  minimum  deckhouse  size  is 
required  a  calculation  is  made   to  balance   the  total 
required  area  and  the   total  available  area. 

10  RKAA=TRAA-DHAA 
RDAA=RHAA-AVAA 
X1=RDAA/(B*LEN) 

IF(X1.GT. .7)   GOTO  1 
IF(X1.GT..2)   GOTO  2 
LRD=-3.22»X1*X1+2.031*X1-.018^ 
GOTO  3 

1  LRD=.463*Xl«Xl+.^6*Xl+.22 
GOTO  3 

2  LRD-X1+.059 

3  LRD=xLRD*LEN 
DECKKT«8.5 

IF(ABS(LRD-LRD1).LT..05*LEN)  GOTO  15 
IF(KND2X.GT.20)  GOTO  1000 
KNDSX-KND2X+1 

LRD1 =LRD 
GOTO  13 

If  more  area  is  required  than  is  available  in  the 

largest  deckhouse  plus  the  hull  to  the  main  deck,  a  raised 

deck  is  added.   The  assumption  is  made  that  a  ship  with 

the  maximum  deckhouse  size  will  result  in  the  least 
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displacement.   This  is  probably  valid  so  long  as  an 

aluminum  deckhouse  is  used. 

The  length  of  raised  deck  required  is  computed  in 

the  above  statements.   The  new  raised  deck  length  is  then 

compared  to  the  former  estimate,  LRD1.   If  the  two  do 

not  agree  within  five  percent  of  the  length  of  the  ship 

a  new  iteration  is  made   A  default  counter  has  been 

included  to  keep  the  computer  from  being  caught  in  a 

nonconverging  loop. 

15  C1«.^*LEN 
C2=.6*LEN 

IF(LRD.LT.C1.0R.LRD.GT.C2)    GOTO  If 
AVAA-AVAA+ ( C2-LRD) *B 
RDAA-RDAA+ ( C 2 -LRD ) *B 
LRD=*C2 

These  statements  check  to  see  if  the  raised  deck 

length  falls  within  0.4  to  0.6  times  the  ship  length. 

If  so  the   calculated  values  are  corrected  so  that  the 

raised  deck  length  io  equal  to  0.6  times   the  ship  length. 

4  IP(LRD.GT.LEN)   GOTO   1000 
THV=THV+RDAA*8.5 
ENCV01>DHV+THV 
RETURN 

A  check  is  first  made   to  insure   that  a  raised  deck 

length  longer  than  the  length  of  the  ship  is  not  required. 

A  corrected  total  volume   is  then  calculated  and  the   program 

returns. 

1000  LEN=-1000. 
RETURN 
END 

Statement  1000  is  a  default  option  warning  subroutine 

XECUTE  to  print  an  error  message. 


3.10.^  Output  List  for  Subroutine  VOLUME 
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ACDA LC(PP) 

AMDA LC(PP) 

AVAA LC(PP) 

CBDA LC(PP) 

CHDA LC(PP) 

CLDA LC(PP) 

CN LC(BB) 

CODA LC(PP) 

CPODA LC(PP) 

CPOHDA LC(PP) 

CSDA LC(PP) 

CSTDA LC(PP) 

DO LC(II) 

D10 LC(II) 

D20 LC(II) 

DAVG LC(BB) 

DHAA LC(PP) 


DHV LC(BB) 

ENCVOL LC(BB) 

LRD LC(JJ) 

ODA. LC(PP) 

OHDA LC(PP) 

OSDA LC(PP) 

OSRDA LC(PP) 

PHDA LC(PP) 

RLDA LC(PP) 

SGDA LC(PP) 

TANKVL LC(PP) 

UTDA LC(PP) 

VMB LC(PP) 

WDA LC(PP) 

XICDA LC(PP) 

XMFDA LC(PP) 


3.10.5  Nomenclature  List 

The  definition  of  all  variables  is  the  same  as 
given  in  the  MAIN  program  nomenclature  except  for  the 
following 1 

AAVOM         available  arrangements  volume  in  hull,  cu  ft 
ACDA  A/C  &  ventilation  deck  area,  sq  ft 

AHVOM         available  hull  volume,  cu  ft 
AM  cross  sectional  area  amidships,  sq  ft 

AMDA  aux.  machinery  spaces  deck  area,  sq  ft 

AREA  projected  lateral  area  above  water,  sq  ft 


ATVOM 
AVAA 

CI 
C2 

CBDA 

CHDA 

CLDA 

CODA 

CPM 

CPODA 

CPOHDA 

CSDA 

CSTDA 

D 

DO 

DIO 

D20 

DECKHT 

DHAA 

DHAMIN 

DHAR 

DHSAA 

DHSY 

DMB 

FO 

P7 

PIO 


available  tankage  volume,  cu  ft 

available  arrangements  area  in  hull,  sq  ft 

0,4  x  length,  ft 

0,6  x  length,  ft 

crews  berthing  deck  area,  sq  ft 

crews  heads  deck  area,  sq  ft 

chain  locker  deck  area,  sq  ft 

CO.  S.R.,  Cabin  &  pantry  deck  area,  sq  ft 

machinery  box  prismatic  coefficient 

CPO  S.R.  deck  area,  sq  ft 

CPO  heads  deck  area,  sq  ft 

crew  special  deck  area,  sq  ft 

commissary  stores  deck  area,  sq  ft 

same  as  EKDDAY  in  MAIN 

depth  at  P.P.,  ft 

depth  at  amidships,  ft 

depth  at  A. P. ,  ft 

height  of  raised  deck,  ft 

deckhouse  arrangements  area,  sq  ft 

minimum  deckhouse  area,  sq  ft 

deckhouse  area  required,  sq  ft 

smallest  deckhouse  arrangements  area,  sq  ft 

smallest  deckhouse  volume,  cu  ft 

same  as  DHM  in  MAIN 

freeboard  at  P.P.,  ft 

flare  factor 

freeboard  at  amidships,  ft 


F20 

FAYG 

HVAW 

HVBW 

KNDEX 

LMB 

LRD1 

NAC 

ODA 

OHDA 

OSDA 

OSRDA 

P 

PHDA 

Q 

RDAA 

RDHA 

RE 

RHAA 

RLDA 

RTAA 

SFCH 

SFCM 

SGDA 

TAAA 

TANKVL 
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freeboard  at  A. P. ,ft 

average  freeboard,  ft 

hull  volume  above  waterline ,  cu  f  t 

hull  volume  below  waterline,  cu  ft 

counter 

same  as  XLH  in  MAIN 

used  to  store  LRD  value 

number  of  accommodations. 

office  deck  area,  sq  ft 

officer  heads  deck  area,  sq  ft 

other  stores  deck  area,  sq  ft 

officer  staterooms  deck  area,  sq  ft 

coefficient 

pilothouse,  chartroom  &  CIC  deck  area,  sq  ft 

coefficient 

raised  dock  arrangements  area,  sq  ft 

area  required  to  be  in  deckhouse,  sq  ft 

same  as  RGBND  in  MAIN 

required  hull  arrangements  area,  sq  ft 

repair  locker  deck  area,  sq  ft 

required  tankage  arrangements  area,  sq  ft 

same  as  SFCHHP  in  MAIN 

same  as  SFCMHP  in  MAIN 

steering  gear  &   windlass  deck  area,  sq  ft 

total  available  arrangements  area,  sq  ft 

tankage  volume ,  cu  ft 


THV 

TRAA 

UTDA 

VACFUL 

VE 

VFUEL 

VFW 

YLO 

VMB 

VREQ 

WDA 

X 

XI 

XICDA 

XMFDA 


total  hull  volume,  cu  ft 

total  required  arrangements  area,  sq  ft 

uptake  deck  area  in  hull,  sq  ft 

volume  of  aircraft  fuel,  cu  ft 

same  as  VEND  in  MAIN 

volume  of  ships  fuel,  cu  ft 

volume  potable  water,  cu  ft 

volume  of  lub  oil,  cu  ft 

volume  of  machinery  box,  cu  ft 

required  tankage  volume,  cu  ft 

workshops  deck  area,  sq  ft 

LMB/LEN 

RDAA/(B*LEN) 

I.e.  &   gyro  room  deck  area,  sq  ft 

messing  facilities  deck  area,  sq  ft 
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3.11  Subroutine  SHEER 

3.11 #1  Introduction 

Subroutine  SKSBR  selects  an  appropriate  sheer  line 
by  specifying  the  freeboard  at  the  forward  perpendicular, 
amidships,  and  at  the  after  perpendicular.   This  sheer 
line  must  satisfy  a  number  of  criteria. 

First,  there  must  be  sufficient  depth  amidships 
to  accommodate  the  engine  room.   The  depth  amidships 
must  also  be  sufficient  to  insure  adequate  structural 
strength*  A  minimum  value  of  L/16  is  used*   Because 
Coast  Guard  vessels  normally  have  engine  rooms  aft  of 
amidships,  a  slightly  greater  depth  amidships  is  used 
rather  than  the  machinery  box  depth.   This  accounts  for 
the  drop  in  the  sheer  line  aft  of  amidships. 

The  next  criterion  which  must  be  satisfied  is  adequate 
freeboard  at  the  forward  and  after  perpendiculars. 
Navy  derived  formulas  were  used.   These  formulas  give 
the  minimum  acceptable  freeboard  based  on  a  deck  wetness 
criterion. 

Once  independent  determinations  of  freeboard  at  the 
forward  and  after  perpendiculars  and  at  amidships  are 
made ,  a  check  is  required  to  insure  that  a  reasonable 
sheer  line  can  be  drawn  between  the  three.  Data  was 
obtained  from  past  designs  and  from  the  standard  merchant 
ship  sheer  curves.   This  data  was  nondimensionalized  by 
subtracting  the  freeboard  amidships  from  the  other  free- 
board values  and  then  dividing  by  the  length  of  the  ship. 
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The  resulting  curves  are  plotted  in  Figure  28, 

The  sheer  line  is  assumed  to  be  practical  if  the 
forward  sheer  fraction  lies  between  0,01  and  0,03  and  the 
after  sheer  fraction  lies  between  0.001  and  0.0075. 
If  this  is  not  the  case,  one  or  more  of  the  freeboard 
values  is  adjusted  until  the  criterion  is  satisfied. 

If  a  raised  deck  is  added  to  the  hullt  a  new  sheer 
line  is  calculated.   It  is  assumed  that  the  required 
freeboard  of  the  main  deck  at  the  forward  perpendicular 
can  bo  one   deck  height  or  8.5  feet  below  that  required 
by  the  deck  wetness  criterion.   This  will  result  in  a 
flattened  sheer  line  as  is  common  for  raised  deck  vessels. 

A  flow  chart  for  this  routine  is  shown  in  Figure  29  • 

3.11.2  Input  List  for  Subroutine  SHEER 

D&B SCDA        H .SCDA 

LEN SCDA        DECKHT SCDA 

3.11.3  Statement  Descriptions 

D10«.0025*LEN+DMB 

IF(D10.LT.LEH/l6)  D10-LEN/16. 

The  minimum  depth  amidships  is  selected  by  these 

statements. 

F0«1.01l827*(100.*H/LEN)-.000636215*LEN+2. 780649 
P0»U3N»F0/100. 

F0-F0-D2CKHT 
DO-FO+H 

The  minimum  depth  at  the  forward  perpendicular  is 

selected  here. 
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SHEER  /  LBP 


d8~l  /  U33HS 


SHEER    SUBROUTINE 
FLOW   CHART 
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1*1 

F20-.01*LEN*(2.125+.00125*LEN) 
D20«F20+H 

Finally  the  minimum  depth  at  the  after  perpendicular 

is  selected. 

7  S0=D0-D10 
RO«SO/LSN 

IF(RO.LT..Ol)   GOTO  1 
IF(R0.GT..03)   GOTO  2 

The  sheer  fraction  forward  is  calculated  and  checked 

against  the  acceptable  limits, 

6  S20=D20-D10 
R20=S20/LEiN 

IF(R20.LT..001)    GOTO   3 
IF(R20.GT..0075)    GOTO  k> 
GOTO  5 

The  sheer  fraction  at  the  after  perpendicular  is 

calculated  and  checked  against  acceptable  limits. 

1  S0=.01*LEN 
D0-SCKD10 
GOTO  6 

If  the  freeboard  at  the  bow  is  too  low  the  depth 

at  the  bow  is  increased. 

2  S0=.03*LEN 
D10=D0-S0 
GOTO  6 

If  the  freeboard  at  the  bow  is  excessive ,   the  depth 

amidships  is   increased. 

3  S20«.001*LEN 
D20*D1Q+S20 
GOTO  5 

If  the  freeboard  aft  is  too  low  the  depth  aft  is 

increased. 

k   S20=.0075*LEN 
D10=*D20-S20 
GOTO  7 
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If  the  freeboard  aft  is  excessive  the  depth  amidships 

is  increased  and  the  sheer  fraction  at  the  bow  is  rechecked. 

5  F0sD0-H 
F20«D20-H 
F10=.D10-H 
RETURN 
END 

Values  of  both  freeboard  and  depth  are  returned 

for  use  in  subroutine  /olume. 


3.11.^  Output  List  for  Subroutine  SHEER 

FO SCDA        DO SCDA 

F10 SCDA        D10 SCDA 

F20 SCDA        D20 SCDA 

3*11.5  Nomenclature  List 

DO  depth  at  forward  perpendicular ,  ft 

D10  depth  at  amidships,  ft 

D20  depth  at  after  perpendicular,  ft 

EECKHT  height  of  raised  deck,  ft 

DMB  machinery  box  depth,  ft 

FO  freeboard  at  bow,  ft 

F10  freeboard  amidships,  ft 

F20  freeboard  at  stern,  ft 

H  draft,  ft 

LEN  length  between  perpendiculars,   ft 

RO  (D0-D10)/LEN 

R20  (D20-D10)/LEN 

SO  D0-D10 

S20  D20-D10 


3.12  Subroutine  WEIGHT 
3.12.1  Introduction 

Subroutine  VOLUME  balanced  the  required  and  available 
volumes  in  the  hull.   This  subroutine  estimates  weights 
so  that  a  balance  can  be  made  between  full  load  weight 
and  displacement. 

Several  of  the  weights  required  have  already  been 
calculated  or  were  given  as  input.  Calculated  weights 
include  the  fuel  weight  and  lub  oil  weight.  All  com- 
munications and  control,  group  k,   weights  and  all  armament, 
group  7,  weights  are  given  as  input.  The  weight  groups 
noted  refer  to  the  NAVSHIPS  Hull  Group  Weight  Classification 
of  1965«   The  categories  in  this  classification  scheme 
are  given  in  Appendix  E.   This  is  the  classification 
system  now  used  by  the  Coast  Guard. 

Also  given  as  input  are  the  weights  of  cargo,  air- 
craft, ammunition  and  aircraft  fuel.  All  other  weights 
which  are  a  part  of  the  full  load  weight  are  estimated 
in  subroutine  WEIGHT.   A  flow  chart  for  this  routine  is 
. shown  in  Figure  30. 

There  are  few  patrol  cutters  of  recent  vintage 
which  can  be  used  for  collecting  weight  data.   Three 
have  been  used  throughout  this  model.   These  are  the 
378'  WHEC,  the  210'  WMEC,  and  the  150'  WPB.   Wherever 
possible  additional  data  has  been  shown.  Light  ship 
weight  data  was  taken  from  the  weight  breakdown  sheets 
prepared  for  each  design.  Load  weights  were  taken  from 
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1*5 

either  load  weight  calculations  or  from  inclining  experiment 

data* 

The  first  hull  weight  group  is  the  hull  structure, 
group  1.  In  order  to  estimate  this  weight  group  more 
accurately,  three  of  the  sub-groups  are  estimated  sep- 
arately.  The  remainder  of  the  sub-groups  are  lumped 
together  as  a  single  group.  In  the  program  this  group 
is  referred  to  as  WT1M,   The  three  groups  subtracted 
are  group  111,  superstructure!  group  112,  foundations 
for  propulsion  plant  machinery;  and  group  113»  foundations 
for  auxiliaries  and  other  equipment.   These  are  identified 
in  the  program  as  WT111,  WT112,  and  WT113»  respectively. 

Figure  31  is  a  plot  of  WT1H  versus  cubic  number. 
The  correlation  here  is  very  good.  Also  included  on  the 
figure  is  the  estimating  line  used  in  the  Navy's  destroyer 
model,  DD07. 

Superstructure  weight,  WT111,  is  shown  in  Figure  32. 
plotted  against  deckhouse  volume.   The  line  labeled  no 
uptakes  is  the  one  used  in  the  program.   The  line  used 
in  DD07  is  also  shown.   This  line  assumes  that  the  deck- 
house is  aluminum. 

The  weight  of  propulsion  plant  foundations,  WT112, 
is  shown  in  Figure  33  plotted  against  the  weight  of 
propulsion  machinery  minus  shafting  and  propellers.  Again, 
a  line  used  in  DDO?  is  shown, 

WT113,  other  foundations,  is  shown  in  Figure  3^ 
versus  the  total  weight  of  groups  3,  **,  5»  and  7*     The 
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correlation  here  appears  reasonable;  however,  there  io 

a  difference  between  the  data  and  the  DDO?  estimate. 

Propulsion  machinery  is  classified  under  group  2. 
In  this  program,  an  estimate  of  the  weight  of  sub-group 
203  (shafting,  bearings,  and  propellers)  was  made  separately. 
The  remaining  sub-groups  were  estimated  as  one  group,  WT2M« 
A  plot  of  WT2M  versus  shaft  horsepower  for  diesel  plants 
is  shown  in  Figure  35  •  Here  more  data  is  available  from 
non  patrol  types.   The  DD07  estimate  for  diesel  plants 
is  poor.  For  C0D0G  plants  it  is  assumed  that  the  diesel 
plant  will  be  7000  horsepower  or  under  and  that  only 
the  endurance  power  will  cause  a  difference  in  weight 
from  that  of  a  378*  WHEC.  A  lower  endurance  power  will 
reduce  the  weight  somewhat. 

Shafting  and  propeller  weight,  WT203,  is  calculated 
as  weight  per  foot  of  shaft  length  versus  shaft  horsepower. 
Twin  shafts  are  assumed  with  a  total  shaft  length  O.56 
times  the  length  of  the  ship.  A  plot  of  this  relationship 
is  shown  in  Figure  36. 

Group  3  includes  the  electrical  plant  weights.   The 
Coast  Guard  intends  to  install  electric  plants  of  relative- 
ly larger  capacity  in  future  ships  so  a  simple  linear 
estimation  using  current  designs  is  unrealistic.  Instead, 
the  weights  of  a  number  of  diesel  generating  sets  of 
various  capacities  have  been  plotted  as  shown  in  Figure  37* 
The  varying  RPM  of  these  plants  resulted  in  three  weight 
ranges. 
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The  weight  of  the  generating  seta  is  only  a  part 

of  the  group  3  weight.  The  weight  of  switchboards,  WT301, 
is  shown  in  Figure  38  as  is  the  remaining  group  3  weight, 
WT3R.   The  correlation  for  those  weight  groups  is  quite 
bad.  Fortunately,  these  weights  make  up  a  small  percentage 
of  the  total  weight  of  the  ship. 

Since  group  4  weights  are  given  as  input,  group  5 
weights,  auxiliary  syetems,  is  estimated  next.   The 
destroyer  model  uses  the  estimating  relation  shown  in 
Figure  39.   The  three  Coast  Guard  vessels  agree  very 
well  with  this  relationship  • 

Group  6,  outfit  and  furnishings  is  the  last  light 
ship  group  estimated  since  group  7  weights  are  given  as 
input.  Group  6  is  broken  into  four  sub-groups  for  estimating 
purposes.   The  first,  WT6l,  includes  groups  600,  602, 
603,  and  605.  These  are  hull  fittings,  rigging  and  canvas, 
ladders  and  gratings,  and  painting,  respectively.  The 
second  group,  WT611,  includes  groups  604,  nonstructural 
bulkheads  and  doors  1  606,  deck  covering*  607,  hull  insulation; 
and  608,  storerooms,  stowages  and  lockers.   Group  three, 
WT6III,  includes  only  group  601,  boats,  boat  stowage 
and  handling.   The  final  group,  WT6IV,  includes  all  the 
remaining  sub-groups  of  weight  group  6. 

WT61  plotted  against  enclosed  volume  is  shown  in 
Figure  40.  WT611  and  WT6IV  are  plotted  against  number  of 
accommodations  as  shown  in  Figure  kl.     WT6III  is  estimated 
in  three  parts.   The  first  is  a  boat  weight,  the  second  a 
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liferaft  weight  and  the  third  a  life  jacket  weight. 

All  of  the  lightship  weights  are  summed  and  then 
multiplied  by  a  design  and  builders  margin.   The  design 
and  builders  margin  is  an  input. 

Crew  weight  is  taken  to  be  0.0737  tons  per  man  or 
165  pounds  per  man.   Personal  effects  are  calculated 
as  0.105  tons  per  man  for  officers,  .0737  for  chiefs 
and  .029  for  enlisted  men. 

The  weight  of  personnel  stores  includes  provisions, 
general  stores,  and  potable  water.  The  weight  of  provisions 
is  based  on  Table  3*   The  data  in  this  table  was  reduced 
to  a  formula  for  which  only  the  endurance  days  of  dry 
provisions  need  be  given.   This  formula  isi 

(.0222+. 00202  x  D)  x  NAC 

where  D  is  the  endurance  days  for  dry  provisions 
NAC  is  the  number  of  accommodations 

Table  3 

days  endurance  pounds/person 

dry  prov.   frozen    chilled      total 

30  97.0  36.8  51.8  185.6 

45  145.5  55.2  52.8  253.5 

60  194.0  73.6  53.8  321.4 

90  291.0  110.4  55.8  457.2 

The  weight  of  general  stores  is  taken  to  be 
.00135  x  D  x  NAC 
and  that  of  potable  water  to  be 

•186  x  NAC  or  50  gallons/man 
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The  full  load  weight  can  then  be  calculated  by  adding 

the  load  weights  to  the  light  ship  weight. 
3.12.2  Input  List  for  Subroutine  WEIGHT 


SHEE LC(CC) 

SHKI LC(CC) 

WACFUL LC(LL) 

WCARGO LC  (LL) 

WTAC LC(LL) 

WTAMMO LC(LL) 

WTFUEL LC(LL) 

WTG2? LC(WW) 

WTG^f LC(WW) 

WTG7 LC(WW) 

WTLO ....LC(LL) 


CN LC(BB) 

0 LC(AA) 

DBMAR LC(WW) 

DHV LC(BB) 

ELKW LC(DD) 

ENCVOL LC(BB) 

JOPT LC(CC) 

LEN LC(BB) 

MTYEE ..LC(CC) 

NCPO LC(AA) 

NENL LC(AA) 

NOFF LC(AA) 

♦only  if  JOPT  «  2 

3.12.3  Statement  Descriptions 

A  detailed  description  of  each  statement  will  not 
be  given  for  this  routine  wince  the  flow  of  the  program 
follows  the  description  given  in  the  introduction  and 
the  relationships  are  those  plotted. 

One  area  that  will  be  discussed  is  WT6III.  Here 
the  number  of  15  man  liferafts  required  is  first  calculated. 

Q.NAC/15.+1. 

RAFTS=AINT(Q) 

where  AINT  is  a  function  which  drops  the  decimal  fraction 
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from  its  argument.  Boat  weight  is  assumed  to  be  1.4  tons 

up  to  25  men,  5  tons  between  25  and  50  men  and  10.1  tons 
above  50  men.  This  weight  includes  launching  equipment. 
Life  jacket  weight  is  taken  to  be  ,002  tons  per  man. 

3.12.4  Output  List  for  Subroutine  WEIGHT 


WFULLD LC(WW) 

WLSHIP LC(WW) 

WT111 LC(W) 

WT112 LC(YV) 

WT113 LC(W) 

WT1M LC(VV) 

WT203 LC(W) 

WT2M LC(W) 

WT300 LC(VV) 

WT301 LC(W) 

WT3R LC(W) 

WT6I LC(VV) 

*  only  if  J0PT  -  1 

3.12.5  Nomenclature  List 

The  definition  of  all  variables  is  the  same  as 
given  in  the  MAIN  program  nomenclature  except  for  the 
following! 

B0ATWT        boat  weight,  tons 
D  same  as  ENDDAY  in  MAIN 

NAC  number  of  accommodations 

Q  dummy  argument 


WT6II LC(VV) 

WT6III LC(W) 

WT6IV LC(YV) 

WTCREW LC(LL) 

WTG1 LC(WW) 

WTG2? LC(WW) 

WTG3 LC(WW) 

WTG5 LC(WW) 

WTG6 LC(VfW) 

WTPE LC(LL) 

WTFS LC(LL) 


RAFTS 
RE 

SFCH 

SFCM 

VE 

VOL 

WT111 

WT112 

WT113 
WT1M 

WT203 

WT2M 

WT300 

WT301 

WT3R 

WT6I 

WT6II 

WT6III 

WT6IV 

WTFDS 


number  of  liferafts 

same  as  RGEND  in  MAIN 

same  as  SFCHHP  in  MAIN 

same  as  SFCMHP  in  MAIN 

same  as  VEND  in  MAIN 

enclosed  volume  x  10""^ 

weight  of  group  lilt  tons 

weight  of  group  112,  tons 

weight  of  group  11 3 ,  tons 

weight  of  remainder  of  group  1,  tons 

weight  of  group  203 »  tons 

weight  of  remainder  of  group  2,  tons 

weight  of  group  300,  tons 

weight  of  group  301»  tons 

weight  of  remainder  of  group  3,   tons 

see  text 

see  text 

see  text 

see  text 

sum  of  weight  groups  3,  kt   5,  and  7 
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3.13  Subroutine  VCO 

3,13,1  Introduction 

The  vertical  center  of  gravity  of  the  ship  is  estimated 
in  this  routine.   The  designer  has  some  latitude  with 
regard  to  the  location  of  the  center  of  gravity 1  there- 
fore, the  value  calculated  by  this  routine  should  be  viewed 
as  a  feasible  rather  tnan  as  the  best  value. 

The  effect  that  the  veg  has  on  the  design  can  be 
investigated  by  using  different  values  of  free  surface 
correction.   The  free  surface  correction  is  added  to  the 
KG.   Therefore,  either  a  positive  or  a  negative  value 
can  be  used  to  represent  a  shift  in  the  KG  location. 

In  order  to  determine  feasible  locations  for  the 
various  centers,  data  for  the  378 •  WHEC,  the  210*  WMEC, 
and  the  150'  WPB  was  used.   This  data  is  shown  in  Tables 
k   and  5»  On  these  tables  the  vegs  of  the  various  weight 
groups  are  shown.   These  centers  are  related  to  various 
hull  depths  by  a  factor  K,  The  hull  depth  used  is  given 
in  the  right  hand  column  of  these  tables.   An  average 
coefficient,  K,  was  chosen  and  is  given  in  the  column 
labeled  Coef .  K. 

The  variables  to  which  the  vegs  are  referred  were 
chosen  arbitrarily  by  observing  which  gave  the  most 
consistent  set  of  coefficients  among  the  three  ships. 

Some  of  the  weight  groups  shown  in  the  tables  are 
really  input  items.   These  weight  groups  are  included 
to  assist  the  designer  in  choosing  realistic  values  to 
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VERTICAL 

CENTERS 

OF   GRAVITY 

Weight 

373' 
WHEC 

210* 

WMEC 

150* 
WPB 

Coef . 
K 

times 

tfTlw 

vcg 
K 

15.73 
.553 

14.72 
.547 

9.55 
.534 

.55 

DAVG 

Willi 

vcg 
K 

37.1 

1.4 

35.3 
1.26 

23.6 
1.47 

1.4 

D10 

WT112 

vcg 
K 

5.73 
.216 

4.4 
.157 

2.5 

.156 

.2 

D10 

WT1 1 3 

vcg 
K 

21.05 
.794 

16.6 
.593 

15.1 

.944 

.87 

Center  groups 
3+4+5+7 

WT2I.; 

vcg 
K 

15.14 
.571 

10.40 
.533 

8.50 
.531 

.54 

D10MN 

WT203 

vcg 
K 

4.02 
.152 

6.0 
.214 

2.6 
.163 

.16 

D10 

WT300 

vcg 
K 

20.35 
.768 

14.7 
.525 

19.0 
1.19 

.7 

D10 

WT301 

vcg 
K 

21.3 
.804 

11.8 
.421 

21.0 
1.31 

.8 

D10 

WT3R 

vcg 
K 

27.5 

1.038 

19.41 
.693 

17.79 
1.11 

1.0 

D10 

WTG5 

vcg 
K 

19.43 
.733 

16.63  ' 
.596 

13.5 
.844 

.71 

D10 

WToI 

vcg 
K 

20.66 
.78 

20.43 
.73 

15.09 
.943 

.75 

D10 

V/ToII 

vcg 
K 

26.48 
1.0 

21.93 
.783 

19.23 
1.2 

1.0 

D10 

WT6III 

vcg 
K 

41.75 
1.575 

33.9 
1.738 

25.9 
1.62 

1.6 

D10MN     . 

WToIV 

vcg 
K 

24.88 
.939 

20.97 
.749 

13.46 
.841 

.84 

D10 

WTG4 

vcg 
K 

31.56 
1.191 

35.0 
1.25 

22.6 
1.41 

1.3 

D10 

WTG7 

vcg 

K 

35.26 
1.33 

24.45 
.873 

Table 

21.3 
1.33 

4 

1.3 

D10 

V2RTICA1 

CENTERS 

OF  GRAVITY 

Weight 

378* 
WHSC 

210' 

WIviEC 

150' 

WPB 

Coef . 
K 

times 

Crew  & 

Personal 

affects 

vcg 
K 

26.26 
.964 

20.43 
1.048 

15.0 
.938 

.97 

Dl  OMN 

Provisions  vcg 
K 

15.0 
.55 

15.0 
.536 

11.0 
.688 

.55 

D10 

General 
Stores 

vcg 
K 

24.05 
.883 

16.16 
.829 

11.0 

.688 

.85 

D10MN 

Ammo 

vcg 
K 

15.99 
.587 

8.4 
.431 

19.0 
1.19 

.5 

D10MN 

Potable 
Water 

vcg 
K 

7.06 
.26 

8.5 
.30 

4.2 
.263 

.26 

D10 

L.  0. 

vcg 
K 

3.58 
.131 

2.5 

.128 

12.0 
.75 

.13 

D10MN 

Fuel 

vcg 
K 

7.67 

.281 

3.17 
.292 

6.6 
.412 

.3 

D10 

Aircraft 

vcg 
K 

39.01 
1.43 

1.4 

D10 

A/C  Fuel 

vcg 

9.56 
.351 

.351 

D10 

DAVG  28.54  26.91  16.36 

D10  27.25  28.0  16.0 

DlOfruN  27.25  19.5  16.0 


Table     5 
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VCG       SUBROUTINE 
FLOW    CHART 


( 


STACr     VC<5 


REAL     AMD 

LABELED 

STATE  nE-Vrs 


CdfA*  .SSr&fivCr 
Die>-  o>0»\N 


CdCS a  .71  #•  Dto 
COiJ.    z    .7S*0'O 


CCslU  = 

— — — — ^^^— — — - 

C6HZ-L 

.Z^o/o 

C6zm* 

.S^XD'onO 

«r<5-2o>3  = 

.(L*  £>io 

£6300  ■» 

.7  *■  Dio 

CCrtol   - 

.0J*Dto 

C02H.Z. 

D/o 

*J0 


C&S4- 


<£<Mom 


ut<&4 


*D«OMA 


CC\0*or\ 


„  VfS 


£C/»fcso»0 


j£i 


ACm  o*x  ACt*  o^f  wwrtt)* 

A-    ,  OClC&£f\, 

fiffiOCCrCF) 


•*t>iOmt> 


£&<n-<l> 


yes 


CCr4*\/*[0-t 


QjAmmos 


AJO 


Figure    42a 


167 


C  &3o  / »  wrj  i  /  +  CSW  ***m  <X*TttJ 

cases  «:str»t*T«ryc6«a»- 


wno)  /u»r<s  i 


C6W  »    ,SS*-t>to 


C4££  *wnS6  -KU67  •*  wrtf?  )*- 


C<Srf>S-((-  Ctt-T-  +.00.20 -&>0>fr 
C  6  Pf  + .  O  0 13  **«>*  C&CrS  +■ 


C 

c 


~r~ 


e^D 


3 


Figure    42b 


163 


input. 

A  flow  chart  for  this  routine   is  shown  in  Figures 
J*2a  and  **2b. 


3.13.2     Input  List  for  Subroutine 

ACCG LC(EE 

ACWT LC(EE 

AMOCG LC(EE 

AMOWT LC(EE 

CACFUL LC(HH; 

CARGOC LC(EE 

CARGOW LC(EE 

D LC(AAi 

diomn LC(II 

DAVG LC(BB; 

DBMAR LC(VW; 

GR4CG LC(EE 

GR^WT LC(EE 

GR7CG LC(EE 

GR7WT LC(EE 

LEN LC(BB 

LRD LC(Jj; 

NOARM LC(FF 

NOELT LC(FF; 

WACFUL LC(LL 

WCARGO LC(LL 

WFULLB LC(WW! 


VCG 

WT113 LC(VV) 

WT1M..... LC(W) 

WT203 LC(VV) 

WT2M LC(VY) 

WT300 LC(W) 

WT301 LC(VV) 

WT3R LC(VY) 

WT6I ...LC(W) 

WT6II LC(VV) 

WT6III LC(W) 

WT6IV LC(YV) 

WTAMMO LC(LL) 

WTCREW LC(LL) 

WTFUEL LC(LL) 

WTG1 LC(WW) 

WTG2 LC(WW) 

WTG3 LC(WW) 

WTG4 LC(WW) 

WTG5 LC(WW) 

WTG6 LC(WW) 

WTG?.. LC(WW) 

WTLO LC(LL) 


WLSHIP LC(WW) 

WT111 LC(W) 

WT112 LC(W) 
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WTPE LC(LL) 

WTPS LC(LL) 


3.13.3  Statement  Descriptions 

Detailed  statement  descriptions  will  not  be  given  for 
this  routine.  In  generp.l  the  statements  consist  of  vcg 
estimates  using  the  relationships  mentioned  in  the  intro- 
duction or  of  calculations  to  combine  weight  groups  by 
summing  the  vertical  moments  of  individual  groups  and 
dividing  by  the  total  weight. 

Note  that  D10  is  different  from  D10MN  only  if  the 
length  of  raised  deck  is  greater  than  one  half  the  ship 
length.  It  must  also  be  remembered  that  input  moment 
arms  are  given  as  fractions  of  the  main  deck  depth  at 
amidships  and  not  as  a  distance  in  feet. 

3.13.4  Output  List  for  Subroutine  VCG 


CCARGO LC(OO) 

CFULD SCDA 

CFULLD LC(OO) 

CGAC LC(HH) 

CGAMMO LC(HH) 

CGCREW LC(HH) 

CGFUEL LC(HH) 

CGG1 LC(OO) 

CGG2 LC(OO) 


CGG3 LC(OO) 

CGG& LC(OO) 

CGG5 LC(OO) 

CGG6 LC(OO) 

CGG7 LC(OO) 

CGLO LC(HH) 

CGPE LC(HH) 

CGPS LC(HH) 

CLSHIP LC(OO) 


3.13»5  Nomenclature  List 

ACMOM  total  vertical  aircraft  moment,  tons 

AMOMOM  total  vertical  ammunition  moment,  tons 

CACFL  vertical  center  aircraft  fuel,  ft 

CFULD  same  as  CFULLD 

CFULLD  vertical  center  of  ship  at  full  load,  ft 

CG111  vertical  center  of  group  111,  ft 

CG112  vertical  center  of  group  112,  ft 

CG113  vertical  center  of  group  113,  ft 

CG1M  vertical  center  of  group  1M,  ft 

CG203  vertical  center  of  group  203,  ft 

CG2M  vertical  center  of  group  2M,  ft 

CG300  vertical  center  of  group  300,  ft 

CG301  vertical  center  of  group  301,  ft 

CG3R  vertical  center  of  group  3R,  ft 

CG61  vertical  center  of  group  61,  ft 

CG6II  vertical  center  of  group  6ll,  ft 

CG6III  vertical  center  of  group  6III,  ft 

CG6IV  vertical  center  of  group  6lV,  ft 

CGG1  vertical  center  of  group  1 ,  ft 

CGG2  vertical  center  of  group  2,  ft 

CGG3  vertical  center  of  group  3,  ft 

CGG4  vertical  center  of  group  b,   ft 

CGG5  vertical  center  of  group  5»  ft 

CGG6  vertical  center  of  group  6,  ft 

CGG?  vertical  center  of  group  7,  ft 

CGGS  vertical  center  of  general  stores,  ft 
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CGOMOM 

CGPV 

CGPW 

CLSHIP 

CPO 

D 

DO 

DIO 

DIOMN 

D20 

ENL 

G4K0M 

G7MOM 

HLEN 

I 

OFF 

WT1M 

WT111 

WT112 

WT113 

WT203 

WT2M 

WT300 

WT301 

WT3R 

WT6I 

WT6II 


171 

vertical  moment  of  cargo,  tons 

vertical  center  of  provisions,  ft 

vertical  center  of  potable  water,  ft 

vertical  center  of  light  ship,  ft 

same  as  NCPO  in  MAIN 

same  as  ENDDAY  in  MAIN 

depth  at  F.P. ,  ft 

depth  amidships  to  main  or  raised  deck,  ft 

depth  amidships  to  main  deck,  ft 

depth  at  A. P. ,  ft 

same  as  NENL  in  MAIN 

total  group  k   moment,  tons 

total  group  7  moment,  tons 

half  the  length  of  the  ship,  ft 

index 

same  as  NOFF  in  MAIN 

weight  group  1M,  tons 

weight  group  111,  tons 

weight  group  112,  tons 

weight  group  113,  tons 

weight  group  203,  tons 

weight  group  Z1A,   tons 

weight  group  300,  tons 

weight  group  301,  tons 

weight  group  3K»  tons 

weight  group  6l,  tons 

weight  group  6ll,  tons 
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WT6III        weight  group  6lII,  tons 

WT6IY         weight  group  6lV,  tons 
XLDMOM         total  load  moment,  ton-ft 
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3.14  Subroutine  COST 

3.14,1  Introduction 

This  subroutine  is  based  exclusively  on  the  pre- 
liminary cost  estimating  procedure  currently  used  by  the 
Coast  Guard.   To  use  the  procedure,  the  weight  of  each 
of  the  seven  light  ship  weight  groups  must  be  known. 
Linear  graphs  for  man  hours  and  material  cost  are  then 
entered  with  these  values  to  determine  the  cost  and  labor 
required  for  each  weight  group. 

Various  relationships  are  also  given  for  determining 
design  costs,  construction  services  cost  and  miscellaneous 
costs.  Miscellaneous  costs  include  initial  outfit  equip- 
ment such  as  life  jackets  as  well  as  retrofit  costs,  spare 
parts  and  administrative  or  resident  inspector's  office 
costs. 

The  curves  used  in  the  procedure  were  converted 
to  formulas  and  the  procedure  was  programmed  as  subroutine 
COST.  A  flow  chart  for  this  routine  is  given  in  Figure 

43. 

All  of  the  material  cost  curves  are  based  on  1959 
prices.   Indices  are  used  to  correct  the  prices  to  the 
year  the  ship  will  be  built.   These  indices  are  based 
on  Bureau  of  Labor  Statistics  data  for  the  particular 
weight  group. 

The  labor  rate  for  all  of  the  weight  groups  is  assumed 
to  be  the  same.   The  labor  rate  is  inputed  as  an  average 
dollars  per  hour  for  the  year  the  ship  will  be  built. 
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Material  cost  indices  and  labor  rate  data  are  shown 

in  Figures  44  and  45.  There  is  a  large  fluctuation  in 
costs  between  shipyards  so  these  numbers  can  only  be 
viewed  as  average  values. 

In  addition  to  the  corrections  for  inflation,  eighty 
percent  is  added  to  labor  cost  for  overhead  and  ten 
percent  is  added  to  all  costs  for  profit. 

Initial  outfit  costs  are  assumed  to  be  two  percent 
of  the  total  cost.  Spare  parts  cost  is  the  sum  of  the 
cost  of  spares  for  propulsion  machinery,  electrical 
machinery,  electronics  and  outfit.   The  cost  of  spares 
as  a  percentage  of  the  weight  group  material  cost  is 
given  be low i 

Group  2 9% 

Group  3 8# 

Group  4  25% 

Group  5  10# 

Retrofit  cost  is  assumed  to  be  4  percent  of  the 
total  cost.  Administrative  socts  are  assumed  to  be  3.5 
percent  of  the  total. 


3.14.2  Input  List  for  Subroutine  COST 


CSTG2 LC(MM) 

CSTG4 LG(MM) 

CSTG7 LC(MM) 

DBMAR LC(WW) 

DOLHR LC(NN) 


G1IND LC(NN) 

G2IND LC(NN) 

G3IND LC(NN) 

G5IND LC(NN) 

G6IND LC(NN) 


176 


— ; .  -  -:  &ws£*o.  -frfm $<?*:■.  si^f-isit^^i-- 


mmmmm 


pat*    F&o/n 


WHOLESALE  -rjPJLKjrrfHLZ* 


YEAR 
Figure   44 


177 


TT 


:?f? 


MM 


lL 


4?' 


::!::- 


l::ru~f 


.  ,i... ».  i— 


rr|.ti: 


£ 


-!_•::):  :r. 


1 


^•g|5?pg| 


TS 


=  E?^ 


■;i- 


gfeti 


t  '   r 


r  -r 


i;.  :, 

a. 


:t:: 


— X,.  -  -     .„ 


::tz- 


^4^ 


gja 


o 

s 


o 

3 


o 


3 


«H  /$ 


31VH     a08VT 


178 
JOPT . .LC(CC )      WTG3 • LC (WW) 

LEN SCDA      WTG4 LC(WW) 

MTYPE LC(CC)      WTG5 LC(WW) 

WTG1 .LC  (WW)      WTG6 LC  (WW) 

WTG2 LC(WW)      WTG7 LC(WW) 

♦only  if  JOPT  -  2 

3.14.3  Statement  Descriptions 
W1«WTG1*DBMAR 

W7=WTG?*DBMAR 

The  design  and  builders  margin  is  evenly  divided 
among  the  weight  groups  for  costing  purposes. 

G1MH=10.*»(AL0G10(W1)*. 75217+2. 97276) 
• 

G7MH=10.»*(AL0G10(W7)*.75217+2. 97276) 

These  formulas  represent  the  relationships  between 

man-hours  and  the  weight  of  individual  weight  groups. 

These  relationships  are  assumed  linear  when  plotted  on 

log-log  paper. 

G1MC=10.**(AL0G10(W1)*. 954243+2. 790484) 
G2MC=CSTG2 

IF(J0PT.EQ.2)    GOTO  1 

G2MC=10.**(AL0G10(W2)*1. 018885+3. 63^328) 
IF(MTYPE.EQ.2)  G2MC-G2MC+1 500000. 
1  G3MC«10,**(AL0G10(W3)*1. 072551+3.667812) 
G4MC=CSTG4 
G5MC»10.**(ALOG10(W5)*. 870906+3. 974192) 

g6mc=io.**(alogio(w6 )*i. 068265+3.603833) 

G7MC-CSTG7 

These  formulas  represent  the  material  cost  curves 
used  by  the  Coast  Guard.     Some   of  these  costs,   CSTG2, 
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CSTG4  and  CSTG7,  are  input.   Group  2  cost  is  an  input 

only  if  the  machinery  plant  is  specified.   Otherwise, 
the  machinery  cost  is  calculated.   The  formula  used 
represents  the  cost  of  diesel  plants.  An  extra  $1.5 
million  is  added  if  a  CODOG  plant  is  to  be  installed. 

G1LC=G1MH*D0LHR*1 . 98 
. 

G7LC»G7MH*DOLHR»l . 98 

The  labor  cost  by  weight  group  is  calculated  here. 

The  value,  1.98,  represents  an  overhead  factor  of  1.8 

times  a  profit  factor  of  1.1. 

G1MC»G1MC*G1IND»1 . 1 
IF(J0PT.EQ.2)  GOTO  2 
G2MC«G2MC*G2IND 
2  G3MC=G3r.iC*G3IND«l .  1 
G5MC-G5KC *G5IND»1 . 1 
G6MC»G6MC*G6lND*l . 1 
G2MC«G2MCn.l 

g4mc»g4i.:c*i.i 

G7MC«G7MC*1.1 

In  these  formulas  the  material  cost  is  corrected 

for  inflation  and  a  factor  of  1.1  is  included  for  profit. 

Note   that  the  costs  given  as  input  are  not  corrected 

for  inflation.      Therefore,   the  cost  of  the  equipment 

must  be   corrected  for  inflation  before  being  inputed. 

LABCST=G1LC+G2LC+G3LC+G4LC+G5LC+G6LC+G7LC 
MATCST=G1MC+G2MC+G3MC+G4KC+G5MC+G6MC+G7MC 

The   total  labor  and  material  costs  are  calculated 

in  these   statements. 
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L0A»1.08»LEN 

FACT0R=(.11^*L0A-1,6)/100, 
DLC«FACTOR*LABCST 
FACTOR=(-.01*LOA+7.5)/lOO. 

DMC=FACTOR*MATCST 

These  statements  calculate   the   design  labor,   DLC, 

and  design  materials,   DMC,   costs  for  the  first  ship  of 

a  class  from  the  beginning  of  the  conceptual  design  stage. 

The   two  values  of  FACTOR  are  based  on  Coast  Guard  curves* 

Since   the  length  overall  is  required  in  these  curves,   its 

value   is  estimated  as  eight  percent  greater  than  the 

length  between  perpendiculars, 

FACT0R«(.0Hf»L0A+7.)/l00. 
CLC=FACT0R*LABCST 
FACT0R=(-.015»L0A+6.5)/l00. 
CMC»FACTGR*MATCST 

The  construction  services  labor  and  material  costs 

are  calculated  here.     Again  the  formulas  for  the   two  values 

of  FACTOR  are  based  on  Coast  Guard  curves, 

TOTCSTxL/iBCST+MATCST+DLC+DMC+CLC+CMC 

0RICST=.02*T0TCST 

SPARES«(.09*G2MC+.08*G3MC+.25*G4MC+.1*G5MC)/1.1 

RETR0=.04*T0TCST 
ADMIN=.035*T0TCST 

The  four  items   that  are   included  in  the  miscellaneous 

cost  are  calculated  here.     Since   profit  is  included  in 

the  individual  weight  group  material  costs  it  is  divided 

out  when  calculating  the  cost  of  spares, 

CSTG1-G1LC+G1MC 

CSTG7-G7LC+G7MC 

Total  weight  group  costs  are  calculated  here 


181 

DCST»DLC+DMC 
CCST-CLC+CMC 

MCST-ORICST+SPARES+RETRO+ ADMIN 

TOTCST-TOTCSI+MCST 

RETURN 

END 

Finally,  the  total  design  cost,  construction  services 

cost,  miscellaneous  cost  and  the  total  first  ship  cost  are 

calculated, 

3.1^.4  Output  List  for  Subroutine  COST 


CSTG6 LC(MM) 

CSTG7 LC(MM) 

DCST LC(MM) 

MOST LC(MM) 

TOTCST LC(NN) 


CCST LC(MM) 

CSTG1 LC(MM) 

CSTG2? LC(MM) 

CSTG3 LC(MM) 

CSTG4 LC(MM) 

CSTG5 LC(MM) 

*  only  if  JOPT  ■  1 

3.1^.5  Nomenclature  List 

The  definition  of  all  variables  is  the  same  as  that 
given  in  the  MAIN  program  nomenclature  except  for  the 
following  1 

ADMIN         administrative  costs,  dollars 
CLC  construction  services  labor  cost,  dollars 

CMC  construction  services  material  cost,  dollars 

DLC  design  labor  cost,  dollars 

DMC  design  materials  cost,  dollars 

FACTOR         see  text 
G1LC  weight  group  1  labor  cost,  dollars 
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G2LC 

G3LC 

G4LC 

G5LC 

G6LC 

G7LC 

G1MC 

G2MC 

G3MC 

G4MC 

G5MC 

G6MC 

G7MC 

G1MH 

G2MH 

G3MH 

(&MH 

G5MH 

G6MH 

G?MH 

LABCST 

LOA 

MATCST 

ORICST 

RE 

RETRO 

SFCH 


weight  group  2  labor  cost,   dollars 

weight  group  3  labor  cost,   dollars 

weight  group  4  labor  cost,   dollars 

weight  group  5  labor  cost,   dollars 

weight  group  6  labor  cost,   dollars 

weight  group  7  labor  cost,   dollars 

weight  group  1  material  cost,   dollars 

weight  group  2  material  cost,   dollars 

weight  group  3  material  cost,  dollars 

weight  group  4  material  cost,   dollars 

weight  group  5  material  cost,   dollars 

weight  group  6  material  cost,   dollars 

weight  group  7  material  cost,   dollars 

weight  group  1  man-hours 

weight  group  2  man-hours 

weight  group  3  man-hours 

weight  group  4  man-hours 

weight  group  5  man-hours 

weight  group  6  man-hours 

weight  group  7  man-hours 

total  construction  labor  cost,   dollars 

length  overall,   ft 

total  construction  material  cost,   dollars 

initial   outfit  cost,   dollars 

same   as  RGEND  in  MAIN 

retrofit  cost,   dollars 

same   as  SFCHHP  in  MAIN 


SFCM 

SPARES 

YE 

Wl 

W2 

W3 

W6 
W7 


same  as  SFCMHP  in  MAIN 
spare  parts  cost,  dollars 
same  as  VEND  in  MAIN 

weight  group  1  with  margin,  tons 

weight  group  2  with  margin,  tons 

weight  group  3  with  margin,  tons 

weight  group  k   with  margin,  tons 

weight  group  5  with  margin,  tons 

weight  group  6  with  margin,  tons 

weight  group  7  with  margin,  tons 
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3.15  Subroutine  SEASPD 

3.15.I  Introduction 

This  routine  is  based  on  the  work  of  Michael  Shen, 
Peter  Gale  and  Kenneth  Walker  as  reported  in  References 
4  and  5,  In  these  references  is  given  a  computer  program 
for  computing  the  average  sea  speed.   This  computer 
program  was  altered  by  dropping  the  option  to  use  a  bow 
sonar  dome  and  by  deleting  the  read  and  write  statements. 
The  program  was  also  rewritten  in  ASA  standard  FORTRAN. 
Except  for  these  changes,  subroutine  SEASPD  is  the  computer 
program  developed  by  these  three  men.   The  interested 
reader  is  referred  to  these  references  for  the  development 
of  the  routine. 

Briefly,  the  routine  assumes  1  (1)  that  the  ship 
is  operating  in  the  North  Atlantic  Ocean  over  a  long 
period  of  time  1  (2)  that  it  is  steaming  into  head  seas 
one  half  of  the  time t    (3)  that  the  remainder  of  the  time 
(when  steaming  at  other  headings)  it  is  not  forced  to 
reduce  speed  due  to  ship  motions.   Part  of  the  time  the 
ship  is  steaming  into  head  seas,  its  speed  will  be  limited 
be  the  sea  state.   In  the  routine  this  is  simplified  to  . 
a  single  sea  state  limited  speed.   The  remainder  of  the 
time  the  ship  is  assumed  to  be  able  to  make  its  maximum 
sustained  speed.   The  percentage  of  the  time  the  ship 
is  limited  in  speed  is  estimated  based  on  average  wave 
heights  in  the  North  Atlantic  Ocean  and  on  ship  character- 
istics. 
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The  average  speed  at  sea  is  calculated  assuming 

that  the  vessel  will  travel  at  its  maximum  sustained 

speed  whenever  it  is  not  limited  in  speed  by  the  sea 

state*  When  the  vessel  is  steaming  at  headings  other 

than  into  head  seas,  its  speed  is  not  limited.   This  is 

assumed  to  be  one  half  the  time*  When  steaming  into 

head  seas  its  speed  will  be  limited  only  part  of  the 

time*  Both  the  limiting  speed  and  the  percentage  of 

time  that  it  applies  are  calculated  in  SEASPD* 

This  routine  has  been  included  to  indicate  the 
effect  that"  changes  in  the  vessels  characteristics  will 
have  on  its  seakeeping  ability*   The  many  assumptions 
made  in  deriving  the  average  sea  speed  make  the  speed 
predicted  practically  meaningless*   However,  the  speed 
calculated  can  be  used  to  compare  two  designs  to  decide 
which  will  give  the  better  sea  keeping  performance* 

A  flow  chart  for  this  routine  is  shown  in  Figure  Jf6* 

3*15.2  Input  List  for  Subroutine  SEASPD 

B LC  (GG)        DPTRY LC(OG) 

CP LC  (GG)        H LC(GG) 

CWP LC  (GG)        LEN LC(BB) 

CX LC(GG)        VSUS LC(CC) 

5* 15. 3  Statement  Descriptions 

For  a  detailed  discussion  of  this  routine  see  Refer- 
ences k   and  5* 


SEASPD       SUBROUTINE 
FLOW       CHART 
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3/4*5  •^u$^*PRi?^or  Subroutine  SEASPD 

i  ic.  iTheAfonly.  output  is  the  value  of  average  sea  speed, 

AVSP. . 
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3.15«5     Nomenclature  List 

The  definitions  used  are   the  same  as  those  given  in 
References  b  and  5  except  for  the   following! 


Value  used 
in  SEASPD 

B 

LEN 

CWP" 

ll   *" 

DPTRt 

KY 

IYY 

f  ^rrr.u  *     "  I  i  loi: 

JY 

tr   ■■•or«vi*r  t    ;>'. 
IY 

vsus ' 
vlr" 


Value  used  in 
references 

BEAM 

SL 

CW 

B 

DISPO 

DKY 

YYI 

YJ 

YI 

wv 

VPX 


3 .  i  t .  <       If:; 


rou  tir.e . 
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3.16     Subroutine   OUTPUT 

3.16.1  Introduction 

All  program  output  is  printed  by  subroutine  OUTPUT 
except  for  error  messages.  A  simplified  flow  chart  for 
this  routine  is  shown  in  Figure  47.  Subroutine  OUTPUT 
is  divided  into  two  main  parts.  The  first  prints  the 
input  data  and  the  second  prints  the  output  generated. 

If  an  error  ties  sage  is  printed,  only  part  of  the 
output  data  is  available  and  the  data  that  is  available 
is  probably  not  correct.   Therefore,  when  an  error  message 
occurs,  only  the  input  data  is  printed.  Variable  II  has 
a  value  of  1  in  this  case.  If  the  program  successfully 
completes  all  proceeding  steps,  II  is  assigned  a  value 
of  2  and  output  data  is  also  printed. 

Subroutine  OUTPUT  consists  primarily  of  write  and 
format  statements.  Some  calculations  are  made  but  only 
to  convert  previously  calculated  values  to  a  desired 
output  form  such  as  determining  ratios  between  weight 
groups  and  full  load  weight. 

Only  two  decision  points  are  used.   One  is  the  decision 
on  whether  or  not  to  print  output.   The  other  checks  the  >*•; 
option  used  for  inputting  vehicle  performance  data  and 
chooses  the  correct  set  of  write  and  format  statements 
to  print  the  input  data. 

3.16.2  Input  for  Subroutine  OUTPUT 

All  values  which  are  printed  are  inputs  to  this 
routine. 


OUTPUT   SUBROUTINE 
FLOW    CHART 
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Figure    47 
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3#l6#3  Statement  Inscriptions 

Detailed  descriptions  of  the  statements  of  this  routine 
will  not  be  given  since  most  of  the  statements  are  self- 
explanatory.   Use  of  the  sample  outputs  given  in  Appendix 
B  will  be  of  use  in  understanding  the  format  statements 
used. 

3.16.4  Output  for  Subroutine  OUTPUT 

No  values  are  generated  in  this  routine  for  use 
elsewhere. 

3.16.5  Nomenclature 

All  the  variables  used  in  this  routine  have  been 
defined  previously  except  for  some  ducsay  variables  whose 
definition  should  be  clear  from  context. 


CHAPTER  IV 
EVALUATION  OP  THE  PROGRAM 

A  number  of  test  runs  were  made  on  an  IBM  370 ,  model 
168,  computer  located  in  the  Information  Processing  Center 
at  M.I.T.  The  test  runs  were  made  using  input  data  for 
four  existing  designs,  the  3?8*  WHEC,  the  210*  WMEC,  the 
150'  WPB,  and  the  Navy's  Patrol  Frigate  design.   The  Cutter 
Model  produces  feasible  designs  for  each  of  these  ships 
except  the  150'  WPB. 

The  design  for  the  150*  WPB  failed  because  a  raised 
deck  length  greater  than  the  ship  length  was  required. 
The  actual  ship  does  not  have  a  raised  deck.   The  probable 
cause  of  this  failure  lies  with  the  available  area  and 
volume  relations  used  in  the  program.   These  relations 
were  developed  for  ships  with  lengths  greater  than  400 
feet  and  probably  have  significant  errors  when  extended 
to  ships  as  short  as  150  feet.  Recommendations  for  cor- 
recting this  error  are  given  in  Chapter  V. 

A  comparison  between  the  results  produced  by  the 
Cutter  Model  and  the  actual  ships  are  shown  in  Tables 
6a  and  6b  for  the  378*  WHEC  and  in  Tables  7a  and  7b  for 
the  210*  WMEC.   The  differences  are  reasonably  small  in 
most  areas)  however,  a  few  do  need  discussion. 

The  largest  difference  with  both  ships  occurs  in 
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Table  6a 

Comparison  of  Cutter  Model  Output  to  the  Actual 
Ship  for  a  378*  WHEC 


& 


Beam 

Draft 

DO 

D10 

D20 

Davg 

Disp. 

ft 

ft 

ft 

ft 

ft 

ft 

tons 

Ship  350.  42.0   14.3  3^.75  27.25  27.5   28,54  2968.0 
Model  350.  43.32  14.6  38.33  27.83  28.18  29.64  3061.5 


WTG1 
WTG2 
WTG3 
WTG4 

WTG5 
WTG6 

WTG7 


WEIGHTS 
Ship     Model 


VCG 


917.46 
383.6 
88.12 
80.2 
268.8 

325.1 
39.0 


962.8 
383.6 
164.8 

80.2 
262.3 
322.3 

39.0 


Ship 


Model 


19.12 

19.42 

11.45 

11.52 

23.49 

21.76 

31.56 

32.52 

19.43 

19.76 

25.63 

25.88 

35.26 

36.30 

WTLSHIP 


2102.2  2215.0 


19.72  19.98 


FUEL  661.8  653.2 

LUB  OIL  13.7  16.0 

PER  EFFECTS  ->  22»°3  7.8 

PER  STORES  107.4  95.1 

CARGO                     0.  0. 

AUMO  37.5  37.5 

AIRCRAFT               3.4  3.4 

A/C   FUEL  19.6  19.6 


7.67 
3.58 

1 26.26 
11.35 

15.99 

39.01 

9.56 


8.35 

3.62 

27.00 

27.00 

14.36 

16.42 

40.08 

9.74 


WTFLLD 


2967.6  3061.4 


16.61 


17.20 


Table  6b 
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Comparison  of  Cutter  Model  Output  to  the  Actual 
Ship  for  a  378*  WH2C 


Hull  Volume 
Deckhouse  Volume 

Total  Arrangements  Area 
Hull  Arrangements  Area 
Deckhouse  Arrangements  Area 

Office  Space 
Messing  Facilities 
Crew  Special 
Officer  S.R. 
Officer  Sanitary 
CPO  S.R. 
CPO  Sanitary- 
Crew  Berthing 
Crew  Sanitary 
CO  S.R.,  Cabin  &  Pantry 
Commissary  Stores 
Other  Stores 
WorkBhops 
Repair  Lockers 
Chain  Lockers 
Uptakes 

Steering  Gear  &  Windlass 
A/C  &  Fan  Spaces 
I.C.  Spaces 
Aux.  Machinery  Spaces 
Pilot  house,  Chartroom  &  CIC 
Passages 
Input  Area 
Excess  Area 


Ship 

277600 
95611 

30863 
19615 
11248 

^93 
2713 

725 
1157 

243 
1040 

250 
3641 

716 

755 

600 

3015 

1161 

148 

47 
462 

736 

369 

312 
i860 
1125 
2567 
5283«» 
1445 


Model 

27373** 
94705 

29664 
18522 
11142 

510 
2713 

787 
1168 

192 
1040 

208 
3168 

720 

753 

652 

3231 

1254 

100 

779 

396 

312 

2180 

1111 

2450 

5900 

0 


cu  ft 
cu  ft 

sq  ft 
sq  ft 
sq  ft 

sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 


*  460  sq  ft  included  in  input  area 

»*  Excess  area  was  used  for  crew  rec.  rooms 
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Table  7a 

Comparison  of  Cutter  Model  Output  to  the  Actual 
Ship  for  a  210'  WMEC 


£t>P 
ft 


Beam     Draft  DO         D10  D20  Davg 

ft  ft  ft         ft  ft  ft 


Ship     200,      34.00     10.5     20.75     19.50     20.90     26.91 
Model  200.     35.02       8.98  21.50     19.50     19.70     24.97 


Disp. 
tons 

990.0 
959.1 


WEIGHTS 

VCG 

Ship 

Model 

Ship 

Model 

WTG1 

379.5 

347.8 

16.4 

13.92 

WTG2 

102.6 

102.6 

9.5 

9.55 

WTG3 

43.7 

75.4 

17.5 

22.28 

WTG4 

26.1 

26.1 

34.2 

33.31 

WTG5 

110.4 

109.6 

16.6 

19.88 

WTG6 

104.4 

106.8 

22.6 

25.93 

WTG7 

8.5 

8.5 

23.4 
17.1 

23.21 

WTLSHIP 

775.2 

776.8 

17.40 

PUEL 

152.3 

137.3 

8.17 

8.40 

LUB   OIL 

3.8 

3.9 

2.5 

2.54 

CREW 

}8.6 

5.5 

20.43 

18.92 

PER  EFFECTS 

3.2 

20.43 

18.92 

PER  STORES 

47.3 

30.4 

10.99 

11.95 

CARGO 

0. 

0. 

AMMO 

2.6 

2.6 

8.4 

8.00 

AIRCRAFT 

0. 

0. 

A/C  FUEL 

0. 

0. 

WTFLLD 


990.0 


959.5 


15.38  15.87 


Raised  Deck  Length 


Ship 
164. 


Model 
120. 


Table  7b 
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Comparison  of  Cutter  Model  Outout  to  the  Actual 
Ship  for  a  210'  WMEC 


Hull  Volume 
Deckhouse  Volume 

Total  Arrangements  Are* 
Hull  Arrangements  Area 
Deckhouse  Arrangements  Area 

Office  Space 

Messing  Facilities 

Crew  Special 

Officer  S.R. 

Officer  Sanitary 

CPO  S.R. 

CPO  Sanitary 

Crew  Berthing 

Crew  Sanitary 

CO  S.R.,   Cabin  &  Pantry 

Commissary  Stores 

Other  Stores 

Workshops 

Repair  Lockers 

Chain  Lockers 

Uptakes 

Steering  Gear  &  Windlass 

A/C  &  Fan  Spaces 

I.C.  Spaces 

Aux.  Machinery  Spaces 

Pilot  house ,  Chartroom  &  CIC 

Passages 

Input  Area 

Excess  Area 


Ship 

124700 
16100 

11242 
9345 
1897 


Model 

123049 
17840 

11021 
8922 
2099 


130 

192 

1452 

1381 

324 

313 

567 

657 

193 

108 

283 

320 

75 

64 

1658 

1254 

384 
406 

285 

398 

480 

393 

1294 

1233 

544 

382 

110 

100 

42 

40 

0 

0 

270 

390 

165 

151 

180 

180 

820 

832 

422 

486 

739 

910 

704 

704 

0 

247 

cu  ft 
cu  ft 

sq  ft 
sq  ft 
sq  ft 

sq  ft 
eq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
eq  ft 
sq  ft 
sq  ft 

sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
sq  ft 
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weight  group  3,  electric  plant.   This  difference  results 

from  the  fact  that  greater  generating  capacity  relative 
to  the  average  load  has  been  built  into  the  model  than  is 
installed  in  existing  ships.   This  increase  in  generating 
capacity  reflects  current  Coast  Guard  policy  for  new 
designs  and  is  not  an  error  in  the  program. 

The  differences  in  weight  group  1  reflect  the  fact 
that  the  hull  dimensions  are  not  exactly  the  same  for 
the  program  designs  and  the  actual  ships.  Weight  group 
1  is  a  function  primarily  of  the  hull  cubic  number. 

The  fuel  weight  is  calculated  in  the  program  based 
on  specific  fuel  consumption!  endurance  range  and  endurance 
speed.  An  attempt  has  been  made  to  approximate  the  actual 
fuel  weight  by  varying  one  or  more  of  these  quantities 
but  an  exact  match  is  difficult  to  obtain. 

The  only  other  significant  weight  difference  occurs 
in  personnel  stores  which  includes  potable  water.  The 
model  results  are  less  than  the  weights  for  the  actual 
ships  due  to  the  fact  that  the  amount  of  potable  water 
per  man  has  been  reduced  from  the  amount  used  on  past 
designs.  This  is  in  accordance  with  current  Coast  Guard 
policy  for  future  designs. 

The  vertical  centers  of  gravity  show  quite  a  bit 
of  variation  when  compared  to  the  actual  ships.  However > 
the  values  in  the  program  are  only  average  values  and 
will  in  general  not  agree  with  the  values  for  any  one 
ship. 


197 
In  the  case  of  the  210'  WMEC,  a  raised  deck  is 

required.   The  model  predicted  a  raised  deck  length  44 

feet  shorter  than  that  of  the  actual  ship.  This  is  another 

indication  that  the  available  area  and  volume  relations 

are  inaccurate  for  smaller  ships. 

Hull  areas  are  given  in  Tables  6b  and  7b.  The  areas 
used  in  the  cutter  model  for  berthing  and  sanitary  spaces 
are  based  on  current  minimum  Coast  Guard  standards.   These 
areas  in  general  are  lower  than  the  areas  on  current 
ships. 

Several  runs  were  also  made  using  data  for  the  Navy's 
Patrol  Frigate  design.  At  the  same  time,  runs  were  rcado  on 
the  Navy's  destroyer  model,  DD07,  using  similar  data. 
The  results  of  these  runs  are  shown  in  Tables  8  and  9* 

The  Patrol  Frigate  design  is  of  interest  because 
it  is  at  the  lower  limit  of  applicability  for  the  des- 
troyer model  and  at  the  upper  limit  of  applicability  for 
the  cutter  model.  Agreement  between  the  values  calculated 
by  the  two  models  and  the  actual  ship  as  shown  in  Table 
8,  is  reasonably  good.  Again,  the  group  3  weight  is  highest 
for  the  cutter  model  because  3  -  2000  KW  generators  are 
used  instead  of  the  3  -  1000  KW  generators  on  the  actual 
ship.  Also,  it  appears  that  the  greater  emphasis  placed 
on  habitability  by  the  Coast  Guard  is  reflected  in  a 
greater  group  6  weight. 

Table  9  gives  the  displacements  predicted  for  greater 
lengths  than  the  actual  ship  length.  As  the  ship  length 
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Comparison  of  Cutter  Model  Output  to  the  Actual 
Ship  and  to  DD07  Estimate  for  the  Patrol  Frigate  Design 


Lbr>   Beam  Draft 
ft* 


ft 


ft 


DO 
ft 


Ship  408.  45.0   14,5 

Model  405.  45.95  15.34  41.07 

DDO?  405.  46.72  15.3  41.00 


D10 
ft 


28.92 
33.99 


D20 
ft 


29.33 
34.80 


Davg   Disp. 
ft    tons 


31.01 
35.30 


3517. 
3706.1 

3766.3 


WTG1 
WTG2 
WTG3 
WTG4 

WTG5 
WTG6 
WTG7 


WEIGHTS 
Ship     Model 


1222.36 
256.1 
155.36 
92.58 

351.3 
268.82 
95.68 


1199.6 
258.0 
211.2 
70.7 
317.5 
333.8 

Q5*7 


DD07 

1331.2 
240.0 
156.6 
70.7 
381.3 
238.0 

85.7 


WTLSHIP 

2442.20 

2637.4 

2666.3 

FUEL 

830.0 

818.7 

845.3 

LUB   OIL 

20.22 

20.2 

0.6 

CREW  &  FE 

20.96 

20.9 

20.9 

PER  STORES 

68.05 

66.6 

84.0 

CARGO 

67.O 

75.0 

80.9 

AMMO 

58.8 

58.8 

58.8 

AIRCRAFT 

10.0 

10.0 

10.0 

A/C   FUEL 

0. 

0. 

0. 

WTFLLD 


3517. 


3707.7    3766.9 


Table  9 
PATROL  FRIGATE 
Pull  Load  Displacement  vs.  Length  Estimates 
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HP 

Displacement 

Cutter 
DD07             Model 

405. 

3766.3 

3707.7 

415. 

3652.2 

3679.5 

425. 

3537.1 

3672.O 

435. 

.    3401.6 

3696.0 

445. 

3458.8 

3764.1 

455. 

3529.2 

3874.7 

465. 

3640.5 

3985.9 
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is  increased ,  the  displacements  predicted  by  the  cuttor 

model  begin  to  diverge  from  those  predicted  by  DD07. 

This  indicates  that  the  cutter  model  begins  to  lose 

accuracy  for  ship  lengths  above  400  feet,  the  assumed 

upper  limit  of  the  cutter  model.  Table  9  also  shows  that 

the  two  models  predict  different  optimums.  This  is  a 

result  of  the  assumptions  built  into  the  two  models. 
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CHAPTER  V 

RECOMMENDATIONS  AND  CONCLUSIONS 

For  the  378 ■  WHEC  and  210*  VfMEC,  the  full  load 
displacements  predicted  by  the  cutter  model  are  within 
four  percent  of  the  actual  full  load  displacements  of 
these  two  vessels •   This  is  better  accuracy  than  is 
obtainable  with  other  currently  available  models  for  the 
size  range  of  Coast  Guard  vessels.   However,  only  a  limited 
number  of  designs  were  used  to  develop  the  cutter  model 
and  these  same  designs  were  used  to  check  it.   The  program 
should  be  used  with  caution  until  other  designs  are  avail- 
able to  confirm  the  relationships  used  in  the  model. 

Although  the  cutter  model  appears  to  give  accurate 
results  in  its  current  form,  there  are  several  areas 
where  significant  improvements  could  be  made.   The  two 
most  important  improvements  are t  (1)  development  of  avail- 
able, arrangement  area  and  volume  relationships  which  are 
more  applicable  to  ships  of  the  size  range  used  in  the 
model  .and  (2)  expansion  of  the  Cr  array  used  in  calculating 
horsepower,  requirements.  The  electrical  generator  sizing 
criteria  should  also  be  given  careful  consideration  to 
determine  if : » the  sizes  which  result  are  what  is  really 
desired*  a-e...  ;■>■•:  •  . 
.;.  ;.  The-,  development  of  area  and  volume  relationships 
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will  require  considerable  developmental  work.  As  a 

minimum,  the  variation  of  area  with  hull  depth,  beam, 

Cp,  Cx,  tankage  volume,  and  machinery  box  length  must 

be  considered.  A  computer  program  should  be  developed 

to  calculate  tankage  volume  and  arrangements  area  for 

systematic  changes  in  these  variables.   The  results  could 

then  be  converted  into  formulas  for  use  in  the  cutter 

model. 

The  expansion  of  the  Cr  array  could  probably  best 
be  accomplished  by  fairing  data  from  another  model  test 
series,  such  as  series  64,  into  the  data  from  Taylor's 
Standard  Series.   Where  necessary  some  data. may  have  to  be 
assumed  that  is  not  available.  An  expansion  of  the 
array  to  speed  length  ratios  of  3*0  and  Cy  values  up 
to  0.008  would  be  very  desirable.  With  the  current 
array  limits,  many  ships  can  not  be  run  which  may  be 
very  desirable.  Also,  if  the  sustained  horsepower  is 
given  as  an  input,  there  are  many  times  when  the  array 
limits  are  exceeded  during  an  iteration  even  though  the 
final  speed  value  may  be  well  within  the  limits  of  the 
array. 

Finally,  the  generator  capacity  predicted  by  the 
model  appears  to  be  excessive.  This  criteria  should 
be  given  careful  consideration.  Existing  378'  VfHBC 
cutters  seldom  exceed  the  capacity  of  one  500  KW  generator 
in  normal  steaming i  however,  the  criteria  built  into  the 
cutter  model  would  result  in  1500  Kff  generator  sets  being 
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installed. 

Some  general  recommendations  on  the  use  of  the  model 
are  also  in  order.  First,  a  ship  synthesis  model  is  only 
a  design  tool.  Complete  reliance  should  never  be  placed 
on  a  modelt  but,  if  checked  by  hand  calculations,  the 
model  will  greatly  inhance  the  engineer's  ship  design 
capability.  Also,  as  a  tool,  the  model  must  be  kept 
current.   This  implies  that  periodic  reviews  must  be  made 
of  the  empirical  relations  used  to  insure  that  they  are 
in  agreement  with  current  policy.  Every  new  design  built 
should  also  be  plotted  on  the  curves  included  in  Chapter 
III  as  additional  data.  Changes  should  be  made  if  the 
new  data  is  not  in  agreement  with  the  old. 

It  is  expected  that  the  model  will  most  often  be 
used  for  designing  new  ships.  Studied  most  often  will 
probably  be  changes  in  length  and  Cp.   These  values  have 
been  arranged  so  that  this  can  be  done  by  changing  only 
one  input  card.   However,  changes  in  any  of  the  input 
variables  can  be  investigated. 

If  the  model  is  used  to  investigate  the  sensitivity 
to  changes  in  the  design  assumptions,  it  will  be  necessary 
to  change  statements  in  the  program  itself.   If  many  of 
these  studies  are  being  performed,  some  of  the  fixed 
values  in  the  program  could  be  read  in  on  input  cards. 

One  constraint  on  all  ship  synthesis  models  is  that 
variations  in  hull  form  can  not  be  investigated.   If 
hull  forms  much  different  from  those  of  current  designs  are 
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desired,  a  new  model  has  to  be  developed.  It  must  be 
remembered  that  a  ship  synthesis  model  is  really  a  multi- 
dimensional linearization  of  a  complex  problem.  If 
large  departures  are  made  from  the  data  base  designs* 
the  errors  will  increase  in  magnitude. 

Finally,  the  cutter  model  could  be  programmed  as 
part  of  a  ship  optimization  routine,  but  it  is  recommended 
that  this  not  be  done.  Such  a  procedure  puts  the  designer 
out  of  a  controlling  position  while  putting  total  reliance 
on  an  optimization  criteria.   Optimization  criteria  are 
difficult  to  define  if  more  than  just  the  cost  of  the 
vessel  is  in  question.  It  is  felt  that  such  questions 
as  the  balance  between  payload  and  cost  are  best  left 
to  the  designer  and  not  to  a  rigid  formula. 
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APPENDIX  A  f 

Description  of  Input 


Input  for  the  cutter  model  can  be  divided  into  two 
main  groups.   The  first  group  contains  the  data  for  the 
Cr  arrays  used  in  estimating  propulsion  power.   The  second 
group  contains  the  data  for  individual  ships  that  are  run. 

Data  for  the  C  arrays  is  read  only  once  regardless 
of  how  many  test  cases  are  run.   This  data  precedes  all 
other  data  and  the  values  never  change.   The  Cr  array 
data  can  be  stored  on  tape  or  disk  if  desired.   A  list 
of  the  values  used  is  given  in  Appendix  D. 

The  second  group  of  data  allows  the  user  to  specify 
the  payload  and  performance  as  well  as  various  other 
parameters  of  the  ship  the  computer  is  to  design.   In 
order  for  the  computer  to  produce  reliable  results,  the 
input  data  must  be  the  best  that  it  is  possible  to  specify, 
Past  designs  are  a  good  starting  point  for  determining 
reasonable  input  values. 

The  minimum  number  of  data  cards  for  the  second 
group  of  input  is  14.   These  cards  are  described  below. 
Where  more  than  one  card  may  be  used  this  fact  is  noted. 

Card  1  —  A  single  number  is  entered  on  this  card 

in  card  column  1.   The  values  that  can  be  used 
are : 

0  -  the  program  will  stop 

1  -  this  is  the  normal  value  to  use 
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2  -  a  fine  hull  form  will  be  assumed  for 

the  first  iteration,  Cv  *  .0025. 


Card  2  —  This  card  contains  a  heading  that  is 
printed  above  the  output.   All  80  columns 
can  be  used  for  any  alphanumeric  description 
desired. 

Card  3  —  A  single  value  is  entered  in  card  column  2, 
The  total  number  of  armament,  aircraft  and 
cargo  payload  data  cards  should  be  entered, 
A  maximum  of  20  can  be  used. 

Card  k   —  This  is  really  a  set  of  cards  containing 
the  armament,  aircraft  and  cargo  payload  data. 
The  number  of  cards  used  is  specified  on  card  3» 
Each  of  these  cards  has  the  following  format » 
(All  data  should  contain  a  decimal  point) 
columns  1--20  Any  alphanumeric  title 

describing  the  remainder  of  the  data 
on  the  card, 
columns  21-25  Cargo  weight,  enter  the 

weight  in  tons, 
columns  26-30  Cargo  vcg,  enter  the  ratio 
of  cargo  vcg  divided  by  the  depth  of 
the  ship  at  amidships, 
columns  31-35  Group  7  weight, ,enter 
weight  in  tons. 
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columns  36-^0  Group  7  vcg,  enter  ratio 

of  kg  to  depth  at  amidships, 
columns  41-^5  Ammunition  weight,  enter 

weight  in  tons. . 
columns  46-50  Ammunition  vcg,  enter  ratio 

of  kg  to  depth  at  amidships, 
columns  jl-55     Aircraft  weight,  enter 

weight  in  tons, 
columns  56-60  Aircraft  vcg,  enter  ratio 

of  kg  to  depth  at  amidships. 
columns  61-65  Electrical  load,  enter  value 

in  KW. 
columns   66-70     Deckhouse   Area,   enter  area 

required  to  be  in  deckhouse  in  sq.  ft. 
columns  71-75  Hull  Area,  enter  all  remaining 

required  area  in  sq.  ft. 
columns  76-8O  Group  7  cost,  enter  acqui- 
sition cost  in  thousands  of  dollars. 
Except  for  the  moment  arms,  the  values  given 
on  these  cards  are  summed  by  category  and  only  this 
sum  is  used  in  calculations.   Therefore  it  does  not 
matter  whether  all  the  data  is  given  on  one  card  or 
on  several.   The  only  requirement  is  that  the  vcg 
values  apply  to  the  weights  which  immediately  precede 
them.   All  values  which  are  included  in  light  ship 
weight  group  7  must  be  given  as  data  on  these  cards. 
If  there  is  no  payload  data,  card  3  should 
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contain  a  1  and  a  blank  card  should  be  included 

for  card  4. 
Card  5  —  A  single  value  is  entered  in  card  column  2. 
This  is  a  count  of  the  number  of  data  cards 
described  under  card  6.  A  maximum  of  20  can 
be  used. 

Card  6  —  This  is  really  a  set  of  electronics  input 
data  cards  up  to  a  maximum  of  20.   The  number 
of  cards  is  specified  by  card  5-  2ach  of  these 
cards  has  the  following  format:  (All  data  should 
contain  a  decimal  point) 

columns  1 — 20  Any  alphanumeric  title 
describing  the  remainder  of  the 
data  on  the  card, 
columns  21-25  Group  4  weight,  enter 

weight  in  tons, 
columns  26-30  Group  4  vcg,  enter  ratio 

of  kg  to  depth  at  amidships, 
columns  31-35  Electrical  load  in  KW. 
columns  36-*+0  Deckhouse  area,  enter  area 

required  to  be  in  deckhouse  in  sq.  ft. 
columns  41-^5  Hull  area,  enter  remaining 

required  area  in  sq.   ft. 
columns  46-50     Group  4  cost,   enter  acqui- 
sition cost  in  thousands   of  dollars. 
Again  the  data  may  be   broken  up  in  any 
way  desired  as  long  as   the   correct  sets   of 
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weight  and  vcg  are  given*   On  all  of  these 

cards  and  .on  the  cards  .described  under  card  k, 
-,arci  the -area  given* under- deckhouse  .area  should 

included  only  the  area .that .is  required  to  be 
in  the  deckhouse.   If  an  area  can  be  either 
in  the  deckhouse  or  in  the  hull,  it  should  be 
listed  under  hull  area. 

All  items  included  in  light  ship  weight 
group  b   must  be  entered  on  these  cards.   This 
includes  standard  navigation  equipment.   If 
there  is  no  electronics  data,  card  5  should 
contain  a  1  and  a  blank  card  should  be  included 
»r  :  for -card  6. 
Card  7  — -  The  number  of  officers,  CPOs  and  crew 

members  are  specified,  on  this  card.   All  numbers 
must'  be 'right- justified  integers. 
'  a- jj columns!' — 5  number  of  officers 
car.'  columns  6--10  number,  of  chief  petty  officers 
•i"'  'columns  11-15  number  of  enlisted  men 
Card  8  —  This ■ card"  contains  cost  indices  input. 
All  values  should  contain  a  decimal  point. 

'  j.i  ■-    i:   -columns  1 5  labor  rate  in  dollars/hour 

«??  : -r<columns  6 — 10  group  1  index  - 

.-•••"    columns   11-15     group  2  index 

*    -  •  columns  16-20     group  3  index 

...«.:  columns  21-25     group  5  index 

"-columns  26-30     group  6  index 
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The  values  of  the  labor  rate  and  cost  indices 

may  be  found  under  Section  3.1^,  Figures    and 

Card  9  —  The  number  of  endurance  days  for  dry 

provisions  and  aircraft  fuel  payload  are  given 

on  this  card.   All  values  should  contain  a 

decimal  point, 

columns  \ 5  endurance  days  for  dry 

provisions 

columns  6 — 10  weight  of  aircraft  fuel,  tons 

columns  11-15  vcg  of  aircraft  fuel, 

enter  the  ratio  of  kg  to  depth  at 

amidships. 

Card  10  —  A  single  value  is  entered  in  card  column 

2.   Two  options  are  available.   If  a  1  is  entered, 

use  card  11a.   This  option  allows  the  user 

to  input  sustained  speed  and  have  the  computer 

calculate  horsepower.   If  a  2  is  entered  use 

card  11  b.  Here  particulars  of  the  machinery 

plant  are  input  and  the  sustained  speed  is 

calculated.   Any  other  value  will  have  an 

unknown  effect. 

Card  11a  —  Use  this  card  type  if  a  value  of  1  is 

entered  on  card  10.   The  maximum  sustained 

speed, .endurance  speed  and  range  and  propulsive 

coefficients  at  endurance  speed  and  maximum 

speed  are  entered  on  this  card.   All  data  should 

contain  a  decimal  point. 


213 

columns  1 5  sustained  speed  in  knots 

columns  6 — 10  endurance  speed  in  knots 

columns  11-20  endurance  range,  nautical  miles 

columns  21-25  propulsive  coef •  at  endurance 

speed 

columns  26-30  propulsive  coef.  at  maximum 

sp*3d. 

Card  lib  —  Use  this  card  type  if  a  value  of  2  is 

entered  on  card  10.   The  particulars  of  a  specific 
machinery  plant  are  entered  on  this  card.   All 
data  should  contain  a  decimal  point. 

columns  1 5  maximum  sustained  horsepower, 

enter  as  HP/1000 
columns  6 — 10  minimum  machinery  box 

depth  in  feet 
columns  11-15  machinery  box  length  in  feet 
columns  16-20  total  weight  of  group  2,  tons 
columns  21-25  specific  fuel  consumption 

at  max.  horsepower  in  lbs/HP-hr 
columns  26-30  specific  fuel  consumption 

at  half  power  in  lbs/HP-hr 
columns  31-35  endurance  speed  in  knots 
columns  36-^5  endurance  range,  nautical  miles 
columns  46-50  propulsion  auxiliaries 

electrical  load  in  KW 
columns  51-55  propulsive  coef.  at  endurance 

speed 
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columns  56-60  propulsive  coef .  at  maximum 

speed 

columns   61-65  acquisition  cost  of  group  2, 

thousands  of  dollars 

columns  66-70  lub  oil  weight  in  tons. 

Card  12  —  This  card  contains  various  input  options 
and  dimensions.   All  values  except  the  first 
should  contain  a  decimal  point. 

columns  1 5  MTYPE,  a  right- justified 

integer  is  entered  here. 

If  MTYPE  =  0  the  program  returns 

to  read  a  new  card  1 
If  MTYPE  =1   a  diesel  plant  is  assumed 
if  the  machinery  particulars  have 
not  been  given 
If  MTYPE  =  2   a  CODOG  plant  is  assumed 
If  card  10  contains  a  value  of  2,  any 
value  of  MTYPE  except  0  may  be  used, 
columns  6 — 10  free  surface  correction, 
enter  value  in  feet.   This  doesnot 
have  to  be  a  positive  number.   A 
negative  number  can  be  interpreted 
as  a  KG  below  the  average  built  into 
the  program . 
columns  11-15  length  between  perpendiculars, 
enter  the  ships  length  in  feet. 
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columns  16-20  C  ,  enter  a  value  for 

prismatic  coefficient 

columns  21-25  Cx,  enter  a  value  for  mid- 
ship section  coefficient 

columns   26-30     GI<l/£eam,    a  value    of   0.1 
is  commonly  used 

columns   31-35     design  and  builders  margin, 
enter  a  value   as  a  decimal,   i.e.    a 
10#  margin   is   entered  as   0.1 

columns  36-^5  AOf,   this  is  the  frictional 
resistance  correction  used  in  horse- 
power calculations.   A  value  of  0.000** 
is  commonly  used. 

Card  12  may  "be  repeated  as  many  times  as 
desired  with  new  data  each  time.  Each  new  card 
is  interpreted  as  a  new  data  case. 

Card  13  —  If  it  is  desired  to  input  an  entirely 
new  set  of  payload  and  performance  data  or  to 
stop,  a  value  of  0  should  be  given  in  column 
5  of  card  12.  This  will  cause  a  new  card  1  to 
be  read. 

Card  I**  —  A  new  card  1  must  be  included.  If  new 
data  will  be  given  all  of  the  data  cards  must 
be  repeated.  A  value  of  0  on  card  1  will  cause 

the  program  to  stop. 
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A  sample  listing  of  the  group  2  data  cards  is  given 

below.   The  data  used  refers  to  a  210*  WMBC. 
Card  1 
1 
Card  2 

210  WKEC  TEST  RUN 
Card  3 

5 

Card  4 

Group  700  4.7  1.64  385.  76. 

Group  701  .6  1.01 

Group  702  2.8   .60 
Group  750  .4   .31 

Anmunition  2.56  .41 

columns     35   40   45   50  75   80 

Card  5 

6 

Card  6 

Group  400       .7  1.8521.99  319.      190. 

Grouo  401      2.4  1.3 

Group  402       .2  I.56 

Group  403  10.4   .92 

Group  404      9.1  2.69 

Group  450      3.3  I.76 

1     •     f     f         t 

columns   25   30   35   40       50 
Card  7 

9    8   57 
•    1    t 

5  10  15  -  columns 
Card  3 

4.75  1.4  1.6  1.19  1.55  1.39 

5  10  15   20   25   30  -  columns 


Card  9  217 

60.   0,   0. 
t     i     i 

5   10   15  -  columns 

Card  10 

2 

Card  lib 

4.5  19.  32.102.6  .'?z     .44  14.  4000.  .71  .55  .6  599.  3.85 
t    t    t    t    i    *    i      tiift    t    t 

5   10   15   20   25   30   35    45  50  55     60  65       70 

Card  12 

1-1.61   200.      .59    .892        .1        0.  .0004 

tftitii  • 

5   10   15   20   25   30   35      45  -  columns 
Card  13 

0  in  column  5 
Card  14 
0  in  column  1 


APPENDIX  B 
SAMPLE  OUTPUT 
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NAYSHIPS  HULL  GROUP  WEIGHT  CLASSIFICATION  OP  I965 

Group  1  -  Hull  Structure 


100  -  Sholl  Plating  and  Planking 

101  -  Longitudinal  and  Transverse  Framing 

102  -  Inner  Bottom 

103  -  Platforms  and  Flats  Below  Lowermost  Continuous  Deck 

104  -  Fourth  and  Lower  Continuous  Decks 

105  -  Third  Deck 

106  -  Second  Deck 

107  -  Main  Deck  and  Hangar  Deck 

108  -  Forecastle  and  Poop  Docks 

109  -  Gallery  Deck 

110  -  Flight  Deck,  Landing  Platforms,  Special  Purpose 

Decks  above  Weather  Deck 

111  -  Superstructure 

112  -  Foundations  for  Propulsion  Plant  Machinery 

113  -  Foundations  for  Auxiliaries  and  Other  Equipment 

114  -  Structural  Bulkheads 

115  -  Trunks  and  Enclosures 

116  -  Structural  Sponsons 

117  -  Armor 

118  -  Aircraft  Fuel  Saddle  Tank  Structure 

119  -  Structural  Castings,  Forgings  and  Equivalent  Weldments 

120  -  Sea  Chests 

121  -  Ballast  and  Buoyancy  Units,  Fixed  or  Fluid 

122  -  Doors  and  Closures,  Special  Purpose 

123  -  Doors,  Hatches,  Manholes  and  Scuttles  -  Nonballistic 

124  -  Doors,  Hatches,  Manholes  and  Scuttles  -  Ballistic 

125  -  Kingposts  and  Support  Frames 

127  -  Sonar  Dome 

128  -  Masts,  Towers,  Tetrapods,  and  Service  Platforms 

150  -  Welding,  Riveting,  and  Fastenings 

151  -  Free  Flooding  Liquids 
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Group  2  -   Propulsion 


200  -  Boilers  and  Energy  Converters  (Non-Nuclear) 

201  -  Propulsion  Units 

202  -  Main  Condensers  and  Air  Ejectors 

203  -  Shafting,  Bearings,  and  Propellers 

204  -  Combustion  Air  Supply  System 

205  -  Uptakes  (Smoke  Pipes) 

206  -  Propulsion  Control  Equipment  (Non-Nuclear) 

207  -  Main  Steam  System 

208  -  Feedwater  and  Condensate  System 

209  -  Circulating  and  Cooling  Water  Systems 

210  -  Fuel  Oil  Service  System 

211  -  Lubricating  Oil  System 

212  -  Nuclear  Steam  Generators 

213  -  Reactors 

214  -  Reactor  Coolant  System 

215  -  Reactor  Coolant  Service  Systems 

21 6  -  Reactor  Plant  Auxiliary  Systems 

217  -  Nuclear  Power  Control  and  Instrumentation 

218  -  Radiation  Shielding  (Primary) 

219  -  Radiation  Shielding  (Secondary) 

250  -  Propulsion  Repair  Parts 

251  -  Propulsion  Operating  Fluids 


Group  3  -  Electric  Plant 


300  -  Electric  Power  Generation 

301  -  Power  Distribution  Switchboards 

302  -  Power  Distribution  System  (Cable) 

303  -  Lighting  System  (Distribution  and  Fixtures) 

350  -  Electric  Plant  Repair  Parts 

351  -  Electric  Power  Generator  Fluids 
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Group  4  -  Communication  and  Control 


400  -  Navigation  Equipment   (Non-Electronic) 

401  -  Interior  Communication  Systems  and  Equipment 

402  -  Gun  Fire  Control  System 

403  -  Countermeasure  Systems   (Non -Electronic) 

404  -  Electronic  Countermeasure  Systems   (ECM) 

405  -  Missile   Fire  Con^ol  Systems 

406  -  ASY/,   Torpedo  Fire  Control  Systems   (Surface  Ships) 

407  -  Torpedo  Fire  Control  System  (Submarines) 

408  -  Radar  Systems 

409  -  Radio  Communication  Systems 

410  -  Electronic  Navigation  Systems 

411  -  Space  Vehicle  Electronic  Tracking  Systems 

412  -  Sonar  Systems 

413  -  Electronic  Tactical  Data  System 

450  -  Communication  and  Control  Repair  Parts 

451  -  Communication  and  Control   Operating  Fluids 
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Group  5  -  Auxiliary  Systems 


500  -  Heating  System 

501  -  Ventilation  System 

502  -  Air  Conditioning  System 

503  -  Refrigerating  Spaces,  Plant,  and  Equipment 

504  -  Gasoline,  HEAP,  All  Liquid  Cargo  Piping,  Oxygen* 

Nitrogen,  and  Aviation  Lubricating  Oil  Systems 

505  -  Plumbing  Installations 

506  -  Firemain,  Flushi^,  Sprinkler,  Washdoisn,  and  Salt 

Water  Service  Systems 

507  -  Fire  Extinguishing  Systems 

508  -  Drainage,  Ballast,  Trimming,  Heeling,  and  Stabiliser 

Tank  Systems 

509  -  Fresh  Water  System 

510  -  Scuppers  and  Back  Drains 

511  -  Fuel* and  Diesel  Oil  Filling,  Venting,  Stowage  and 

Transfer  Systems 

512  -  Tank  Heating  Systems 

513  -  Compressed  Air  System 

5l£  -  Auxiliary  Steam,  Exhaust  Steam,  and  Steam  Drains 

515  -  Buoyancy  Control  System 

516  -  Miscellaneous  Piping  Systems 

517  -  Distilling  Plant 

518  -  Steering  Systems 

519  -  Rudders 

520  -  Mooring,  Towing,  Anchor  and  Aircraft  Handling  Systems 

and  Deck  Machinery 

521  -  Elevators,  lioving  Stairways,  Stores  Strikedown,  and 

Stores  Handling  Equipment 

522  -  Operating  Gear  for  Retracting  and  Elevating  Units 

523  -  Aircraft  Elevators 

524  -  Aircraft  Arresting  Gear,  Barriers,  and  Barrioades 

525  -  Catapults  and  Jet  Blast  Deflectors 

527  -  Diving  Planes  and  Stabilizing  Fins 

528  -  Replenishment  at  Sea  and  Cargo  Handling 

550  -  Auxiliary  Systems  Repair  Parts 

551  -  Auxiliary  Systems  Operating  Fluids,  Gases,  and 

Stabilizer  Fluids 
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Group  6  -  Outfit  and  Furnishings 

600  -  Hull  Fittings 

601  -  Boats,  Boat  Stowage,  and  Handling 

602  -  Rigging  and  Canvas 

603  -  Ladders  and  Gratings 

6oh   -  Nonstructural  Bulkheads  and  Nonstructural  Doors 

605  -  Painting 

606  -  Deck  Covering 

607  -  Hull  Insulation 

608  -  Storerooms,  Stowages  and  Lockers 

609  -  Equipment  for  Utility  Spaces 

610  -  Equipment  for  Workshops,  Laboratories,  and  Test  Areas 

611  -  Equipment  for  Galloy,  Pantry,  Scullery  and  Commissary 

612  -  Furnishings  for  Living  Spaces 

613  -  Furnishingst  Offices,  Control  Centers,  Machinery  Spaces 
61^  -  Furnishings  for  Medical,  Dental,  and  Pharmaceutical 
615  -  Radiation  Shielding  for  Nuclear  Support  Facilities 

650  -  Outfit  and  Furnishings  Repair  Parts 

651  -  Outfit  and  Furnishings  Operating  Fluids 


Group  7  -  Armament 

700  -  Guns  and  Gun  Mounts 

701  -  Ammunition  Handling  Systems 

702  -  Ammunition  Stowage 

703  -  Special  Weapons,  Handling  and  Stowage 

?oh   -  Rocket  and  Missile  Launching  Devices  (Surface  to  Air, 
Surface  to  Surface ,  and  Sub-Surface ) 

705  -  Rocket  and  Missile  Launching  Devices  (Antisubmarine 

Warfare ) 

706  -  Rocket,   Missile,   and  Components  Handling  Systems 

707  -  Rocket,   Missile,   and  Components  Stowage 

708  -  Torpedo  Tubes 

709  -  Torpedo  Handling  and  Stowage 

710  -  Mine  Handling  Systems  and  Stowage 

711  -  Small  Arms  and  Pyrotechnic  Stowage 

712  -  Air-Launched  V/eapons  Handling  Systems 

713  -  Air-Launched  Weapons  Stowage 
720  -  Cargo  Munition  Stowage 

7^5  -  General  Arrangement  -  Armament  Drawings 

750  -  Armament  Repair  Parts 

751  -  Armament  Operating  Fluids 
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